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Dedication 


A Marialuisa 



Preface 


This book is the result and the consequence of a series of seminars and lectures 
I gave on the subject of anyons in many places. Everything started in September 
1990 at the meeting “Condensed Matter and High Energy Physics (Fractional 
Statistics)” in Chia (Cagliari, Italy) organized by Professors S. Fubini and L. 
Alvarez-Gaume, who invited me to speak. Then, in the spring of 1991, I had 
the pleasure of giving a series of lectures on anyons at the Centro Brasileiro de 
Pesquisas Fisicas in Rio de Janeiro (Brazil) at the invitation of Professors J.A. 
Mignaco and M.A. do Rego Monteiro. Finally I had the opportunity of presenting 
part of the material of this book in a series of seminars at the Center for Theoretical 
Physics of M.I.T. (Cambridge, USA) and at the Institute for Theoretical Physics 
of S.U.N.Y. at Stony Brook (New York, USA). All of these, together with the 
many talks I gave on the subject, have been extremely valuable and challenging 
experiences for me. In this book I have tried to collect together and include as 
much as I have been able of all the suggestions, comments and observations that 
were made during my lectures, and also I have tried to answer the questions that 
were posed to me on those occasions. 

Obviously I do not presume to have done a perfect job; I am the first to 
realize that my book does not cover the whole subject and that many interesting 
topics have been left out from my exposition. This is due partly to lack of time, 
and partly also to my lack of familiarity with those topics. However, as a partial 
justification, I would like to mention that many excellent reviews on these subjects 
are available, and hence a further duplication would have been superfluous. Among 
these reviews, to which I refer the interested reader, I would like to mention: 

- D.P. Arovas, in Geometric Phases in Physics edited by A. Shapere and F. 
Wilczek (World Scientific, Singapore 1989); 

- A.P. Balachandran, E. Ercolessi, G. Morandi and A.M. Srivastava, Hubbard 
Model and Anyon Superconductivity (World Scientific, Singapore 1990); 

- G.S. Canright and S.M. Girvin, Science 247 (1990) 1197; 

- S. Forte, Rev. Mod. Phys. 64 (1992) 193; 

- P. de Sousa Gerbert, Int. Jour. Mod. Phys. A6 (1991) 173; 

- G.A. Goldin and D. H. Sharp, Int. Jour. Mod. Phys. B5 (1991) 2625; 

- R. Iengo and K. Lechner, Phys. Rep. 213 (1992) 179; 

- J.D. Lykken, J. Sonnenschein and N. Weiss, Int. Jour. Mod. Phys. A6 (1991) 
5155; 
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M. Stone, Int. Jour. Mod. Phys. B4 (1990) 1465; 

X.G. Wen and A. Zee, Santa Barbara preprint, NSF-ITP-89-155 (November 
1989) 

F. Wilczek, in Fractional Statistics and Anyon Superconductivity edited by F. 
Wilczek (World Scientific, Singapore 1989); 

F. Wilczek, Int. Jour. Mod. Phys. B5 (1991) 1273. 
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1. Introduction 


Since the early days of quantum mechanics, it has been well-known that the phys¬ 
ical behavior of collections of identical particles is influenced not only by conven¬ 
tional interactions but also by the statistics of the particles. These are classified 
into bosons and fermions depending on whether they obey the Bose-Einstein or the 
Fermi-Dirac statistics respectively. The particle statistics determines the structure 
of the many-body wavefunctions, that turn out to be symmetric under permu¬ 
tations of identical bosons, but antisymmetric under permutations of identical 
fermions. 

The distinction between symmetric and antisymmetric wavefunctions, or equiv¬ 
alently between bosons and fermions, is as old as quantum mechanics, and has 
many important and far reaching consequences. For example, the Pauli exclusion 
principle according to which two identical fermions cannot occupy the same quan¬ 
tum state, naturally follows from the antisymmetry of the fermionic wavefunctions. 
By contrast, the probability for bosons to be in some state is enhanced when that 
state is already occupied, so that the ground state of a system of identical bosons 
is realized when all particles condense into the lowest energy state. 

For many years till quite recently, bosons and fermions have been thought to 
be the only logical possibilities. This is certainly correct for particles moving in at 
least three space dimensions, but for systems in two dimensions the situation is 
more interesting. Here the quantum statistics is not limited to the Bose-Einstein 
and the Fermi-Dirac cases, but rather a continuous interpolation between these two 
extremes is possible. Therefore, in two space dimensions we do not have only bosons 
and fermions, but also particles with any statistics in between. These particles are 
called anyons and are the subject of this book. 

According to the spin-statistics connection, bosons are particles with integer 
spin while fermions are particles with half-integer spin. Extending this idea, one 
may conjecture that anyons are characterized by a fractional spin, or more gener¬ 
ally by fractional quantum numbers. As we will see hereinafter, this is indeed the 
case. However, the appearance of unusual quantum numbers is not an entirely new 
feature of an anyonic theory. In fact, long ago T.H.R. Skyrme showed that in the 
non-linear sigma model of pions - particles obeying Bose statistics - there are soli- 
ton configurations which carry the quantum numbers of the nucleons - particles 
obeying Fermi statistics (Skyrme 1961). Actually these solitons are true fermions 
(Finkelstein and Rubistein 1968), despite the fact that they arise from bosonic 
fields. Similarly, magnetic monopoles in three space dimensions can be fermions 
even when the elementary fields are bosons (Hasenfratz and’t Hooft 1976), and 
can carry fractional (Jackiw and Rebbi 1976) or even irrational electric charges 
(Witten 1979). 
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Introduction 


These examples show that one can trasmute the statistics and form composite 
objects which obey different rules from their constituents. In three dimensions, 
however, there is not much choice since only bosons and fermions are possible. On 
the contrary in two dimensions, as we have mentioned, there is a whole variety 
of different possibilities, and thus one can form objects obeying exhotic fractional 
statistics. It was only quite recently that actually these new possibilities for spin 
and statistics were understood. In the late 70’s, J.M. Leinaas and J. Myrheim 
clearly stated the principles involved in the concept of fractional statistics (Leinaas 
and Myrheim 1977), and identified its origin in the peculiar topological properties 
of the configuration space of collections of identical particles. In fact, this space 
is multiply connected in two dimensions but only doubly connected in three or 
more dimensions. Furthermore, in two dimensions the fundamental group of this 
configuration space is the braid group (Artin 1926, 1947) which replaces the more 
familiar permutation group when discussing statistics. A few years later, G.A. 
Goldin, R. Menikoff and D.H. Sharp reached similar conclusions using a completely 
different method based on the rigorous study of the unitary representations of 
current algebras and diffeomorphism groups (Goldin et al. 1980, 1981) l . 

Most of the great interest that anyons have attracted in the past few years 
derives from the (unexpected) applications of these ideas to certain two-dimesional 
condensed matter systems, most notably those exhibiting the fractional quantum 
Hall effect (see for instance (Prange and Girvin 1990)). In this case a series of 
new states of matter emerge as incompressible quantum liquids (Laughlin 1983) 
around which the low-energy excitations are localized quasi-particles with unusual 
fractional quantum numbers, i.e. anyons. Furthermore, it is also very likely that 
anyonic excitations with fractional statistics exist in films of liquid 3 He in the 
A-phase (Volovik and Yakovenko 1989). The application of anyons to the theory 
of high temperature superconductivity has also been considered quite extensively 
(for reviews see (Wilczek 1990; Lykken et al. 1991)), but their actual relevance in 
this context is quite controversial and doubtful. 

In this book I will consider the subject of anyons from many different points of 
view. In Chapter 2 I will give a general introduction to fractional statistics in two 
dimensions and discuss its relation to the braid group. In Chapter 3 I will show that 
that anyonic statistics can be implemented on ordinary bosons or fermions with 
the addition of special topological interactions, and in particular I will consider 
the connection between anyons and Chern-Simons field theories. In Chapter 4 
I will show how to describe anyons with multi-valued wavefunctions carrying a 
representation of the braid group. Chapter 5 will be devoted to a discussion of 
the second-quantized theory of anyons and to the spin-statistic connection. In 
Chapter 6 I will consider the case of anyons in a magnetic field, which is relevant 
for the applications to the fractional quantum Hall effect, and derive the energy 
spectrum and the wavefunctions. In Chapter 7 I will address the important problem 
of the statistical mechanics of systems of many anyons and discuss their partition 
functions, their virial coefficients and their magnetic moments. Finally, in Chapters 

1 In this book we will not discuss this approach and we refer the reader to the original 
literature as well as to the review article (Goldin and Sharp 1991). 
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8 and 9 I will consider the relevance of anyons to the theory of the fractional 
quantum Hall effect and comment on the surprising relation between anyons and 
conformal field theory. 

My exposition certainly does not cover the whole subject of anyons; in par¬ 
ticular I will not consider at all the mechanisms of anyon superconductivity, the 
non-abelian extensions of fractional statistics and the gravitational anyons. More¬ 
over, the topics I treated in Chapters 7 and 9 are still the subject of active current 
investigations, and thus my presentation is not conclusive. The list of references 
I have included, contains the bibliography which is directly relevant to my expo¬ 
sition; though quite extended, it is by no means complete and exhaustive, and I 
sincerely apologize for all omissions, which however are completely unintentional. 



2. Introduction to Fractional Statistics 

in Two Dimensions 


So many gods, so many creeds, 

So many paths that wind and wind, 


Ella Wheeler Wilcox 


In two space dimensions there are possibilities for spin and statistics which do 
not occur in higher dimensions: the spin is not quantized in integer or half-integer 
units, and particles are not necessarily bosons or fermions. Indeed they may obey 
any statistics and thus are called anyons (Wilczek 1982a,b, 1990). These new pos¬ 
sibilities for spin and statistics arise from the peculiar topological properties of the 
configuration space of collections of identical two-dimensional particles (Leinaas 
and Myrheim 1977; Goldin et al 1980, 1981). However anyons are not simply topo¬ 
logical fantasies or objects of purely mathematical interest; on the contrary they 
might play an important role in certain physical phenomena of the real world. Of 
course, since we are living in at least three space dimensions where particles can be 
only bosons or fermions, anyons are not real particles. However there exist certain 
condensed-matter systems (for example thin layers at the interface between dif¬ 
ferent semiconductors) that can be regarded effectively as two-dimensional. Their 
localized excitations (if they exist) are quasi-particles subject to the rules of a 
two-dimensional world. It is these quasi-particles that may be anyons and may be 
observed in certain cases. For example the collective excitations above the ground 
state of systems exhibiting the fractional quantum Hall effect (for a review see 
(Prange and Girvin 1990)) have been identified as localized quasi-particles of frac¬ 
tional charge (Laughlin 1983), fractional spin and fractional statistics (Arovas et 
al. 1984; Halperin 1984), and thus they can be regarded as anyons. Furthermore, 
anyons are conjectured to play a role also in the theory of high temperature su¬ 
perconductivity (Chen et al. 1989), even though in this case no conclusive word 
can be said at the moment (Lyons et al. 1990; Kiefl et al. 1990; Spielman et al. 
1990). 

The fact that in two dimensions the spin may be arbitrary (that is not only 
integer or half-integer in units of ft) is not too surprising. In fact, the rotation 
group in d = 2 is 50(2), which is an abelian group. Therefore, there are no 
commutation relations to quantize, and no restrictions on the possible eigenvalues 
of the spin come from the algebra. If we believe in some connection between spin 
and statistics, the arbitrariness of the spin leads us to conceive the idea that also 
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the statistics may be arbitrary in two dimensions 2 . As we will see, this is indeed 
the case. 

The notion of statistics is usually related to the sign that a many-body wave- 
function acquires when any two particles are interchanged. If the wavefunction 
gets a plus sign, we say that it describes a bosonic system; if it gets a minus sign 
we say that it describes a fermionic system. Actually, we can give a more gen¬ 
eral definition of statistics: let ^(1,2) be the wavefunction describing two identical 
hard-core particles with definite angular momentum, and let us assume that when 
we move particle 2 around particle 1 by an azimuthal angle A(p (see Fig. 2.1), the 
wavefunction changes according to 

^( 1 , 2 )-^'( 1 , 2 ) = ^"^( 1 , 2 ) . ( 2 . 1 ) 





2 





V>'(1,2) 


Fig. 2.1. When particle 2 is moved around particle 1 by an angle A<p, the wavefunction 
^ changes into tj/. 

The phase acquired by the wavefunction depends on a parameter v which is 
usually called statistics . The meaning of v and of (2.1) becomes more clear if we 
consider the exchange of the two particles. This can be realized in two ways: 

i) Moving particle 2 around particle 1 by an angle A(p = 7r, and then performing 
a rigid translation of the center of mass to reach the initial spatial configuration 
(see Fig. 2.2); 

it) Moving particle 2 around particle 1 by an angle A<p = — 7r and then performing 
the rigid translation of the center of mass (see Fig. 2.3). 

In the first case the wave function acquires a phase exp(i7ri/) according to (2.1), 
whilst in the second case, it gets a phase exp(— mu). This simple example shows 
that there is a dramatic difference between two and three or more dimensions as 
far as statistics is concerned. In fact in d > 3 there is no intrinsic difference between 
case i) and case ii) since, as is clear from the figures, we can always deform the 
transformation in i) into the one in it) in a continuous way: for example we can 
lift the path in i) into the third dimension, then fold it back onto the plane, and 
finally superpose it to the path in ii). That we cannot distinguish between i) and 
ii) is a reflection of the fact that the azimuthal angle A(p is not a well-defined 

2 See Chapters 3 and 5 for a discussion on the spin-statistics connection for anyons. 
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2 



Fig, 2,2. The exchange of two particles is realized by moving particle 2 counterclockwise 
around particle 1 by an angle A(p = it. 



Fig, 2.3. The exchange of two particles is realized by moving particle 2 clockwise around 
particle 1 by an angle A(p = —it. 

quantity in d > 3. The consequences of this observation are very important: if t) 
and it) are the same physical operation, the wave function should behave in the 
same way, which means that in d > 3 one must have 

e i7ri/ = e -i7ri/ . (2.2) 

Clearly (2.2) can only be true if v = 0,1 (modulo 2). This simple and intuitive 
example shows that in d > 3 the statistics cannot be arbitrary. Under the exchange 
of two particles, the wavefunction picks up either a plus sign if v = 0 (bosonic 
statistics) or a minus sign if v — 1 (fermionic statistics). There are no other 
possibilities. 

The situation changes drastically in two dimensions where it is not any more 
possible to deform continuously the path in i) into the one in it) since by assump¬ 
tion the particles cannot go through each other. Hence in d = 2, i) and it) are 
two topologically and physically distinct operations - a reflection of the fact that 
the azimuthal angle Aip is a well-defined quantity in two dimensions. The equality 
(2.2) does not necessarily hold any more and the statistical parameter v can be 
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arbitrary, at least in principle 3 . From this example we learn another important 
fact: in d = 2 it is not enough to specify the initial and final configurations to 
completely characterize a system; it is also necessary to specify how the different 
trajectories wind or braid around each other. In other words the time-evolution of 
the particles is important and cannot be neglected in d = 2. As we will see later, 
this fact implies that in order to classify and characterize anyons, the representa¬ 
tions of the permutation group must be replaced by those of the more complicated 
braid group. 

One of the essential features of anyons is the violation of the discrete symmetries 
of parity P and time reversal T. This is particularly evident in our two-particle 
example. In fact, under a parity or time reversal transformation, the phase acquired 
by the wave function ^(1,2) changes according to 


—> e =Fil ' ir 



for cases i) and it) respectively, and for i/ / 0,1 a violation of P and T occurs. 
This breaking of the discrete symmetries is the hallmark of anyonic statistics and 
is the signal that experiments look for to demonstrate the presence (or absence) 
of anyons in two-dimensional phenomena (Halperin et al. 1989). 

Now we try to be more rigorous and study the problem of anyonic statistics 
from a more formal point of view, following essentially the presentation of Y.-S. Wu 
(Wu 1984a,b). Let M $ be the configuration space of a collection of N identical 
hard-core particles in d dimensions and let q and q f be two arbitrary points in 
Mf;. According to the standard path-integral formulation of quantum-mechanics 
(see for instance (Feyman and Hibbs 1965; Laidlaw and Morette-de Witt 1971)), 
the amplitude for the system to evolve from the configuration q at time t to the 
configuration q ' at time t f is given by the kernel 


1 = W,t'\<I.‘) = [ Vq AS:'* 

Jg(t)=g;g(t / )=g / 


(2.4) 


where C(q , q) is the lagrangian density for the iV-particle system and the symbol 
fq(t)= q ‘ q ( t ')=q' denotes the sum over all paths connecting q at time t to q ' at time 

t!. The kernel K(q\t g,t) evolves the single-valued wavefunction ^(q^t) according 
to 



Without any loss of generality, we choose q = q f and hence describe loops in 
Mtf. Two loops are considered equivalent (or homotopic) if one can be obtained 
from the other by a continuous deformation. All homotopic loops are grouped 


3 We will see that actually there are restrictions on v coming from the topology of the 
two dimensional space. For example for particles moving on a torus, v can be only a 
rational number. 



9 


into one class and the set of all such classes is called the fundamental group 4 
and is denoted by 7Ti. Thus an element of 7Ti ( M $) is simply the set of all loops 
in Mfj which can be continuously deformed into each other. On the other hand, 
loops belonging to two different elements of 7Ti (Af$) cannot be connected by a 
continuous transformation. With this in mind we can organize the sum over all 
loops in (2.4) into a sum over homotopic classes a € 7Ti ( Mfj ) and into a path- 
integral in each class. Therefore (2.4) may be rewritten as 




^2 K a (q,t';q,t) 

qc€*-i(m£) 

E / 

a e*i(Mi l y 9 ° (t)=9,9a( * )=9 


e ji J* dT ^llc(r),q a (r)] 



This formula can be interpreted as a decomposition of the amplitude K(q,t f ;q,t) 
into a sum of subamplitudes K a (q,t f ;q,t) to which only homotopic loops con¬ 
tribute. With such a decomposition, it is clear that in principle we can assign 
different weights to the different subamplitudes K a {q^t fm ^q^t) provided that we 
preserve the conventional rules for the composition of probabilities. Thus, instead 
of (2.6) we can write 


K(q,t;q,t') 



X(a) f 

Jq Q (t)=q;q*(t')=q 



e'k I! dT £ fa«( T )’S«( T )] 



where x(a) are some complex numbers. For (2.7) to make sense as a probability 
amplitude, the weights x( a ) cannot be arbitrary. In fact, since we want to mantain 
the usual rule for combining probabilities 



the weights x( a ) must satisfy 


x(«i)x(«2) = x(«i • «2) (2.9) 

for any ai and 0 : 2 . Equation (2.9) can also be read as the statement that x(a) 
must be a one-dimensional representation of the fundamental group 7Ti (Af^). 

To see which representations are possible, we have to specify better what are 
Mfi and its fundamental group. To this end, let us consider a system of N identical 

4 In the set 7Ti one can define a product * in a very simple and natural way: if ai and 
G !2 are two classes with representatives qi and £2 respectively, then ai * a 2 is the class 
whose representative is the loop qi £2 (that is the loop qi followed by the loop 52 )* It 
can be shown that this product furnishes 7Ti with a group structure. 
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hard-core particles moving in the euclidean d-dimensional space, IR d . A configu¬ 
ration of such a system is clearly specified by the N coordinates of the particles, 
that is by an element of (Si d ) N . However because of the hard-core assumption 5 , 
any two particles cannot occupy the same position, which means that from (JR d ) N 
we have to remove the generalized diagonal 

A = {(ri,...,rjv) e (JR d ) N : rj = rj for some I ^ jj . (2.10) 


Furthermore our particles are identical and indistinguishable, and thus we should 
identify configurations which differ only in the ordering of the particles. In other 
words we should divide by the permutation group Sn • Therefore we conclude that 
the configuration space for our system is 



( 2 . 11 ) 


To find the fundamental group of such a space is a standard problem in algebraic 
topology, which was addressed and solved in the early 60’s (Fadell and Neuwirth 
1962; Fox and Neuwirth 1962; Fadell and Van Buskirk 1962). Here we simply quote 
the results and refer the reader to the specialized literature for their derivation. It 
turns out that the fundamental group of Mfj is given by 


ail < 2i2 > 

where Bn is Artin’s braid group of N objects which contains the permutation 
group Sn as a finite subgroup (Artin 1926, 1947). 

Even from this formal point of view we see again that there is a crucial difference 
between two and three or more dimensions, exactly as we discovered in the simple 
example at the beginning of this chapter. To have a more explicit understanding 
of (2.12), let us reconsider our two-particle example in the light of what we have 
just observed, following the discussion of (Leinaas and Myrheim 1977). Let us 
start with the case of two dimensions. Instead of assigning the position vectors ri 
and 7*2 for the two particles, it is more convenient to introduce the center-of-mass 
coordinate 

R=\{ri+r 2 ) SIR 2 , (2.13) 

and the relative coordinate 

r = r*i — 7*2 G IR 2 — {0} . (2.14) 


(We have removed the origin because of the hard-core requirement n ^ 7*2.) Since 
R is invariant under the permutations of S 2 , we can write 

= JR 2 x r\ (2.15) 


5 It will be clear later on why this assumption is important to discuss non trivial statistics. 
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where r 2 is “some” space describing the two degrees of freedom of the relative 
motion. We now argue that r 2 has the topology of a cone. Since two configurations 
which differ only in the ordering of the particle indices are indistinguishable, r and 
—r must be identified. The space r 2 is then the upper-half plane without the origin 
and with the positive x-axis identified with the negative one. The results of such 
identification is a cone without tip (see Fig. 2.4). 



Fig. 2.4. The manifold rf describing the relative coordinate of two identical particles in 
two dimensions has the topology of a cone without the tip. 

According to the decomposition (2.15), any loop in can be classified by the 
number of times it winds around the cone r 2 . Two loops q and q f with different 
winding numbers are homotopically inequivalent: It is not possible to deform one 
into the other since the vertex of the cone has been removed. Thus the spaces r 2 
and IR 2 x r 2 are infinitely connected, and 

7 n (M|) =7T! (IR 2 x r 2 2 )=TL = B 2 . (2.16) 

It is important to realize that if the vertex of the cone were included - that is if 
we allowed particles to occupy the same position - the configuration space would 
be simply connected since any loop, even when winding around the cone, could 
be contracted to a point by deforming and unwinding it through the tip. Thus, if 
we do not impose the hard-core constraint on the particles, we can describe only 
bosonic statistics. 

Let us now turn to the case of two particles in three dimensions. After intro¬ 
ducing the center-of-mass coordinate R G IR 3 , we can decompose the configuration 
space M 2 as 

Ml = IR 3 x r 3 (2.17) 

where the space describes the three degrees of freedom of the relative motion. 
These are the length and the two angles of the relative coordinate r, and as before 
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r and — r are identified. It is easy to realize that rf is just the product of the semi¬ 
infinite line describing |r| times the projective space V 2 describing the orientation 
of ±r/|r|. In turn V 2 can be described as the northern hemisphere with opposite 
points on the equator being identified (see Fig. 2.5). 



Fig. 2.5. The projective space V 2 can be represented as the northern hemisphere with 
opposite points on the equator being identified. 

The space V 2 is doubly connected and admits only two classes of loops: Those 
which can be shrunk to a point by a continuous transformation and those which 
cannot. In Fig. 2.6 we exhibit a typical contractible loop q\ and a typical non- 
contractible loop q 2 . 

There are no other classes because the square of a non-contractible loop is 
contractible as Fig. 2.7 shows in a specific example. 

Therefore from the decomposition (2.17) and the topology of r|, we deduce 
that 

TTi (Mi) = TTi (m 3 x r|) = Z 2 = s 2 . (2.18) 

Thus only bosons and fermions can exist, the former corresponding to contractible 
loops and the latter to non-contractible loops. 

We have seen that at the heart of anyonic statistics there is the braid group Bn 
in place of the permutation group Sjn which is responsible for ordinary statistics. 
In fact, there are just two one-dimensional representations of Sjn (the identical 
one and the alternating one, corresponding respectively to bosonic and fermionic 
statistics), whereas the braid group Bn admits a whole variety of one-dimensional 
representations whose labelling parameter will be identified with the statistics v . 

A few remarks are in order at this point. So far we have been focusing on 
one-dimensional representations of 7Ti (M$) because we are considering the sim¬ 
plest possible case of scalar quantum mechanics, when the wavefunctions are 
one-component objects. Of course we could consider more complicated situations 
where wavefunctions are multiplets; in such cases of course appropriate higher¬ 
dimensional representations of 7Ti (Mjv) would be necessary. In the past few years 
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Fig. 2.6. Examples of contractible (qi) and non-contractible ( 52 ) loops. Notice that even 
though q 2 is depicted as an open curve, it is really a closed loop because the two end 
points on the equator are identified in V 2 . 



Fig. 2.7. The square of a non-contractible loop, q f, is contractible. In fact, the con¬ 
tractible loop passing through x can be considered as a deformation of qrf. 

higher-dimensional representations of the braid group B jy have been seen to play 
important roles in a variety of different contexts. Among these we can mention the 
study of the monodromy properties of the iV-point conformal blocks in conformal 
field theories (Moore and Seiberg 1989); the theory of links, knots and Jones poly¬ 
nomials (Witten 1989); the exactly solvable models and the Yang-Baxter equation 
(see for instance (Yang and Ge 1989); the theory of non-abelian anyons (Moore 
and Read 1991). On the other hand, higher-dimensional representations of the 
permutation group were considered long ago in the context of parastatistics. 
In this respect it should be clear that anyonic statistics (which is typical of two 
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dimensions) and parastatistics (which can be defined in any dimension) are two 
completely different concepts, because they are based on two different structures, 
the braid group and the permutation group respectively. Furthermore it can be 
rigorously proven that parastatistics cannot be realized in a non trivial way within 
a local theory for d > 2 (for an extensive discussion of this result see (Frohlich and 
Marchetti 1988, 1989; Frohlich et al 1989)). More simply this means that particles 
obeying parastatistics and having local interactions can be described as ordinary 
bosons or fermions with new degrees of freedom. Indeed the famous puzzle in the 
quantum numbers of the baryonic resonance Zl ++ (for which the idea of paras¬ 
tatistics was introduced in physics) can be solved either with the assumption that 
quarks obey parastatistics (Greenberg 1964) or equivalently with the assumption 
that fermionic quarks carry a new quantum number, the color (Han and Nambu 
1965; Nambu 1966). 

To fully appreciate the new features of anyonic statistics as compared to ordi¬ 
nary statistics, let us now have a closer look at the braid group. The braid group 
of N strands, Hjv, is an infinite group which is generated by N — 1 elementary 
moves <7i,..., satisfying 


<7/ <7/+l Ol = CT/+1 <7/ CT/+1 (2.19) 

for I = 1,..., N — 2, and 

ai gj = gj gj (2.20) 

for \I — J\ > 2. The inverse of gj is denoted by a/ -1 , the identity by 1, and the 

center of B jy is generated by (<7i... gn-\) N • To describe the elementary moves <7/, 
it is convenient to introduce a pictorial representation as follows. Given N vertical 
strands, the generator gj acts on them by simply braiding the J-th strand around 
the (I + l)-st in a definite way, as shown in Fig. 2.8. 



Fig. 2.8. Pictorial representation of the elementary move ai . 

The inverse generator gj~ 1 is represented by the move shown in Fig. 2.9. 
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Fig. 2.9. Pictorial representation of 07 



Fig. 2.10. Pictorial representation of the braid of 07 ^1 • 


A generic braid is a word in the generators 07 and their inverses aj" 1 , which 
however can be rearranged using the relations (2.19) and (2.20). For example, for 
N = 3 the braid represented in Fig. 2.10 corresponds to al 07 07 (we adopt the 
convention that the generators on the right act first). 
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We remark (and Fig. 2.10 shows it explicitly) that in general a/ 2 ^ 1. If 
ai 2 = 1 for all J, then the braid group reduces simply to the permutation group 
Snj which is a finite group. Diagrams like those in Figs. 2.8, 2.9, or 2.10 can also 
be interpreted as describing the time evolution of identical particles if we regard 
the strands as world-lines. To see this we can consider an easy example with three 
particles. Let us suppose that at time t their configuration is the one displayed in 
Fig. 2.11. We can describe it by listing the azimuthal angles of all possible pairs 
of particles measured with respect to some arbitrary reference axes. In the specific 
case of Fig. 2.11 we have 

<Pl2(t) = 0 , 

<Pi3(t) = V , (2.21) 


<P23(<) = £ • 

The symbol <pu denotes the azimuthal angle of particle J with respect to particle 
/, namely 

-l 


<Pij = tan 


m2 m2 

X J X I 

m 1 — ml 

Xj Xj 


( 2 . 22 ) 


(x},x 2 ) being the cartesian coordinates of the J-th particle. 



Fig. 2.11. Configuration of three particles at time t. 

Let us then suppose that at time t f the particles reach the positions shown in 
Fig. 2.12. Now the winding angles are 

<pi2(0 = £ + 7r , 

= T) + TV , (2.23) 

<P23 (<') = 7T . 

The configuration at time t in Fig. 2.11 and the one at time t f in Fig. 2.12 
are the same even if particles 1 and 3 have been interchanged. This is because all 
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Fig. 2.12. Configuration of the three particles at time t f 

particles are indistinguishable! Notice that (pu(t f ) — <pu{t) is not in general an 
integer multiple of 7r; however it is always true that 

^2 <pu(t f ) - ^2 ¥*•*(*) = n r ( 2 - 24 ) 

i<j i<j 

where n is an integer (in our example n = 3). This can be interpreted by saying 
that to complete a loop in configuration space an integer number of exchanges is 
always necessary. We can see explicitly this fact in our example if we place Fig. 
2.12 on top of Fig. 2.11 and draw lines to connect the initial and final positions of 
each particle. In a side-view this is represented by the braiding of Fig. 2.13, which 
is <7 1 02 o\ in terms of the elementary moves. 

We see that indeed there are n = 3 generators, i.e. n = 3 exchanges. Notice 
that if we considered the braiding of 02 Oi (see Fig. 2.10) instead of o\ 02 01 , 
the final configuration would look the same as in Fig. 2.12, but the winding angles 
would be different (specifically (p\ 2 {t*) would be £ + 37T instead of £ + 7r). If we 
just look at the particle positions on the plane, there is no way of distinguishing 
between these two cases. They correspond to the same permutation of the initial 
positions, but if oj ^ 1, they are dynamically different (compare Fig. 2.10 and 
Fig. 2.13). Again we rediscover that the actual braiding of the particles, and not 
simply their permutations, are the important things to be considered in discussing 
anyonic statistics. 

Let us now come back to the problem of finding the one-dimensional represen¬ 
tations x of the braid group. These are given by 




(2.25) 
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Fig. 2.13. Pictorial representation of the braid <j\ <?2 

for any K = 1,..., TV — 1, where v is a real parameter defined modulo 2, which 
will be identified with the statistics. Since in general ok 2 ^ 1, v is an arbitrary 
number. In the elementary move ok all winding angles <pu remain constant except 
for (Pk,k- f l which changes by 7T. Thus we can rewrite x( a ic) as 


x(°k) = exp [— i v A<p KiK+ i] = exp 



(2.26) 


where we have introduced 

= 7T Sl,K$J,K +1 (2.27) 

to represent the change of the winding angles in the move ok * Using this notation, 
(2.26) can be easily generalized to an arbitrary braiding a, and hence we can write 


X(a) = exp 



(2.28) 


where the increment ') — has been written as an integral over an 

evolution parameter r. Notice that the functions are in general very com¬ 

plicated and can be specified only when the dynamics of the particles is fully taken 
into account. Therefore at the moment, we have to regard (2.28) just as a formal 
definition which however is extremely useful. In fact, substituting (2.28) into (2.7), 
we obtain 
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If we now define 




(2.30) 


we see that the kernel K(q i t;q 1 t f ) is decomposed with respect to C! into subam¬ 
plitudes each of which is weighted equally as if we were describing bosons, that 
is 


K{q,t'\q,t) = ^ f 


q a (t)=q;q a (t , )=q 


Vq a 


ki: 


dr £' 


(2.31) 


In other words, instead of dealing with anyons governed by the Lagrangian £, we 
can deal with bosons whose dynamics is dictated by the new Lagrangian £' given 
in (2.30). Thus, we can trade anyonic statistics for some kind of “ficticious” force 
and describe anyons as ordinary particles (bosons for example) with an additional 
statistical interaction. Notice that this statistical interaction is very peculiar and 
intrinsically topological in nature - in our example it is actually a total derivative. 
Its addition to the lagrangian £ clearly does not change the equations of motion, 
which are a reflection of the local structure of the configuration space, but does 
change the statistical properties of the particles, which are instead related to the 
global topological structure of the configuration space. In the next chapter we 
will give an explicit derivation of such statistical interaction, we will clarify its 
topological character and see that it can be locally realized as a gauge theory with 
a Chern-Simons kinetic term. 

So far we have considered N indistinguishable particles moving on the plane, 
and have seen that there exist infinitely many one-dimensional representations of 
the braid group Bn labeled by a continuous statistical parameter v on which a 
priori no restrictions are made. We now briefly consider the case in which N indis¬ 
tinguishable particles move on a compact Riemann surface E. The corresponding 
braid group is 


Bn(E) = 7Ti 



(2.32) 


and has been classified for all kinds of surfaces E (Birman 1969; Scott 1970). As 
far as fractional statistics is concerned, it turns out that the topology of E plays 
a very important role in that it restricts the possible values of v , the number of 
particles N and the number of components of the wavefunctions. 

Let us first consider the case of the sphere (E = S 2 ). The braid group Bn(S 2 ) 
is generated by aj with I = 1,..., AT — 1 which satisfy the Artin relations (2.19) 
and (2.20) plus an additional constraint: 


<7i (72 


_2 

a N-l 


(72 (7i 



(2.33) 
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This expresses the fact that a closed loop taking one particle around all other 
particles can be continuously deformed and shrunk to a point on the back side of 
the sphere. If we use (2.25) inside the relation (2.33), we obtain 

e ~i 2(JV— 1) «/ ir _ j ^ (2.34) 

which immediately restricts v to be rational. If v = p/q with p and q coprime 
integers, (2.34) is equivalent to a restriction on the number of particles, namely 

N = qn + 1 (2.35) 

where n is a non-negative integer. Clearly for bosons (u = 0) or fermions (u = 1) 
there are no restrictions coming from (2.34). On Riemann surfaces of higher genus, 
the restrictions are even more complicated. A Riemann surface of genus g, S gj can 
be described as a sphere with g handles (see for instance (Farkas and Kra 1980)) 
each of which is characterized by a pair of non-contractible curves (homology 
cycles). Hence it is not difficult to realize that in this case besides the exchange 
generators <7/, the braid group B^{S g ) has also “topological” generators which 
take the particles along the 2 g homology cycles around the handles. Following the 
standard notation, we call these generators t^j and pij where I = 1,..., N — 1 and 
i = 1, The generators < tj , Tij and pij satisfy the Artin relations (2.19) and 
(2.20) and 

7t,/+i = oi nj aj l , 

P*j+i = °i Pi,i a 7 l > 

Pi,I+l T »J+1 P*J +1 T i,J+ i = a i > 

9 

(T\ (72 * * • • * * &2 &1 = Pi,l ^*,1 Pi, 1 

i=l 

In particular the last equation is the higher genus generalization of (2.33). We refer 
the reader to the specialized literature (Scott 1970; Einarsson 1990, 1991; Imbo 
and March-Russell 1990) for a proof and a discussion of these relations; here we 
limit ourselves to a brief analysis of the consequences of (2.36) on the quantum 
mechanics of particles of fractional statistics. 

Let us consider the case of a one-dimensional unitary representation of the 
braid group B^{1 7 ff ), as is appropriate for scalar wavefunctions. Every generator 
is represented as a phase, and from (2.36c) it follows immediately that 

l=e- i2u * , (2.37) 

i.e. v = 0,1. Therefore, on a compact Riemann surface of genus g > 1 only 
bosons and fermions are possible if the multi-particle wavefunctions carry a one¬ 
dimensional representation of B^(S g ); this is in contrast with the case of the plane 
where any statistics is allowed with scalar wavefunctions. 

However anyons are possible even on a compact surface provided that wave- 
functions with many components are considered. In this case one has to look at 


(2.36a) 

(2.36b) 

(2.36c) 

(2.36d) 
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higher dimensional representations of Bj^(S g ) which lead to the concepts of gener¬ 
alized fractional statistics and generalized anyons (Einarsson 1990, 1991; Imbo and 
March-Russell 1990; Wen et aL 1990a). These representations are characterized by 
the fact that the exchange moves aj are represented as a phase factor times a 
fc-dimensional unit matrix, namely 



—it/* 



(2.38) 


while the topological moves t^j and pij are represented as unitary k x k matrices. 
The interpretation of such a representation is that the generators aj act on a 
wavefunction with k components by multiplying all components by the same phase 
factor, while the generators t^j and pij mix among themselves the components 
of the wavefunction. Using (2.38) in (2.36a,b), we find 


x(t.\/+i) = x(n,i) = Ti , 

x(Pi,i+ 1 ) = x{Pi,l) = Pi 


Furthermore, (2.36c) implies 



n Pi T ; 


-l 



— 2il^7T 


(2.39) 


or 

Ti pi = e -2 '"* pi n . 

If we now use this last relation in (2.36d), we easily get 

e —i2(JV— l)u* _ e ~2igi/* 


(2.40) 

(2.41) 


which is the higher genus generalization of (2.34). This equation restricts v to be 
a rational number. If v = p/q with p and q coprime integers, then 

N = qn + 1 — g . (2.42) 

For g = 0, we recover the sphere condition (2.35). When v — p/q it is possible 
to derive explicit expressions for the k x k matrices r* and pi . Consistency with 
(2.36) requires that k = q 9 (Einarsson 1990,1991). Thus, on compact surfaces, only 
rational statistics are allowed, and the number of particles and the dimensionality 
of the braid group representations are not arbitrary. For an extensive discussion 
of anyons on compact surfaces and on the torus (g = 1) in particular, we refer the 
reader to the recent review by R. Iengo and K. Lechner (Iengo and Lechner 1992). 
Throughout the rest of this book, we will always consider the case of particles 
moving on the plane, where, as we have seen, no restrictions apply. 

We conclude this introductory chapter by summarizing the main points of our 
discussions. The configuration space of identical two dimensional particles has a 
non trivial topology and is an infinitely connected space. Its fundamental group 
is the braid group whose representations are labeled by an arbitray parameter v 
called statistics. This unusual statistics can be implemented on ordinary particles 
(for instance bosons) by the addition of a topological statistical interaction. 



3. Fractional Statistics in the Chern-Simons 

Gauge 


In Chapter 2 we have seen that fractional statistics can be implemented on ordinary 
bosons by the addition of fictitious statistical interactions (the same of course 
would be true on ordinary fermions). In this chapter we are going to discuss this 
scenario in more detail and find that the use of Chern-Simons gauge fields provides 
an efficient method to realize anyons. Before this however, we consider useful to 
discuss the dynamics of a charged point-particle interacting with an infinitely long 
magnetic solenoid (flux tube). This system has been called cyon (Goldhaber 1982) 
and is the prototype for anyons. When the motion along the solenoid is ignored, 
the dynamics is essentially planar and the cyon is subject to the laws of the two- 
dimensional world. In particular it may acquire fractional statistics. 


3.1 The Cyon System and Its Symmetries 


Let us consider a non-relativistic particle of mass m and electric charge e that 
moves in the magnetic field B created by an infinitely long and thin solenoid 
passing through the origin and directed along the z-axis. If we neglect the (free) 
motion along the solenoid, the relevant dynamics takes place in the (x,?/)-plane 
and is governed by the non-relativistic Lagrangian 

C = \mv 2 + -v * .A(r) (3.1.1) 

2 c 


where r = (x,?/) £ IR 2 denotes the particle position, v = r its velocity and A the 
vector potential for the solenoid configuration. In a convenient symmetric gauge, 
A is given by 



27r 



(3.1.2) 


where i and j are the unit vectors along the x- and the y-axis respectively, and 
# is the flux of the solenoid. The magnetic field B (which in two dimensions is a 
pseudo-scalar) associated to (3.1.2) is that of a vortex of strength # localized at 
the origin, 

B = VAA = <f <5 (2) (r) (3.1.3) 

where the 6-function arises from the singularity at x = y = 0 of the vector potential 
(3.1.2). Clearly we have 
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d 2 r B = $ , (3.1.4) 

confirming the interpretation of # as the flux of the solenoid. Following (Goldhaber 
1982) we will call this particle-vortex system, a cyon. The canonical momentum 
p can be easily derived from the Lagrangian (3.1.1); it differs from the kinetic 
momentum mv by a term due to the vector potential A, namely 

dC e 

p = = mv + -A 

ov c 

The Hamiltonian associated to (3.1.1) is 


(3.1.5) 



H = p v-C= =lmv 2 . (3.1.6) 

Notice that the magnetic field and the vector potential are invisible when the 
Hamiltonian is written in terms of the kinetic momentum. In fact, as we see in 
(3.1.6), the Hamiltonian of the cyon is numerically equal to that of a free particle! 
The only effect of the presence of the vortex is to enforce the non-trivial relationship 
(3.1.5) between the canonical and the kinetic momenta. The classical equations of 
motion are the same as those of a free particle, but non-trivial features are present 
in the quantum theory. 

The Lagrangian (3.1.1) is rotationally invariant: If the dynamical variables are 
rotated in the plane, C does not change. The constant of motion which is associated 
to this rotational symmetry is the canonical orbital angular momentum J c , 


J c = rAp = rA mv + -r A A 

c 


= r A mv + 


e# 

27TC 



e# 

27TC 


(3.1.7) 


where J is the gauge invariant kinetic angular momentum. As clearly stated in 
(Jackiw and Redlich 1983), the conserved canonical angular momentum J c has a 
conventional spectrum: Its eigenvalues are always integers in units of h. This is 
true despite the fact that the algebra of the two-dimensional rotations is abelian 
and in principle an arbitrary constant could be added to the angular momentum 
operator, thus obtaining arbitrary eigenvalues. 

The difference between J c and J is entirely due to the magnetic flux # in the 
solenoid. Both in the absence and in the presence of the canonical angular 
momentum is always represented by the quantum mechanical operator 



(3.1.8) 


where ip is the polar angle on the plane. When it acts on single-valued wavefunc- 
tions with angular dependence exp(im^>), it becomes simply 


J c = hm , m € TL . 


( 3 . 1 . 9 ) 



25 


Thus the spectrum of J c is always the conventional one (Jackiw and Redlich 1983). 
Of course when # = 0 and hence J c = J, also the kinetic angular momentum has 
only integer eigenvalues. However when # ^ 0, this is not true any more. In fact 
using (3.1.7) and (3.1.8), we see that the kinetic angular momentum operator is 


J = J C - 


e$ 
27 vc 



e$ 
27 vc 


(3.1.10) 


which, when acting on single-valued wavefunctions with angular 
exp(im^>), becomes 


J = k 



m € TL . 


dependence 


(3.1.11) 


Therefore the spectrum of J consists of integers shifted by — e$fhc . 

Usually in quantum mechanics we never consider the kinetic angular momen¬ 
tum, since the conserved quantity is the canonical one. However for the present 
situation, A.S Goldhaber and R. Mackenzie (Goldhaber and Mackenzie, 1988) lu¬ 
cidly pointed out that the integer canonical angular momentum is divided into 
two pieces: One which is localized near the cyon and is in general fractional, and 
one which is located at the spatial infinity and is also fractional. Furthermore they 
argued that this diffused angular momentum is irrelevant in describing local phe¬ 
nomena, and identified the piece localized on the cyon with the kinetic angular 
momentum. 

We can reformulate their argument as follows. First of all let us observe that 
the following identity holds 


J c = J — r A A 

c 


= J — ^ J d 2 x x • E(t , x ) B(t, x) + - J d 2 x V • [ E(t , x) x A A] . 


(3.1.12) 


In this formula E(t, x) represents the electric field created by the moving charge 
which satisfies the Gauss law V-£?(t,x) = e 6( 2 )(x —r(t)) where r(t) is the particle 
position at time t. 6 When we substitute (3.1.3), we find that the second term in 
the right-hand-side of (3.1.12) identically vanishes, so that for the cyon system we 
are left with 

J c = J + ^ J d 2 x V • [E(x) x A A] . (3.1.13) 

The difference between J c and J is just a surface term which however cannot 
be neglected tout-court - indeed as we have seen its value is e$/2ixc ! We can 
understand and clarify the role of such surface term if we imagine slowly turning 
on the flux through the solenoid from a zero initial value to the final value 
When # = 0, J c = J and both are integers. When # is turned on, only J c remains 
constant and thus integer. Its conservation is in fact a consequence of the rotational 


6 The reader may verify that by computing the divergence inside the last integral of 
(3.1.12) one produces the | r A A term of J c upon using the Gauss law, plus a term of 
the form J d 2 x E(tjx) • V(sc A A) which exactly cancels the first integral in (3.1.12). 
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invariance of the system, for which J c is the conserved Noether charge. When the 
flux reaches its final value a piece of J c (specifically the surface term in (3.1.13)) 
is radiated away to infinity by the gauge fields. For phenomena on a finite length 
scale, after a finite time, this “dissipated” piece of the angular momentum is no 
longer relevant, and only J is left on the cyon. Thus even though the total canonical 
angular momentum remains integer, the cyon retains just a part of it which is in 
general fractional. 

Following (Wilczek 1982a,b, 1990), we call this (kinetic) angular momentum 
the spin of the cyon. More precisely we have 



J(m = 0) 

k 


e# 

he 


(3.1.14) 


In general $ is neither integer nor half-integer. 

If some connection between spin and statistics exists, we should expect that 
the cyon be in general an anyon. To establish their statistical properties, let us 
consider two identical cyons with a wavefunction ^>(1,2), and assume that the 
magnetic flux and the electric charge are tightly bound on each particle. Now let 
us slowly move one cyon around the other by a full loop and neglect both charge- 
charge and vortex-vortex interactions. As it was discovered by Y. Aharonov and D. 
Bohm (Aharonov and Bohm 1959), even though the presence of a vortex does not 
effect the classical dynamics of a charged particle, it has profound consequences 
on its quantum mechanics. Specifically in our case when, say, particle 1 is moved 
around the vortex 2 on a closed loop r , the wave function acquires a phase 



(3.1.15) 


Using Stoke’s theorem, (3.1.15) can be rewritten in terms of the magnetic flux, 
namely 



(3.1.16) 


In the two-cyon system, when the particles are rotated around each other, there 
are actually two contributions to the phase: One due to the motion of the first 
charged particle in the vortex field of the second, and one due to the motion of the 
second charged particle in the vortex field of the first. The total phase acquired 
by the wavefunction ^>(1,2) under a full 2ir rotation is then 


exp 



(3.1.17) 


If we compare (3.1.17) with (2.1) in Chapter 2 and recall that in this case A<p = 2i r, 
we can easily deduce that the statistics of the cyon is 

2e# 


v = — 


he 


(3.1.18) 
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Moreover the spin $ and the statistics v appear to be related in the conventional 
way 

v = 2s . (3.1.19) 

Thus in general the cyon is an anyon and the standard spin-statistics connection 
is satisfied! In particular, a bosonic particle with a flux # = (where = hc/e 
is the fundamental flux unit) behaves effectively as a fermion. We conclude this 
section by mentioning that the cyon system possesses not only spatial rotational 
symmetry, but also a conformal symmetry (Jackiw 1990). In fact, one can check 
that under the following transformations 

Sr = v , 

Sr = tv — \-r , (3.1.20) 

Sr = t 2 v — tr , 


the Lagrangian (3.1.1) changes only by a total derivative and thus (3.1.20) are 
symmetries of the system. The conserved Noether charges associated to (3.1.20) 
are respectively 

H = \mv 2 , 

2 

D = tH — i(r * v + v * r) , (3.1.21) 

K = —t 2 H + 2 tD + \mr 2 . 

Using the canonical commutation relations and keeping in mind the non-trivial 
relationship between the canonical momentum p and the velocity v (see (3.1.5)), 
one may show that the charges H , D and K close the following algebra 


[D,H] = -ihH , 
[D J K]= ihK , 

[J5T, K\ = 2i HD . 

If we define 

L ±=2K {K ~ H)± 'k D ' 
L a = -±(K + H) , 

then the commutators (3.1.22) can be rewritten as 

[L+j L—\ = 2Lq , 
[Lq,L±] = ^L± . 


(3.1.22) 


(3.1.23) 


(3.1.24) 


In these we recognize the commutation relations of 50(2,1) in the standard Cartan 
basis. We recall that 50(2,1) is the finite subalgebra of the Virasoro algebra which 
describes the conformal transformations in two dimensions. We will see in Chapter 
9 that there is a remarkable even though mysterious connection between anyons 
and conformal field theory. 
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3.2 Chern-Simons Construction of Fractional Statistics 

In the previous section we have seen that a particle which carries both a charge e 
and a flux # possesses a spin $ 

(3.2.1) 

which is neither integer nor half-integer in general, and is an anyon of statistics 
v = 2 $. In this section we are going to present a simple construction of fractional 
statistics which generalizes these features to a many-body system by using Chern- 
Simons gauge fields. The use of the Chern-Simons theory to realize fractional 
statistics was first suggested by F. Wilczek (Wilczek 1982a,b) and subsequently 
clarified by many others (see for example Arovas et al 1985; Hansson et al 1988). 
The essential feature of the Chern-Simons formulation is that fractional statistics 
is implemented by means of a local, long-range, abelian gauge interaction taking 
place in a (2+l)-dimensional space-time. In fact rather than affixing by hand an 
externally prescribed vector potential A to an ordinary particle, which is then 
transmuted into an anyon, one promotes A to be the space component of a dy¬ 
namical (2+l)-dimensional Z7(l) gauge field A a whose action is taken to be 

Scs = J d 3 x C C s = Y c f d * x e a ^A a d^A y (3.2.2) 

where e a/37 is the completely antisymmetric tensor density. Here and in the fol¬ 
lowing we adopt the convention that Greek indices take the values 0,1,2 and are 
contracted with the metric rjap = diag(l, —1, —1). The space components are in¬ 
stead labeled by Latin indices taking the values 1,2. The three-vector x a = ( ct,x l ) 
is denoted for short also as x, so that d 3 x = cdt d 2 x. 

The expression (3.2.2) is known as Chern-Simons action (Siegel 1979; Schon- 
feld 1981; Jackiw and Templeton 1981; Deser et al 1982a,b; Hagen 1984). It 
is gauge invariant even if the Lagrangian contains an undifferentiated gauge 
field A a . This is so because the gauge potential is contracted with the current 
e <*£7 dp A y = 1/2(where as usual Fp y = dp Ay — dyAp) which is conserved 
because of the Bianchi identity. In other words under gauge transformations 

A a —► Aa + d a A (3.2.3) 

where A is a space-time dependent parameter, the Chern-Simons Lagrangian 
changes only by a total derivative 



S£cs = ~ e a ^ d a [{dpA) Ay\ , (3.2.4) 

z c 

so that the action (3.2.2) remains invariant. 

The Chern-Simons term as written here is special to (2+1) dimensions, but it 
can be easily generalized to any odd dimensional space-time and appears in a host 
of different contexts: anomalies (see for example (Alvarez-Gaume and Ginsparg 
1985)), supergravity (see for example (Ceresole et al 1986)), string theory (see for 
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example (Green and Schwarz 1984)) and so on. The Chern-Simons term provides 
also an example of topological field theory (for a review see (Birmingham et al 
1991)) since even in a curved space-time, the action has the same form as in 
(3.2.2) without any additional metric insertions (the e-symbol is indeed a good 
tensor density). We refer the reader to the extensive literature on these numerous 
applications of the Chern-Simons term, and concentrate only on its use to construct 
anyons. 

To this end, let us couple the gauge field A& to a matter system consisting 
of N non-relativistic point particles of mass m and charge e 7 , whose coordinates 
r/(t) serve as dynamical variables (capital Latin indices take the values 1,..., N 
and label the particles). The current 

N 

j a ( x ) = ^2 ev ?( t ) $ (2) (* - r i ( t )) = (cp,j) (3.2.5) 

7=1 

where vf (t) = (c, v/(t)), clearly satisfies the continuity equation 

dj* = dtp + V • j = 0 . (3.2.6) 

The meaning of j a becomes more evident if we write explicitly the time and space 
components, namely 


N 

p( x ) = ^2 esW ( x - r i( t )) > 

7=1 

N 

i( x ) = ]C e v i( t ) s{2) ( x - *v(0) 

7=1 


(3.2.7) 


We recognize in p and j respectively the conventional charge density and current 
density for N non-relativistic point particles located at rj(t) and moving with 
velocities ^ 7 (t) = f*/(t). 

Then we couple the conserved current j a to the gauge field A& in the standard 
minimal way 





d 3 x j a (x)A a (x) 




dt 



[«/(<) -A(t,n(t)) - >lo(<,r/(t))] 


(3.2.8) 


(recall that d 3 x = cdtd 2 x ), and take as kinetic term for the N particles, the 
non-relativistic action 

7 Notice that this charge is not the standard electric charge since the gauge field A a does 
not have the Maxwell action of the standard electrodynamics but the unconventional 
Chern-Simons term (3.2.2). 
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matter 


= /* (e 


(3.2.9) 


The total action for our system is therefore 

S = ^matter + S'mt + ScS = J dt C 

where the total Lagrangian £, spelled out in detail, is 

N r- 

C = Vj -mv 2 j + -«/ • A(t,r/(t)) - -A^t^r^t)) 

1=1 c c 

J d 2 x [E(t, x) A A(<, x) + Ao (t, x)B(t, x)] 

Here we have introduced the Chern-Simons magnetic field 


(3.2.10) 


(3.2.11) 


B = V A A = diA* - d 2 A' = -F 12 , 


(3.2.12) 


and the Chern-Simons electric field 


E l = 


— d t A l - diAo = F oi 

c 


(3.2.13) 


(We recall that space indices are lowered and raised with the metric rjij = 
diag(—1,-1).) 

By varying C with respect to r/(£), we obtain the equations of motion for the 
matter variables 


mv}(t) = e ^(«,r/(«)) + ^vj B(t,n(t)) 


(3.2.14) 


where e 12 = —e 21 = 1. These are the standard Lorentz force equations for particles 
of mass m and charge e moving in an electric field E and in a magnetic field B. 
By varying C with respect to A& we get instead the equations of motion for the 
gauge fields. However, because of the peculiar nature of the Chern-Simons term 
which contains only first-derivatives of A&, more than true equations of motion 
these are actually identities that relate the fields E and B to the matter current 
j and p ; indeed one gets 


or in components 



(3.2.15) 

1 t 

E' = —e t 3 j 3 , 

(3.2.16a) 

KC 


B = —— p . 

(3.2.16b) 


K 


These equations are very important. First of all, they tell us that given a matter 
configuration (that is given p and j), the field strengths of the Chern-Simons fields 
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are not arbitrary but are uniquely prescribed according to (3.2.16). This is very 
different from the cyon system of Section 2 where the magnetic field B was an 
externally given field and could be chosen to be that of a vortex. Secondly, if we 
use (3.2.16) for a single particle (N = 1 and subscript I suppressed), it is easy to 
prove that 


E % {x) + — e x *v*{t) B(x) = — e* J (j 3 (x) — rf(t) p(x)) 

= —e* J (t^(£) — t^*(t)) £^(® — r(£)) = 0 
/c c 


(3.2.17) 


In going from the first to the second line, we have used the appropriate explicit ex¬ 
pressions (3.2.7). Thus, a particle moving in a Chern-Simons background does not 
experience any self-interaction since the Lorentz force is vanishing. However, due 
to the non-trivial relationship between the velocity v and the canonical momentum 
p, z.e. 

p = mv + -A(r ) , (3.2.18) 

c 

the presence of the Chern-Simons electric and magnetic fields is not at all im¬ 
material. For an N particle system, the equations of motion are explicitly given 

by 

2 N 

mv} = — ^2 e' 3 (v^jit) - v}(t)\ S^Xr^t) - rj(t )) . (3.2.19) 

KC j =1 

Notice that for J — I the coefficient of the ill-defined £^( r /(^) — r /(^)) on the r.h.s. 
of (3.2.19), identically vanishes, and hence no ambiguities and no self-interactions 
are present in the system (Jackiw 1990). 

The vanishing of the Lorentz force does not come as a total surprise (De Sousa 
Gerbert 1990). In fact as is well-known, (3.2.14) can also be written in covariant 
notation as 

d a T^ ttei = -±f (3.2.20) 

rv C 

where 7^f tter is the energy-momentum tensor of the matter system and is the 
force density given by 


f = = {E i j i , + e ik j k B) . (3.2.21) 


Using (3.2.16), we easily see that f 13 = 0 identically, and thus 7^f tter is conserved. 
The vanishing of and hence the conservation of T^f tter are necessary because 
the Chern-Simons action, being topological and metric independent, does not con¬ 
tribute to the energy momentum tensor. So 7^f tter is the full energy momentum 
tensor of the whole system which better be conserved! 

The most remarkable consequence of the Chern-Simons interaction is however 
the relation (3.2.16b) between the magnetic field B and the charge density p which 
we rewrite here for convenience 



(3.2.22) 
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If we integrate this equation over a small two-dimensional disc Cj which includes 
only the J-th particle, the left-hand-side yields the magnetic flux attached to that 
particle while the right-hand-side, upon using the explicit expression (3.2.7), gives 
its charge, namely 

r r N 

<Pj = I d 2 x B = —— I d 2 x\^6^ (x — rj(t)) =—— . (3.2.23) 

JCi K Jci K 

Therefore, whenever a particle possesses a charge e, it necessarily possesses also a 
flux # = —e//c. The Chern-Simons dynamics automatically binds charge and flux 
to the particles, and thus realizes in a natural way what had to be imposed as an 
external requirement in the cyon system. We can picture the situation as follows: 
each particle possesses both a charge and a flux, and moves in the background fields 
created by the other particles; the Chern-Simons electric and magnetic fields con¬ 
spire to cancel all self-interactions and the resulting Lorentz force is zero. However, 
since a charge interacts with all other fluxes and vice versa, we expect interesting 
quantum effects from the Aharonov-Bohm mechanism. 

The Hamiltonian arising via canonical procedure from the total Lagrangian 
(3.2.11), is 



(3.2.24) 


If p = —k B is imposed as a constraint, Aq may be set equal to zero by choosing 
the Weyl gauge, and hence H becomes numerically equal to the Hamiltonian of N 
non-interacting particles 

N /i \ 

h = J2 ■ ( 3 - 2 - 

As in the cyon system, the non-trivial dynamics resides entirely in the relation 
between the canonical and the kinetic momenta. The condition Ao = 0 does not 
fix the gauge completely since we are still free to perform any time-independent 
gauge transformation without leaving the Weyl gauge. This residual freedom is 
removed by imposing the subsidiary condition 




(3.2.26) 


Now the constraint p = —k B = Ke tJ diAj with p given in (3.2.7), can be unam¬ 
biguously solved and one finds 


Aj{v\ , ... t * If} 



(3.2.27) 


Thus, by eliminating the Chern-Simons gauge field A a through its field equations 
(or more properly through the constraints (3.2.16) after gauge fixing), we end 
up with a system of particles moving in the effective non-local vector potential 
(3.2.27). The Hamiltonian (3.2.25) can then be rewritten as follows 
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A,(r, ’- r ' r 




(3.2.28) 


where the prime on H denotes that we used the solution given in (3.2.27). We 
remark once again that this effective vector potential is non-local since it depends 
on the positions of all particles, and in particular when N = 1 it vanishes. If we 
compute the magnetic field associated to Aj(ri, ...r^), we get 



eij -£r Ai i( r 1 »- p *) = 


--^5 (2) (r/ -rj) 
K J*i 


(3.2.29) 


Each particle sees the (N — 1) others as vortices carrying a flux # = — e//c. This is 
in agreement with our previous derivation (3.2.23) based on the use of the Chern- 
Simons constraints. 

The non-local potential Aj(ri, ...r^r) seems badly divergent at first sight. How¬ 
ever this is not the case because the would-be singularities occur only at coincident 
points (rj = rj for some I and J) which do not belong to the configuration space 
of the system (see (2.11)). 

It is interesting to observe that the solution (3.2.27) can be written also as 




2 ™ dr] 
e d 


- rj 


27 r/c dr} 

Jr I 


(3.2.30) 


where <pu is the winding angle of particle J with respect to particle I, such that 


<pu = tan 


l ( x 2 i~ x 2 j \ 


(3.2.31) 


(see (2.22)). The Lagrangian corresponding to (3.2.28) is then 


£ ' = E 

7=1 

N 

= £ 

7=1 

N 

= E 



where in deriving the last term we used the property 



(3.2.32) 


(3.2.33) 


If we now observe that 
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d_ 

dt 


<PIJ 



£' becomes simply 



(3.2.34) 


(3.2.35) 


We managed to cast our Lagrangian in the form of (2.30) and so we immediately 
deduce that our particle are generically anyons of statistics 



2'kUck 


e$ k 
2n he 2 it he 


(3.2.36) 


where the second two equalities follow from the charge-flux relation (3.2.23). Thus 
we can turn bosons into fermions, fermions into bosons or in general anyons into 
other anyons by suitably adjusting the coupling constant k of the Chern-Simons 
term. Notice that the statistics v in (3.2.36) is one half of the one computed in the 
previous section for the cyon system (Hansson et al . 1988; Wen and Zee 1989a). 
It is now interesting to see how this factor of 1/2 comes about. To this end let 
us consider again the total Lagrangian (3.2.11) and observe that under a gauge 
transformation with parameter A 


A a —► A a + dA , 


(3.2.37) 


its variation is 


6C = J d 2 x (-ij“ + d a A 


(3.2.38) 


where j a is the matter current (3.2.5) and the e F term originates from the Chern- 
Simons action (see (3.2.4)). The Noether current associated to the gauge transfor¬ 
mation (3.2.37) is 

J a = -j a - . (3.2.39) 

0 * 

Both terms in (3.2.39) are clearly conserved (because of the continuity equation 
on the matter current and the Bianchi identity on the Chern-Simons current) and 
hence the total action S = J dt C is invariant against U(l) gauge transformations. 
Using the Chern-Simons constraint (3.2.22), we can rewrite the Noether current 
(3.2.39) as 

J a = ^-j a . (3.2.40) 

2c 

If we integrate its time component over a small disc Cj containing only the J-th 
particle, we deduce that its conserved U(l) Noether charge is 


Qi 


= X/* J0 = U/* ' = 5 


(3.2.41) 


The charge associated to the J7(l) Chern-Simons field is not e but 1/2 e! The use 
of the Chern-Simons action action modifies the standard electrodynamic relation 
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by introducing a factor of 1/2. In fact, if we had considered the ordinary Maxwell 
term instead of the Chern-Simons one, the U(l) Noether current would have been 



(3.2.42) 


and correspondingly the J7(l) Noether charge (that is the standard electric charge) 
would have been 



(3.2.43) 


The difference originates from the fact that the Maxwell Lagrangian is invariant 
under gauge transformations, whilst the Chern-Simons Lagrangian varies into a 
total derivative. 

Let us now compute the statistical phase induced by the Aharonov-Bohm effect 
on our particles. Let us consider a two-body wave function ^>(1,2) and suppose 
that particle 2 is moved on a full loop around particle 1. Under such operation, 
due to the interaction between the charge of particle 2 and the flux of particle 1, 
the wave function acquires a phase 



(3.2.44) 


where q is the Noether charge (3.2.41) (we suppressed the index on q because 
the particles are assumed to have the same charge). To get the complete phase, 
we must multiply the exponent of (3.2.44) by two to count also the interaction 
between the charge of particle 1 and the flux of particle 2. In conclusion we obtain 
that the total phase acquired by ^>(1,2) is 


exp 




(3.2.45) 


Using (2.1) and remembering that Atp = 2n for a full loop T, we deduce that 



(3.2.46) 


in perfect agreement with our previous calculation (see (3.2.36)). 

In particular a bosonic particle with a flux # = 1/2 ( where = hc/q = 

2 hc/e is the fundamental flux unit of the Chern-Simons theory) is a fermion. Notice 
that the Chern-Simons flux unit Sq = hc/q is twice the ordinary electrodynamic 
flux unit = hc/e because of (3.2.41). However both in the cyon system and in 
the Chern-simons theory, one can say that a boson is converted into a fermion by 
affixing 1/2 unit of flux! 

The system described by the total Lagrangian (3.2.11) is invariant under rota¬ 
tions of the spatial coordinates 


Sr} = —e tJ rj , 


(3.2.47) 


and of the gauge potentials 
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6A Q (t,x) = -e 3k x 3 d k Ao(t,x) , 

6Ai(t,x) = —e* k x?dkAi(t,x) — e'iAj(t,x) . 

The corresponding Noether charge is the total angular momentum 

N ^ - 

Jtot = ^ r i A Pj + ^ I d 2 x ( x A -A(x)) B(x ) . 


(3.2.48) 


(3.2.49) 


The first term in (3.2.49) is the canonical angular momentum for the N particle 
system, whilst the second term is the explicit contribution of the Chern-Simons 
fields. Using the constraint (3.2.22), we can rewrite (3.2.49) as 

N N 

Jtot =y]^Ap / --y'r/A -A(p/) 

1=1 1=1 

N 

rj A mvj = J . 

7=1 

Thus Jt 0 t coincides with the total kinetic angular momentum of the iV-particle 
system. Furthermore, we observe that, using (3.2.7), J tot can be written in terms 
of the current density as follows 

Tfb £ 

Jtot = J = — / d 2 x x A j(x) . (3.2.51) 

Contrary to the cyon system where we had to argue that the relevant part of the 
angular momentum for long time-scale phenomena was the gauge invariant kinetic 
piece, here in the Chern-Simons theory the conventional derivation shows directly 
that J is the conserved Noether charge for rotations. There is no contradiction in 
this, since, as we have already remarked, for one particle the Chern-Simons gauge 
fields are irrelevant (the effective vector potential (3.2.27) vanishes for N = 1). 
In this formulation in fact, the statistical gauge interaction only occurs among 
couples of different particles. 

From (3.2.50) we expect that the spectrum of J tot = J consists of in general 
non-integer eigenvalues. To see this explicitly, we can use the relation between pj 
and mvj in the effective theory described by the Hamiltonian H l given in (3.2.28), 
and write 



(3.2.50) 


N 

J = ^ rj A mvj 

i=i 



= J c --^2 e* J r}i4j(ri, ...r N ) . 
c 


N 


1=1 


(3.2.52) 


Using the explicit expression of the vector potential given in (3.2.27) and keeping 
in mind that J c has only integer eigenvalues, the spectrum of J is 
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integer + |z/^iV(iV — 1) | = integer + t u (3.2.53) 

where the statistics v = —e^/2'KTic as in (3.2.46) is defined modulo 2, and the 
symbol {x}p denotes the periodic extension of x over a period of 2 (see Fig. 3.1) 
(Dunne et al 1992b). 



Fig. 3.1. The fractional orbital angular momentum l v as function of the statistics v . 


In general the eigenvalues of the angular momentum are neither integers nor 
half-integers (in units of fi). Notice that for v = 1, t v —\ = 0. In fact for any iV, 
N(N — l)/2 is always an integer and hence it can be reabsorbed into a suitable 
shift of J c . Furthermore, for N = 1 l u = 0 for any v. Thus t v has nothing to 
do with the “spin” of the anyons; on the contrary it is really a genuine fractional 
orbital angular momentum. The non-relativistic Chern-Simons theory, as we have 
presented it here, provides a framework within which the spin is an absolutely 
conserved quantity, so that it is perfectly consistent to ignore it, or simply set 
it to zero. Therefore we can say that the non-relativistic Chern-Simons theory 
describes spinless particles with fractional statistics. (We should not be surprised 
at this fact; after all it is customary to consider spinless fermions in non-relativistic 
quantum mechanics!) Since there is no spin for our particles, of course there is no 
spin-statistics connection. In Chapter 5 we will see that there is a natural way to 
incorporate the spin in the theory using the second-quantized formalism. Then, 
we will find that the usual spin-statistics relation is satisfied. 

We conclude our presentation of the Chern-Simons theory by mentioning that 
the system described by (3.2.11) is invariant also under the following conformal 


transformations (Jackiw 1990) 

Srj = vj , 

SA a = cF 0a + d a (c Ao) 


(3.2.54a) 


Sr i = tvj — -rj 

& 




8A a = ct F 0a + d a (ct Ao) - 



(3.2.54b) 
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Srj = t 2 vj — trj , 

SA a = ct 2 Fo a + d a (ct 2 Aq)- tx x F ia - d a (tx x Ai) 


(3.2.54c) 


To be more precise, the Lagrangian (3.2.11) is not invariant under (3.2.54) but its 
variation is a total derivative so that these transformations are true symmetries of 
the system. The conserved Noether charges associated to (3.2.54) are respectively 


H = 


+ J d 2 x A 0 (x) ( kB(x) + p(x)) , 


D 


1 N 

= tH — -m ^ (rj * vj + vi * rj) 


+ 


/a 


1=1 


tA 0 (x) - —x • A(x) ) ( kB(x ) + p(x)) 


(3.2.55) 


K = -t 2 H + 2tD + Y^ 


+ 


J d 2 x (< 2 Aq(x) — • A(x) \ ( kB(x) + p(x)) 


Upon using the constraint (3.2.22), the gauge field contributions to H , D and K 
disappear and one is left with the straightforward iV-particle generalization of the 
conformal generators given in (3.1.21) for the cyon system. From the canonical 
commutation relations it follows that 


[L+,L_] —2Lo ? 
[L 0 ,£±] = =F L ± 


(3.2.56) 


where, as in (3.1.23), 

L ± = - H) ± ID , 

U = -^K + H) . 

This is the algebra of the conformal group 50(2,1) in two dimensions, which is a 
symmetry group of the Chern-Simons theory. 
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There is a very simple way to check if a many-body wavefunction describes bosons 
or fermions: If under an exchange of any two particles acquires a plus sign, 
it describes bosons; on the other hand if it acquires a minus sign, it describes 
fermions. The statistics is therefore encoded in the boundary conditions of In the 
previous chapter we have seen that fractional statistics could be realized by means 
of an effective non-local interaction among particle of ordinary statistics. Now 
we are going to demonstrate that there is an equivalent description of fractional 
statistics where complicated boundary conditions replace the effective interaction 
(Wu 1984b). The two descriptions are related to each other by a (singular) gauge 
transformation and thus can be used interchangeably. 

Let us consider again the Hamiltonian (3.2.28), which we rewrite here for con¬ 
venience 


H ' = E ^ (*/ - s A '( r " 


(4.1) 


where 




2ttk 




3*1 
e d 


r i r j 

1*7 - rj | 


2ttk dr 




(4.2) 




<pu = tan 


-1 


m2 m2 

X I X J 

m 1 _ rr* 1 

X I X J 


This is the first-quantized Hamiltonian for a system of N anyons of statistics 


v = 


2ttkc 


(4.3) 


which are described by a single-valued, symmetric wavefunction such that 


H* = 


(4.4) 


We want to emphasize that I s single-valued and “bosonic” according to our 
definition, because 


-0(ri, * * *, r j, * * *, r/, * * *, r N ) = V , ( r i? * * * 5 *7> * * •, r J, * * *, »7r) * 


(4.5) 
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The non-standard statistics v is explicitly reproduced by the non-local potential 
Aj. However, from (4.2) we see that Aj is a pure gauge and can be removed by a 
gauge transformation. The transformed wavefunction is 


iu^2<pu 

ip = e J < J ip , 


and the transformed Hamiltonian is 





e 


i<j 



so that the eigenvalue problem (4.4) becomes 




H ip = E ip 



H is a free Hamiltonian, but the new wavefunctions are multi-valued and satisfy 
the “twisted” boundary conditions 




= e i7rl/ 




This equation, which follows from (4.6) and tpjj = (fjj + ir, explicitly shows 
that ip carries an abelian representation of the braid group Bn. Extending our 
initial definition, we may therefore say that ip is an “anyonic” wavefunction. It is 
customary call this the anyon gauge description of fractional statistics. 

This result can be retrieved in another way which actually shows how general 
the anyon gauge description is, because it does not rely on the specific features of 
the model. We saw in Chapter 3 that the Lagrangian associated to H f is 



(4.10) 


This is also the general form of the Lagrangian describing anyons of statistics v 
as we discussed in Chapter 2. The kernel K(q / ,t / ;q,t) J which is built out of is 
(see (2.29)) 8 


K(q\t f ;q, 




and propagates the single-valued wavefunctions according to 


(4.11) 


8 


In Chapter 2 we considered for simplicity only closed paths in the configuration space, 
but the construction can be clearly generalized to arbitrary open paths. 
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(4.12) 


Since the “statistical” interaction is a total time-derivative, (4.11) can be rewritten 
as (Wu 1984b) 


oo -iu i * ¥>jj(t) 

K(q',t';q,t) = ^ e '< J K^ n} (q',t';q,t) e << J (4.13) 

n/j=—oo 

The sums over the integers njj are a consequence of the multi-valuedness of the 
angles pj j, and each contribution in these sums corresponds to paths with winding 

numbers njj in the configuration space. Finally, is computed only from paths 
with fixed winding numbers, using the free Lagrangian £. 

Given (4.13) and (4.12), it is natural to redefine the wavefunction and introduce 

= e 1<J 't/j . (4-14) 

Furthermore, if we understand the angles ipu as defined modulo 2n so that the 
different windings njj are automatically taken into account, we can simplify the 
notation and rewrite (4.12) as follows 

i > (q',t')= Jdq' K 0 {q',t'\q,t)i){q,t) . (4.15) 

The wavefunction propagates with the free kernel Kq, but of course it is multi- 
valued. The transformation from ^ to ^ in (4.14) and that from K to Kq in (4.13) 
are obviously the same as those in (4.6) and (4.7) which allowed us to eliminate 
the effective non-local long-range interaction of the Chern-Simons construction, 
and can be performed in general whenever the Lagrangian has the form (4.10). 

The wavefunction in (4.14) can also be written in a different way which 
is more useful for applications. Let us introduce for each particle the complex 
coordinates zj and zj in the following way 


*i = + ix 2 j , 

Zi = x) - ix] , 


(4.16) 


where ( xxj) are the Cartesian coordinates of the J-th particle. It is then easy to 
see that 

ZIJ = Zi - Zj = I Zi - Zj I e 1VlJ , 

zu = zi - zj = |zi - zj\ e~ 1<p,J 
Therefore, (4.14) can be written as 
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* ■(afe) 

= n ( z u ) u , 

i<j 

or equivalently as 

* ■ a.- 

= II(^/j) ” P f (z U ...,Z N ;Z U ...,Z N ) . 

I<J 

In these equations P and P ' denote single-valued functions of the particle positions. 

The generic form of the wavefunction in the anyon gauge is the one presented 

in (4.18) and (4.19), namely a prefactor (either JJ (zuY or n ( zij ) u ) times a 

i<j i<j 

single-valued function. 

It is interesting to observe that these prefactors appear as the obvious general¬ 
ization of the Jastrow-type prefactors that are usually considered for fermions and, 
as we will see in Chapter 9, they can be obtained as expectation values of vertex 
operators in conformal field theory (Fubini 1991; Fubini and Liitken 1991; Stone 
1991a; Dunne et al . 1991a,b; Moore and Read 1991; Cristofano et al . 1991a,b,c). 
Furthermore we point out that the Laughlin wavefunction for the ground state of 
the fractional quantum Hall effect (Laughlin 1983) has precisely the form (4.18) 
with v an odd integer, and that the wavefunction for the excitations above the 
ground state has the same form as in (4.18) with v a rational number. We will 
return to these points later (see Chapter 8) after introducing the second-quantized 
formalism for anyons, and solving the eigenvalue problem (4.4) or (4.8) in some 
simple models. 


(4.19) 


(4.18) 



5. Non-relativistic Chern-Simons Field 

Theory 


In the previous two chapters we have shown that a system described by the Hamil¬ 
tonian 



A,(ri1 ' r " 



2 


5 



where 


A x I {r 1 ,...r N ) = 


2ttk 


E ^ 


r t — r 


J±l 


\ri - rj | 


and by single-valued wavefunctions ip such that 



tp(r i,-- 



■ ,ri, ■■ ■ ,r N ) = tp(r i,- • • ,r/, • • • ,rj, ■ ■ ■ ,r N ) , 



is equivalent to the system described by the free Hamiltonian 




acting on multi-valued wavefunctions ip such that 

• ,rj,- •• ,r N ) = e 1 ™ V>(ri, • • • ,r/, • • •, rj, ■ ■ ■, r N ) , (5.5) 

where 

v = 

The two Hamiltonians (5.1) and (5.4) are simply related to each other by a gauge 
transformation which on the plane allows to remove completely the effective non¬ 
local vector potential Aj(r i, ...r^) 9 . The solution of either one of the two equiv¬ 
alent eigenvalue problems 

H* ip = E ip , (5.7a) 

Hip = Eip , (5.7b) 

is highly non-trivial. In the first case the difficulty is due to the explicit presence 
of a long-range and non-local interaction, while in the second case the complex¬ 
ity arises from the complicated boundary conditions (5.5) that the multi-valued 

9 On a compact surface, like for example a sphere or a torus, one must carefully treat 
zero modes as well as possible topological components of the gauge fields (Iengo and 

Lechner 1990,1991,1992). 


2'kTick 


(5.6) 
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wavefunctions have to satisfy. A complete solution to (5.7) is known only for the 
two-body problem (Leinaas and Myrheim 1977; Wilczek 1982a,b; Jackiw 1990; 
Johnson and Canright 1990; Vercin 1991); since there are no self-interactions the 
one-body problem is trivial and only partial results are available at the moment for 
the iV-body problem with N > 3. In the next chapter we will report on some recent 
progress in finding the solution to (5.7) for any N when external fields are present. 
Here instead, we point out that an equivalent approach to the iV-body problem 
with identical particles is provided by second quantization where the Schroedinger 
eigenvalue equations (5.7) are replaced by a 2+1 (non-relativistic) quantum field 
theory. 

To show this, let us consider a non-relativistic (complex) matter field <f>(t,x) 
of mass m, which for definiteness we take to be bosonic 10 , and let us minimally 
couple it to an abelian gauge field A a with a Chern-Simons kinetic term (see 
Chapter 3). This system is described by the following action (Jackiw 1990; Jackiw 
and Pi 1990; Ezawa et al . 1991a,b; Ezawa and Iwazaki 1991) 




1 



(Di 2 + D 2 2 ) 4> + ^e a ^A a d 0 A, 



where D a = d a + ie A a is the covariant derivative. (To avoid unnecessary redun¬ 
dancy in the formulas, from now on we set c = h = 1, so that in particular x° = t . 
As far as indices are concerned, we keep instead the same conventions of Chapter 
3.) 

Varying the action So with respect to A a , we obtain 

e a ^F^ = — j n (5.9) 


where the current j a is explicitly given by 


j° = <^<!> = p , 

ji = i Di * ~ ( D <^ *) 


(5.10) 


p and j are respectively the number-density and the current-density operators and 
satisfy the continuity equation 


d t p + v-j = 0 . (5.11) 

As is clear from (5.9), the Chern-Simons field strength = dpA^ — d^Ap is 
completely determined by the particle currents. Now we are going to demonstrate 
that also the Chern-Simons potential A a itself is not an independent degree of 
freedom. 

The a = 0 component of (5.9) is 


10 Of course a similar and equivalent discussion can be done for fermionic matter fields. 
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B = --p (5.12) 

K 

where B = V A A = —F 12 is the Chern-Simons magnetic field. This equation 
clearly represents the second-quantized version of the quantum mechanical con¬ 
straint (3.2.22) 11 . In the transverse gauge diA % = 0, we can invert (5.12) without 
ambiguities and solve for the vector potential A . Formally we get 



where G is the Green function for the Laplacian A = V • V, satisfying 


(5.13) 


AG(x - y) = S(x - y) . 

As is well known, the explicit solution to (5.14) is 


G(x -y) = 


27r 


In \x-y\ , 


(5.14) 


(5.15) 


so that A * can be written as follows 



tj 


d 


dxi V 2ttk 



y In \x — y\ p(y) 


2ttk 



d 


y 


(5.16) 


where 

<p(x -y) = tan" 1 (5-17) 

(see also (3.2.30) and (3.2.31)). Since the function (p(x — y ) is multi-valued, partic¬ 
ular care must be used in moving the derivative d/dx % out of the second integral 
in (5.16) and thus displaying A as a gradient. In fact, to eliminate the ambiguities 
deriving from the multi-valuedness of (p(x — y), one must fix a branch-cut in the 
y-plane starting at x . Whatever choice is made for this cut, the resulting range 
of integration of y will depend on ar, and hence extra contributions are produced 
in moving d/dx % outside the y integral. Therefore, it is not correct in general to 
write 

V<p(x — 
d 2 y <p(x 

and consequently it is not true in general that A is a pure gauge and can be re¬ 
moved by a gauge transformation. However there is a special situation where (5.18) 
is certainly correct and A is a pure gauge, namely the case in which the number- 
density operator p(x) is a sum of 5-functions, as is appropriate for a collection 




ii 


Notice that here p is a matter density whereas in Chapter 3, p was a charge density; 
the latter is simply obtained from the matter density upon multiplication by e. 
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of non-relativistic, localized point-particles (Jackiw and Pi 1990). In the quantum 
field theory described by (5.8), we can always choose to remain consistently in such 
a case, provided that we do not consider the model as the non-relativistic limit of 
some relativistic quantum field theory. Indeed it is well-known that in a relativistic 
theory the eigenvalues of p{x) are not sums of 5-functions because particles are 
not point-like but extended (see for instance (Forte 1992a,b)). In this case then, 
particular attention must be paid to the order in which the non-relativistic limit 
and the integrals (5.16) or (5.18) are computed (Dunne et al . 1992b). 

The space components of (5.9) are 

Fio = diAo — doAi = —e *jp , (5.19) 

which together with the continuity equation (5.11), allow to solve for the scalar 
potential Ao yielding 

A 0 (x) = ~ J d 2 y G(x - y) V Aj(y) . (5.20) 

From (5.16) and (5.20), we explicitly see that the Chern-Simons gauge field A a is 
entirely determined by the matter configuration, that is by p and j. 

The use of the Green function G(x — y) in expressions like (5.16) and (5.20) 
must be accompanied by a proper analysis of possible divergent and singular terms. 
For example, it is important to realize that 

(5.21) 

is ill-defined at x = y. Therefore for a complete specification of the theory we 
must supplement (5.14) and (5.15) with a regularization prescription. Every such 
a regularization that preserves the antisymmetry under space reflections leads to 
a vanishing result for (5.21) at x = y. For definiteness it is convenient to use the 
following prescription 


^G<-> - • 


where the regulated Green function G a (x) is 




z\ e 


-z 2 / a 


When the regulator parameter a is removed, we clearly get 

lim G a (x) = G(ar) = In lari . 

a—►O 27T 

On the other hand, it is easy to see that 


(5.22) 


(5.23) 


(5.24) 


lim e^^G a (x) = 0 

X -.0 dxJ v ' 


(5.25) 



47 


for any a. If (5.22) is systematically used, all ambiguities are eliminated and the 
theory is completely specified. 

Let us now quantize the action (5.8) by imposing equal-time commutation 
relations for the bosonic field <j>, namely 


{<j>{x),<j)^{y)\ = 6{x-y) , 
[<t>{x),<t>{y)\ = [<fi{x), <jS{y)\ = 0 


(5.26) 


(For simplicity we dropped the time arguments.) Since the gauge field A is a func¬ 
tion of the number-density operator p = <f>^<f>, the commutation relation between 
A and <f> is not trivial. In fact from (5.26) and (5.16), it is easy to derive that 

[^(x), <j>{y)\ = G ( x - vMv) ■ ( 5 - 27 ) 

As pointed out earlier, one should really use (5.22) so that no ill-defined objects 
at x = y appear. In particular, from (5.25) we conclude that 


[A 1 {x), <j){x )] = 0 


(5.28) 


Therefore there are no ordering ambiguities in the quantum Lagrangian (5.8) or 
in the Hamiltonian arising from it. The latter turns out to be 



n 0 {x) = - l -n'{x)-n{x) 

where II(x) is the momentum operator associated to <£, namely 


(5.29) 


II(x) = [V — ieA(x)\<f>(x) = D<t>{x) , (5.30) 

and iJt (ar) = ( D<f>(x ))* is its hermitian conjugate. For later convenience we observe 
that, after an integration by parts, the Hamiltonian (5.29) can also be written in 
the following way 





d 2 x + ieA * V<f> - V<f>* <f>) + e 2 A 2 


(5.31) 


Let us now derive the first-quantized Schroedinger equation (5.7a) from our 
second-quantized theory. First of all let us observe that the number operator 



(5.32) 


and the Hamiltonian Ho in (5.29) commute. Therefore they can be simultaneously 
diagonalized as follows 

H 0 \E,N) = E\E,N) , 

N\E,N)=N\E,N) . 


(5.33) 
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Furthermore, we assume the existence of a vacuum state |0) annihilated by <f>(x ), 
Ho and N, 

4>( x ) |o) = (o| <f>^(x) = o , 

H 0 |0> = N |0) = 0 , 

but such that <f>^(x) acts non trivially on it. The states for a fixed number N of 
particles in configuration space are 


|ri,---,r w ) = <^ t (r 1 )---<^ t (r w )|0) , 

and consequently the iV-body wavefunctions are 

,r N ) = (ri,---,r N \E,N) = (0| <j>{r x ) • ■ ■ <f>(r N )\E,N) . 

If we project the first equation in (5.33) onto the states (5.35), we get 

( ri ,---,r N \Ho \E,N) = (0\<l>(r 1 )--- ( l>(r N ) H 0 \E,N > 

= (0| [ <t>(r 1 )---<t>(r N ),H 0 ] \E,N) 

= E (0\<f>(r 1 )---<t>(r N )\E,N) 

= E tpirt,---,r N ) . 


(5.35) 


(5.36) 


(5.37) 


The commutator in the second line of (5.37) can be computed without any problem. 
Here we present only the explicit calculation for N = 2, and then infer from it the 
generalization to any N. First of all let us compute the commutator [<j>{ri), Hq\. 
Using (5.31) and the canonical commutation relations (5.26), one easily gets 

[(f>( r i),H 0 ] = — — Ai<f>(ri) + — \A(ri) ■ Vi^rj) + V/ • {A{ri)<j){r /))] 


+ 7T- A ( r i) 2 <t>( r i) 


2m 


+ 


i cP* # 

2m k J 


d 


0x3 


■G(x - rj)(f>(rj) 


X “ ^7^ (x) ^ (x) ) 


+ 


2mn 


J d 2 x 


a 


dx3 


9 

■G(x - r I )4>(r I )A t (x)p(x) 


d 


+A t (x)-Q-r G (x - rj)0(rj)/>(*) 


(5.38) 

where Aj = V/ * V/ = (9/9x})(9/9x}), and the /c-dependent terms arise from 
(5.27). The right hand side of (5.38) can be reorganized with the help of (5.20), 
yielding 

[<t>(ri),H 0 \ = - J—Arfiri) + eA°(ri) <^(r/) 


+ 


2m 
e 


2m/c 2 


— f d 2 x VG(x — rj) * VG(ar — rj) p(x) <f>(rj) . 

IK J 


(5.39) 
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The last term in this equation is a quantum correction arising from the reordering 
of the integrand in the last line of (5.38). 

Since both A 0 and A involve expressions that contain the operator <f >t standing 
on the left and therefore annihilating the vacuum (0|, we immediately conclude 
that 

(0\[cf > (r I ),H o ) = -^-Ai(0\<f>(r I ) . (5.40) 

Let us now return to the calculation of (5.37) for N = 2. Using (5.40) for I = 1 
and (5.39) for I = 2, after some simple algebra, we obtain 

i 2 

(0p(riMp 2 ),i/ o p,2> = - 2 ^£(V/-i eA(n,r 2 )) 2 {Q\<t>{ ri )<t>{r 2 )\E,2) 

(5.41) 

where 

AHn,r 2 ) = ^^G(r/-rj)j . (5.42) 

Using the explicit expression for the Green function given in (5.15), we can rewrite 
this equation and get 





_ r 3 
r I T J 

ri-rj^ 


Therefore for N = 2, (5.37) becomes 


1 



2 

£(V/ -ieA/(ri,r 2 )) 2 

7=1 


V>(ri,r 2 ) = E V>(pi>p 2 ) 


(5.43) 


(5.44) 


This is the first-quantized Schroedinger equation for two particles in the non-local 
potential (5.43). The generalization to any N is now straightforward and the re¬ 
sulting system is described precisely by (5.1). In this way we have explicitly shown 
the equivalence of the field theory formulation and the first-quantized approach to 
the problem of N identical particles. 

It is now interesting to rewrite the first-quantized Hamiltonian (5.1) in complex 
notation. As we explained in Chapter 4, for every particle we can introduce the 
complex coordinates 

zi = x\+ i xj , 
zi = x\- ixj , 

and the corresponding derivative operators 



d = i / a . a \ 

dzj 2 \dx\ 1 dxj) ’ 

= JL = 1(JL + -JL\ 

dzj 2 \dx l j *dxj) 


(5.46) 
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Then we define 




4 ™ 1*1 ZI - ZJ 


Ai = ~\ ( A) + iAj) - ^ * - , 


(5.47) 


4 ™ ~ ZJ - ZJ 


where the second equalities follow from the explicit expression of A/, in such a 
way that the covariant derivative operators in complex notation turn out to be 


Di = di + ieAi = di + 




“ z, - ZJ 


Dj = dj + ieAj = dj — 




4 ™ ~ ZJ - ZJ 


If we recall that 


d - = d r = 7T S^ 2 \z) , 


it is straightforward to check that 


[£>/,£>/] = 5(2) ( 2 -f ~ 


J*i 


Consequently the first-quantized Hamiltonian 


1 N 

H'= -^J2(Vi-ieA I (r 1 ,---,r N )) 


(5.48) 


(5.49) 


(5.50) 


can be rewritten as follows 


H' = V (Dj Dj + Dj Di) 

= --Y i D i D i + — Y^ S^(zj - ZJ) 

mK I?J 

=^ + ^E* (2) (^) • 

I<J 


The total Hamiltonian H 1 differs from the normal ordered term 



Y.DiDj = 


7=1 


_2 

m 



Y,d\d, 


(5.51) 


(5.52) 


by a repulsive hard-core interaction among all pairs of particles. This 5-function 
interaction could be canceled by considering the action 
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S = S ° + 2 f d3x 


(5.53) 


where So is as in (5.8). The action S gives rise to the following quantum Hamilto¬ 
nian 




i7 f (x) • II(x) - | : (<f>1(x) <f>(x)) 2 : 


(5.54) 


Using (5.10) and (5.12), the quartic interaction with coupling constant g that we 
have just introduced, may be written also as 


^ : B(x)p(x) : . (5.55) 

Z6 

Displayed in this form, it is clear that it describes a magnetic-field - charge-density 
interaction. Going through the same procedure as we did before for g = 0, we can 
derive from (5.54) the Schroedinger equation and get the following first-quantized 
Hamiltonian 



_2 

m 


t,D\D I + 


7=1 



(5.56) 


The quartic interaction of the field theory is equivalent to a 5-function potential 
in the quantum mechanical problem. If we fine-tune the coupling g and choose 




ran ’ 


(5.57) 


then the hard-core interaction in (5.56) disappears and the Hamiltonian H f be¬ 
comes simply the normal ordered one as in (5.52). It can be seen that the condition 
(5.57) on the coupling constant g plays a crucial role in studying non-trivial static 
solutions of the classical field theory described by (5.53). In particular it leads 
to the existence of self-dual solitons (Jackiw and Pi 1990; Dunne et al. 1991c; 
Dunne 1992). We mention also that the Hamiltonian (5.52) has been considered 
in (Girvin et al. 1990) as an example of an exactly solvable model for fractional 

statistics and was inspired by supersymmetry considerations (Dj and D\ resemble 
supersymmetry charges) as well as by the use of the Atiyah-Singer index theorem 
for particles in magnetic fields. Moreover, the relevance of (5.54) for the theory of 
the fractional quantum Hall Effect has been discussed in detail in (Ezawa et al. 
1991a,b; Ezawa and Iwazaki 1991). 

Let us now turn again to the second-quantized theory with g = 0. We have 
already remarked that in a non-relativistic context, the eigenvalues of the particle 
number-density operator p are sums of 5-functions, and therefore (5.17) is a cor¬ 
rect formula. This implies that the space components of the gauge field A a are a 
gradient 

A(x) = VA(x) , 

M x ) = J d2 y v( x - y) p(y) ■ 

Now we want to show that also Aq is a pure gauge, namely that 


(5.58) 
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A 0 (x) = -d t A(x) . 


(5.59) 


To show this, we first perform an integration by parts and rewrite (5.20) as follows 


(*) = lj d2 y e ' 3 ^ G ( x ~ • 


(5.60) 


Then, upon using (5.14), we get 


A o( x ) = — 2 ^ J d 2 y v v( x - y) ■ i (y) 

J d 2 y V' <p(x - y) ■ j (y) 

J d 2 y cp(x - y) V' ■ j{y) 


27T/C 

e 


(5.61) 


2ttk 
e 


2ttk 


J d 2 y <p(x - y) d t p{y) = - d t A(x) 


In deriving this result we have used the continuity equation (5.11) and have 
dropped surface terms in the integration by parts because, by assumption, the 
densities are localized. 

Now we can combine (5.58) and (5.61) into a single covariant equation 


■^at(*^) — &cx.A (*e) j 


(5.62) 


which displays the Chern-Simons field as a pure gauge. Therefore it can be removed 
with a gauge transformation of parameter A, namely 


a 


> A' a = A a + d a A = 0 . 


(5.63) 


Under such singular transformation, covariant derivatives turn into ordinary 
derivatives, and the action becomes simply 


Sq = Jd 3 x 2 ^# (di 2 + & 2 2 ) <j> 


(5.64) 


where the new matter fields <f> are defined as 

4>(x) = e~ ieA ( x) <j>{x) , 
<ft{x) = <ft(x) e ieA( - x ^ . 


(5.65) 


_ +* 

The action (5.64) is that of a free (complex) scalar field <£; however such a field 

does not obey ordinary commutation relations. Indeed, using the equal-time com¬ 
mutator 

which follows from (5.26) and (5.58), it is easy to show that for x y 


(5.66) 
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<f>(t,x)4>(t,y) = e ieA( - t ’ x ' ) <f>(t,x)e ieA(t ' y) (i){t,y) 

= e" ie t A ^ x '»+ A ( t ^e il '^y- X )(f>(t,x)cf>(t, y) 

= e - ie I A(t,*)+A(t f y)]^^(y-*)^( tj 3.) (5.57) 

= e -i e [A(t,®) + A(t,y)] e i^(y-*) e i e A(t,y)^ f)y)e i e A(t,®)^ fa . ) 

where 1 / = e 2 /(27r/c) as in (5.6). If y is moved counterclockwise around x, then 


cp(y - x) - cp(x -y) = 7r , (5.68) 

as is clear from Fig. 5.1, and so we have 

4>{t, y) = e i7ri V>(t, y)<j>(t, x) (5.69) 

Thus we have proved that the matter field <f> obeys anyonic commutation relations 
of statistics v . A calculation similar to (5.67) leads to 

x)^(t, y) = y)<j>(t, x) + 5(x - y) . (5.70) 

If x ^ y, we obtain 

= e _1,ri/ <£ t (f,y)<£(f,a;) , (5.71) 

which exhibits again anyonic statistics v\ however we remark that for x = y the 


canonical commutator remains unchanged since the phase proportional to v in this 
case is vanishing. 




y?(a; - y) 


X 


Fig. 5.1. (p(y — a;) — (p(x — y) = 7T . 

The non-relativistic field theory described by the action So in (5.8) or by the 
action S in (5.53) is invariant under the usual space-time symmetries: time trans¬ 
lations, space translations, space rotations and Galilean boosts. To each of these 
symmetries one can associate a conserved Noether current and a corresponding 
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conserved Noether charge according to the standard procedure. We refer the reader 
to (Jackiw and Pi 1990) for a complete treatment of such symmetries and the ex¬ 
plicit derivation of the Noether currents. Here we simply focus on the rotational 
invariance and the associated angular momentum, which turn out to be (see also 
(3.2.51)) 

J = m J d 2 x e t ^x t j^ (5.72) 

where the current j is defined in (5.10). Using the explicit expression of y, we 
obtain 


J = — J d 2 x e tj x % [^(x) (dj + ieAj(x)) <f>(x) - (dj - ieAj(x)) <f>*(x) ^(x)] 
= — Jd 2 x e^x 1 [<f>\x)dj<f)(x) — dj<f)^(x) ^>(x)] + e J d 2 x e t ^x t Aj(x)p(x) 


where p = <f>^<f). 

Upon setting <f> = e l9 p 1 / 2 and <f >* = e~ l9 p 1 / 2 , we see that 


(5.73) 


[(f>Hx)dj(f>(x) - dj<f> ] (x) <f>(x)\ = dj0 p . 


(5.74) 


Therefore, if we assume that the matter fields have a configuration which vanishes 
rapidly enough at spatial infinity 12 , then the first integral in (5.73) vanishes and 
the angular momentum is simply 


J = e 



(5.75) 


Now we have to insert in this formula the expression for Aj given in (5.13) which 
we rewrite here for convenience 


Aj(x) = -A 3 (x) = ~ J d 2 y e jk ^G{x - 

= j d?y e jk^_ G a {x _ y )p ( y) 

The last equality follows from (5.6) and (5.22). Combining everything together we 
obtain 



J = -27 rv J d 2 x J d 2 y p(x)e t3 x t e 3k -^G a (x - y)p(y) 
= 2%v J d 2 x J d 2 y p(x)x l -^rG a (x - y)p(y) 

Let us evaluate this expression for the matter density given by 


a —>0 


(5.77) 


12 This is certainly the case for non-relativistic configurations, where p is a sum of 8- 
functions. 
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N 

p( x ) = 5Z^ (2) (* - r dt)) ■ (5.78) 

1=1 

As we have already remarked, this is equivalent to consider N point-like non- 
relativistic particles located at r/(i). Inserting (5.78) into (5.77), we have 


J = 2,7ns 


N r r n 

Y J d2x J d2y ^ (2) (* - r i( t )) xt fai Ga ( x - y) s(2) (y - rj(t)) 


a —*0 


N 


271 - 1 / Y, r i-£i Ga ( r i - r i) 
ij =i dr 


a —*0 


N 


= 2 ™YY, r i 

1=1 J^I 


: d 


dr 


jG a {ri - rj) 


a —*0 


(5.79) 

where the last step is due to (5.25), that is to the absence of self-interactions in the 
regulated theory. We can proceed even further if we observe that the regulator a 
can now be removed from (5.79), so that using the explicit expression for G given 
in (5.15), we obtain 


N 


J = 


i=i jjti 1 


-rj) 


a —>0 


N 


7=1 J^I 


r r — r 


rj - rj 


j _ 


v-N{N 


- 1 ) 


(5.80) 


This result coincides with the one we have derived in Chapter 3 using the first- 
quantized theory (cf. (3.2.53)). Notice that for point-like particles there are no self 
interactions and the angular momentum J comprises only the orbital part, as is 
clear from the fact that (5.80) vanishes for N = 1. We can therefore conclude that 
in a non-relativistic theory with Chern-Simons gauge fields, there is a fractional 
orbital angular momentum and particles carry no spin. 

The situation is drastically different in a relativistic field theory (Forte and 
Jolicoeur 1991; Forte 1992a,b) or even when the model emerges from a relativistic 
theory in the non-relativistic limit. In this latter case, one should consider extended 
particles in the limit as their size approaches zero; then a delicate question arises 
about the order in which the size of the particles and the regulator parameter a 
are taken to zero. If we perform the zero-size limit first, and then compute J and 
remove the regulator at the end, we clearly obtain the result (5.80); however if we 
proceed in the reverse order and first compute J for extended densities and then 
take the limit of zero-size, we get a different result (Dunne et aL 1992b) . Now we 
show explicitly how this may happen. 

Let us consider a matter density of the form 
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N 

p( x ) = pi ( x ) ( 5 - 81 ) 

1=1 

where pi(x) is a spherically symmetric function centered around the point r/(t), 
which in the limit of point-like particles eventually reduces to a 6-function. For 
example we could take 


Pl (x) = 

U 7T 


(5.82) 


such that 


d 2 x pi(x) = 1 , 

lim pi(x) = 8^ 2 \x — r/(t)) 


(5.83) 


The parameter b can be considered as related to the size of the particle, and will 
be set to zero at the very end of the calculation. Upon inserting (5.81) into (5.77), 
we get 


J = 2'KV 


n r r 8 

Y, d 2 x d 2 y p I (x)x t — i G a {x - y)pj(y) 
7,7=1 d J ax 


a —*0 


b—>0 


(5.84) 


= Lu + S v , 


where 


L u = J d 2 x J d 2 y p I {x)x t ^- i G a {x-y)pj{y) 

N r r f) 

YlJ d2x J d2y pi ( x ) xt faiG a ( x ~ v)pi(y) 


a —>0 


b—>0 


S u = 2'KV 


(5.85) 


6-*0 


The reason for this decomposition and the choice of the symbols will become ap¬ 
parent in a moment. The calculation of L u in (5.85) presents no difficulty, yielding 
as before 

L v = v)-N{N - 1) . (5.86) 

Comparing this with (5.80), we can conclude that in this calculation the order in 
which the limits a —> 0 and b —> 0 are taken, does not influence the result. However 
the situation is different for S u in (5.85). To see this, first of all let us observe that 
using the obvious property 

_ „) = - V) , (5.87) 


S u can be rewritten as follows 
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/=1 J 



d 


d l x j dry pi{x){x % - y % )—Xi a {x - y)pi{y) 


a —*0 


. (5.88) 


6->0 


Then, using (5.82) and shifting the integration variables, we obtain 



o 


o 



(5.89) 

We stress once again that if we had computed (5.88) by taking first the limit b —> 0 
and then the limit a —> 0, we would have gotten a zero result. 

It is now clear what is the meaning of L u in (5.86) and of S u in (5.89). L u is 
a fractional orbital angular momentum, whereas S u is a fractional spin. which 
is the second quantized version of the quantum mechanical angular momentum 
l u (see (3.2.53)), originates from the interaction among different particles and is 
proportional to the number of pairs, N(N — l)/2. On the contrary S*,, which 
does not vanish for N = 1, comes from self-interactions of extended anyons and is 
proportional to the number of particles, N . The combination 


J = L v + S v = \v N(N - 1) + \v N (5.90a) 

is often presented in the literature as 

J=\v N 2 , (5.90b) 

and is regarded as a manifestation of a non-standard composition rule for the 
angular momentum of particles with fractional statistics. Indeed, if one anyon has 
angular momentum v j 2 - in agreement with the spin-statistics connection - N 
anyons have a total angular momentum N 2 v/ 2. However we should realize that 
in N 2 v/ 2 there is an orbital contribution and a true spin contribution, and if we 
keep this in mind, (5.90b) does not come as a surprise at all. 

We can summarize our results as follows. For N point-like particles of statistics 
v (see Fig. 5.2), the total angular momentum is 

J = L V = \v N(N-l) , (5.91) 

z 


and no spin is present. 

For N extended anyons in the point-like limit (see Fig. 5.3), the total angular 
momentum is 

J = L v + S v = \v N(N - 1) + \v N 

£ £ 

where S u is the spin and the standard spin-statistics relation is verified. 


(5.92) 
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Fig. 5.2. Point-like anyons can be represented as particles with infinitesimally thin flux 
tubes. The fractional orbital angular momentum L u originates from the interaction be¬ 
tween the flux tubes of different particles. 



Fig. 5.3. Extended anyons can be represented as particles with extended flux tubes. In 
this case self-intersections occur and give rise to the spin Si/. 



6. Anyons in a Magnetic Field 


The general theory of fractional statistics we have presented in the previous chap¬ 
ters, shows that anyons provide a fascinating model in which to study topological 
effects in quantum mechanics and in quantum field theory. To complete the anal¬ 
ysis and see whether or not anyons are relevant for physics, one needs to develop 
the many-body theory, and then study the thermodynamic properties of systems 
of very large numbers of anyons as we customarily do for many-boson and many- 
fermion systems. However a deep understanding of the many-body problem for 
anyons is still lacking, despite the significant progress that has been recently made 
in various models using analytical, approximate and numerical methods. 

It has been known since the early days of fractional statistics that the two- 
body anyon problem is completely soluble (Leinaas and Myrheim 1977; Wilczek 
1982a,b; Wu 1984a) and that in the three-anyon problem with a harmonic potential 
a set of exact wavefunctions can be easily found by solving a partial differential 
equation (Wu 1984b). The latter approach has recently been extended to the N- 
anyon problem in the presence of a harmonic potential (Chou 1991a). On the 
other hand, the (related) problem of anyons in an external magnetic field has been 
analyzed by M.D. Johnson and G.S. Canright, who gave detailed results for the 
two-body problem as well as very partial results for the iV-body problem (Johnson 
and Canright 1990). The external magnetic field breaks the free particle degeneracy 
and reveals a spectral structure reminiscent of the familiar Landau levels of charged 
particles in an external magnetic field (see for example (Landau and Lifschitz 
1977)). This iV-body problem has recently been studied in more detail (Dunne et 
al. 1991a,b) in connection with the surprising relationship between anyonic theories 
and conformal field theories which has been pointed out by many authors (Fubini 
1991; Stone 1991a; Fubini and Liitken 1991; Moore and Read 1991; Cristofano 
et al. 1991a,b,c). Finally, we mention that numerical studies have recently been 
performed in the three- and four-anyon problems with a harmonic potential (Sporre 
et al. 1991a; Murphy et al. 1991; Sporre et aL 1991b), and that perturbative 
calculations for statistics infinitesimally close to the fermionic or bosonic ones are 
now available for many models (Chou 1991b; Me Cabe and Ouvry 1991; Comtet 
et al. 1991; Karlhede and Westerberg 1991; Khare and Me Cabe 1991; Dasnieres 
de Veigy and Ouvry 1991; Sporre et al. 1991c; Chou et al. 1992). 

In this chapter, I am going to present a systematic discussion of the iV-body 
problem for anyons in the presence of an external magnetic field, following closely 
the article I wrote in collaboration with G. Dunne, S. Sciuto and C. TVugenberger 
(Dunne et al. 1992a). The study of anyons in magnetic fields is motivated by 
the following considerations. The only known example of objects which can be 
described as anyons are the quasi-particles aiid the quasi-holes of the fractional 
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quantum Hall effect (see for instance (Prange and Girvin 1990) and Chapter 8). 
These have the important feature of feeling the underlying electrons exactly as 
a charged particle would feel an external uniform magnetic field. Thus in their 
physical realizations, anyons experience effective magnetic fields (not to mention 
the strong external magnetic field which drives the Hall effect itself). To find the 
solution to the problem of anyons in magnetic fields is therefore interesting and 
important also for practical purposes (for related works see also (Grundberg et al 
1991b; Karlhede and Westerberg 1991; Cho and Rim 1992). 


6.1 The Eigenvalue Problem 


The Hamiltonian for a collection of N identical anyons (labeled by capital Latin 
indices /, J,...) of statistics v , mass m and electric charge e, moving on a plane 
in the presence of an external constant, uniform magnetic field B (perpendicular 
to the plane) is given by 13 

1 N 

H* = - 2m — i4z‘( r i»— r Jv) -ieA(r/)) 2 . (6.1.1) 

rn i=i 

We choose the external vector potential A ( rj ) in the symmetric gauge, namely 

A i (r) = - |e”V , (6.1.2) 

and without any loss of generality we assume B to be positive. The non-local 
effective gauge potential 




(6.1.3) 


describes the statistical interactions (see (3.2.27)). The associated statistical mag¬ 
netic field Bj is 


Bi (t*i,... ,r^r) = V/ A A/(ri, ...r N ) = -2nv S 2 (r 7 - rj) , (6.1.4) 


so that each particle sees the (N — 1) others as vortices of strength (—27ri'), 
and a statistical phase itv under interchange of two particles is acquired via the 
Aharonov-Bohm effect. Note that since the statistical phase is 7r v we may restrict 
v to values in the interval [0,2). 

The non-local potential Aj is pure gauge except at the locations of the particles, 



<Pl J = 



(6.1.5) 


13 In this chapter we set h = c = 1. 
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and therefore can be removed from (6.1.1) by going to the anyon gauge (see Chap¬ 
ter 4). The resulting Hamiltonian is 




( 6 . 1 . 6 ) 


7=1 


which acts on the redefined wavefunctions 


iv X) Vu 

ip(ri,...,r N ) = e ’< J ip (ri,..., r N ) 


(6.1.7) 


Thus in the anyon gauge, the Hamiltonian reduces to the usual sum of one-body 
Hamiltonians, the statistical interactions disappear and become entirely hidden 
in the boundary conditions which make the wavefunctions 'ip multi-valued . Conse¬ 
quently, the many-body Hilbert space is not a simple direct product of one-body 
Hilbert spaces. Throughout the rest of this chapter we shall work with these multi- 

a * a * 

valued wavefunctions, and for simplicity we drop the * sign from 'ip and H. 

Using the complex notation 


z = x 1 +ix 2 




the Hamiltonian becomes 




—d/d/ + 

m 




( 6 . 1 . 8 ) 


(6.1.9) 


where J c is the total (canonical) angular momentum operator 



Jc = £ (2/d/ - 2/d/) 
i=i 


( 6 . 1 . 10 ) 


We note that J c commutes with H so that we can solve the following eigenvalue 
equations 


H'ip = E'ip 

Jc Ip = J Ip 


( 6 . 1 . 11 ) 


For N = 1 these equations form the Landau problem whose solution is known since 
the early days of quantum mechanics (Landau 1930). Let us now briefly review 
some elementary facts concerning this solution. 

The properly normalized one-body Landau wavefunctions are 





nl 


3 T 3 

7T2J+ 1 (n+7)! 2 n 


eB 


eB 


exp - 


( 6 . 1 . 12 ) 


where n and j are integers such that n > 0 and j > —n. They satisfy (6.1.11) with 
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E = to + - 

j=j ■ 


(6.1.13) 


where u = eB/m is the cyclotron frequency and the U n are generalized Laguerre 
polynomials. We refer the reader to Appendix A where we list some properties 
of the Laguerre polynomials that are useful for our discussion. Here we remark 
in particular that when j is negative, which can happen only for excited states 

(n > 1 ), 

. _/ eB. ( eB. 




3L n+j 


(6.1.14) 


Thus an antianalytic prefactor can appear only for excited states. The energy spec¬ 
trum is structured in discrete levels similar to those of the harmonic oscillator and 
known as Landau levels (see for example (Landau and Lifschitz 1977)). These lev¬ 
els have an infinite degeneracy given by all states with the same quantum number 
n and arbitrary angular momentum j > — n. However for samples of finite area 
this degeneracy is finite since in this case the angular momentum is bounded from 
above. 

Given the wavefunctions (6.1.12), the solution of the iV-body problem is im¬ 
mediate for bosons or fermions: in fact, all many-body wavefunctions are simply 
obtained by forming respectively all totally symmetric or antisymmetric products 
of N one-body wavefunctions. In particular, it is easy to show that the ground 
state for N particles of ordinary statistics has energy E and angular momentum 
J given by 


bosons: 


E = (o 


N 


J — d 


4* j 


fermions: 


E = uj 


N 
2 * 


(6.1.15) 


J=-N(N-l) + d + , 


where is a non-negative integer. These are of course the expected results. 

Let us now turn to the iV-body problem for particles of fractional statis¬ 
tics. It is convenient to extract from the wavefunction in (6.1.11) the factor 


exp (-f E |2/| 2 

V 7=1 


so that the eigenvalue problem reduces to 


N 


Hip = \ E-u-Jil> , 

A» A» 

J c ip = Jip , 


(6.1.16a) 

(6.1.16b) 


A A 

where the new Hamiltonian H and the new wavefunction ^ are defined by 


H = -—didi + —2/d/ , 

' m m 

/=i L J 


m 


$ = ex P ( S l 2/ l 


(6.1.17) 


•4) . 
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^ A 

The wavefunction $ must satisfy some physical requirements. First of all, ^ 
must vanish at points of coincidence if v ^ 0 (cf also (Grundberg et aL 1991a)). 
This hard-core requirement excludes the possibility of overlapping particles: indeed 
as we have seen in Chapter 2, the possibility of intermediate statistics in two 
space dimensions is due to the non-trivial topology of the configuration space 
of N identical particles with the coincidence points excluded and with permuted 
configurations identified, whose fundamental group is the braid group Bn (see 

A» 

(2.11) and (2.12)). According to our construction the wavefunction ^ must satisfy 
the multi-valued boundary conditions (6.1.7) which indeed characterize an Abelian 
representation of Bn- 

A very general form of ^ with the above properties is 

= n (zi j) aiJ (zu) blJ <f> (zi ,--.,z n ;z u ...,zn) + symm (6.1.18) 

i<J 

where zjj = (zj — zj ), “symm” means symmetrization with respect to permu¬ 
tations of the index I in the variables 2 /, 2 /, and the parameters ajj and bjj 
satisfy 

a u ~ hj = v [mod 2] , a/j + bjj > 0 , (6.1.19) 

for any pair (/J). Finally, <f> is a single-valued function which we choose to be 
regular and non-vanishing at all coincidence points. 

A few remarks are in order. First of all, we observe that contrary to naive expec¬ 
tations, symmetrization does not automatically produce a bosonic wavefunction 
for v ^ 0,1 as one can see from the simplest example 

(z 12 r + symm = (z 12 y + (z 21 y = {1 + e™") (z 12 r . (6.1.20) 

Secondly, we want to emphasize that the restriction on the function <f> in (6.1.18) 
to be regular and non-vanishing, amounts to the requirement that <f) be expandable 
as Taylor series in zjj and zjj for all pairs (/J), and that the leading term in this 
expansion be non vanishing. Moreover, since <f> is by assumption single-valued, it 
is clear what is the meaning of the two conditions (6.1.19): the first one enforces 

A» 

the right statistics on the second excludes possible divergences (notice however 
that ajj and bjj are not necessarily positive). 

We want to stress that (6.1.18) is not the most general form of a wavefunction 
that is not singular at coincident points and has the right symmetry properties un¬ 
der exchanges. In fact, all functions with unusual angular dependence are excluded 
from our analysis. For example for N = 3 the function 

(^12 223 213 )" V \ z 12\ 2 + 1223 \ 2 + 1^13 | 2 ( 6 . 1 . 21 ) 

has the right braid group properties, is not singular for zjj —> 0 but is not of 
the form (6.1.18) because the square root term is not expandable as Taylor series 
in zij \ Clearly (6.1.21) is not a solution to our eigenvalue problem (6.1.16a) but 
is an example of those acceptable wavefunctions which are left out from our dis¬ 
cussion and which could give rise to missing states. At the present time it is not 
known how to solve the eigenvalue problem for N > 3 anyons in its full generality. 
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Indeed, the exact analysis we are going to present in the following, based only 
on wavefunctions of the form (6.1.18), is the most extensive one available using 
only analytic methods. Further considerations about the missing states and the 
full solution of the many anyon problem can be obtained using semiclassical (Il¬ 
luminati et al. 1992; Sporre et al. 1991c), numerical (Sporre et al. 1991a; Murphy 
et al. 1991; Sporre et al. 1991b; Chin and Hu 1992) or perturbative (Chou 1991b; 
McCabe and Ouvry 1991; Comtet et al. 1991; Khare and McCabe 1991; Dasnieres 
de Veigy and Ouvry 1991; Chou et al. 1992) methods as we will see in Section 4, 
Let us now return to our exact analysis. Under our assumptions, we are now 
going to show that for v ^ 0,1 either all bjj or all ajj vanish, if one requires that 

A A _ A A. 

ip be an eigenfunction of the Hamiltonian H. The most singular term in H^ comes 
from the Laplacian ^ djdj and is explicitly given by 

I 


2 ^ 

m 


E 



I<J 


,Z N ]Z X ,... >zn) -l-symm 


( 6 . 1 . 22 ) 

where for I > J we define ajj = ajj and bjj = bjj. Of course to satisfy the 
eigenvalue equation (6.1.16a), this term must cancel. Let us suppose that the set 
A of non-integer aj/s is not empty and let us call a a p an element of A such that 


fla/3 < fl/J ) Va/J € A 


(6.1.23) 


A moment thought reveals that because of (6.1.23), symmetrization cannot pro¬ 
duce accidental cancellations and therefore (6.1.16a) implies that the expression 
in square brackets in (6.1.22) must vanish, that is 


b<xi — bpj =0 VJ^a , VJ^/3 . 

Upon inserting this into (6,1.19), one gets 

a a i = v [mod 2] , apj — v [mod 2] . 

The argument can be repeated, until proving that all bjj vanish and all ajj = 
v [mod 2]. 

Since the regularity of the wavefunction at coincident points excludes the pos¬ 
sibility that some ajj be equal to v — 2, v — 4 and so on for v € [0,2), we find that 
the wavefunction (6.1.18) is actually of the form 

fa = W{zijY P{z x ,...,z N ]z u ...z N ) , (6.1.24a) 

I<J 

where possible extra powers coming from ajj = i^+2, i/+4,.., , have been absorbed 
in P and the same has been done for symmetrization. Hence P is completely 
symmetric under any permutation of the particle labels. 

If one starts by supposing that some bjj is not an integer, the proof is exactly 

A» 

as above, and one finds that ^ must have the following expression 
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^II(^) 2 1 ' p ( z Ii'"'> z n\zU"^ z n) , (6.1.24b) 

I<J 

where, once again, the regularity for zjj —> 0 has been required. We will call 
the states in (6.1.24a) and (6.1.24b) Type-I and Type-II states respectively to 
emphasize the different structure of their multivalued prefactors. 

We remark that for v ^ 0,1 the states t/>i and do not overlap and form two 
distinct families. In fact the radial behavior of the function P in (6.1.23) at points 
of coincidence is \zjj\ e with l even and non-negative. This prevents transforming 
Type-I states into Type-II states by writing for instance 

P= H(zu ) 2 \zu\- 2v P • (6-1-25) 

I<J 


A A 

On the contrary xfti and V'ii states do overlap for v = 0,1. It is important to realize 
that for v — 0, all bosonic wavefunctions can be written as in (6.1.23) with P a 
polynomial in 2 /, 2 /. Also all fermionic (z/ = 1) wavefunctions can be written in 
the representation (6.1.23): in this case however P is not necessarily a polynomial 
in 2 /, 2 /, as is clear from the following example. Consider three fermions described 
by the wavefunction 

^ = z \2 z 3 + ^23^1 + Z 31 z 2 * (6.1.26) 

This is an eigenfunction both of energy and angular momentum, that it is com¬ 
pletely antisymmetric under particle exchanges and that it can be written as in 
(6.1.24a) with 

P = —— + —— + — — , (6.1.27) 

^13^23 ^12^13 ^12^32 

which is finite at all coincidence points but is not a polynomial 14 . The special 
case in which P is a polynomial also for fermions is when all particles are in the 
first Landau level (Dunne et al. 1991a), 

In the next sections we are going to study in detail the eigenvalue problem 
(6.1.17) with wavefunctions ^ of the form (6.1.23). However already at this point, 
it is clear that the simplest of all such wavefunctions are obtained by setting P = 1, 
namely 

$ 0) = II (*")" (6.1.28a) 

I<J 

for Type-I, and 

tn = II (^) 2_l/ (6.1.28b) 

I<J 

for Type-II. They are energy and angular momentum eigenstates with eigenvalues 



Jj (0) = ^-N(N - 1) , 


(6.1.29a) 


14 Moreover, its value in 212 —> 0 for instance, depends on arg( 2 i 2 ). 
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and 

4°’ = » (^N(N - 1) + f) 

-C = ~N(N - 1) , 

respectively. Therefore for v ^ 0,1 and form two distinct and preferred 
“base” states upon which one can build two families (Type-I and Type-II) of 
regular functions as we will see in the next section. 


(6.1.29b) 


6.2 Step Operators and Wavefunctions 


A possible strategy to construct exact many-body wavefunctions for anyons in a 
magnetic field is to start from the “base” states 


for Type-I, and 


tr =n <*»>" 

(6.2.1) 

I<J 



(6.2.2) 


I<J 


for Type-II, and act on them step operators which increase the energy or the 
angular momentum. For simplicity and to avoid repetitions, we limit our discussion 
to Type-I wavefunctions; Type-II wavefunctions can be treated in the same way 
with zj —> 2 /, v —> 2 — v. 

To construct the appropriate step operators, let us observe that the Hamil- 
tonian H in (6.2.17) and the angular momentum J c in (6.2.10) can be written 
as 


N 


H = u> ^ a\aj , 

7=1 

(6.2.3a) 

N 

j c =y: _ a \ a i ) » 

7=1 V 

(6.2.3b) 

where the operators 


*■?-+ 

II 

1 

to M 

0 

II 

S* 

(6.2.4a) 

b\ = zi ^di , 6/ = di , 

(6.2.4b) 

satisfy [a/,a^] = [ 67 , = Sjj , all other possible commutators being zero. 

From (6.2.3) and (6.2.4) we see that the a\ and b\ act as “creation” opera¬ 
tors while the a/ and bj act as “annihilation” operators (this will be made more 
precise below). Therefore our general strategy will be to construct multi-anyon 
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wavefunctions by acting with these creation operators on the base states in (6.2.1) 
and (6.2.2). We must, however, respect the statistics and hard-core requirements 
for the resulting wavefunctions as discussed in the previous section. To respect the 
statistics we may only use symmetric (in the particle labels I) combinations of the 
step operators — the statistics being encoded entirely in the base states (6.2.1) and 
(6.2.2). To respect the hard-core requirement we ask that wavefunctions vanish at 
coincident points (for v ^ 0). Therefore we will have reject those combinations of 
step operators which produce singular (at coincident points) wavefunctions. 

The reason we called aj annihilation operator (for Type-I states) is that it 
vanishes on the “base” state, namely 

a/ $ 0) = 0 (6.2.5) 

for all /. Moreover, from (6.2.3a) it is immediate to see that 

H $ 0) = 0 . (6.2.6) 

Also bj is somehow an annihilation operator, but its action on the base state is 
not so simple; indeed one finds 



N 

However, the symmetric combination bj does annihilate 

7=1 



i.e. 


(6.2.7) 


N 

Y2 bi $ 0) = 0 , (6.2.8) 

7=1 

and this justifies our assertion. When v ^ 0,1 (i.e. for generic anyons), higher 
powers of bj (as well as products of a\ and bj) produce singularities when acting 

^ /q\ 

on 'tpl , which cannot be removed by symmetrization. Thus these operators cannot 

be used and only the “creation” modes (a* and h|) have to be considered. When 
v is integer (v = 0,1) the situation is different. In fact, for v = 0 the operator 

(h/) n annihilates and for v = 1 it produces regular functions when acting on 

This property is what singles out the integer (0 or 1) values of the statistics 
parameter v. In conclusion, to produce acceptable wavefunctions from the base 
state (6.2.1) for non-integer i/, we can use only the creation operators a\ and b\, 
while for i/ = lwe may also use bj 15 . In particular one may show that, for v = 1, 
symmetric combinations of a*, b\ and bj acting on (6.2.1) produce all fermionic 

15 Note that also the combinations b\ n bi (appropriately symmetrized) produce regular 
wavefunctions when acting on However, these have the form of polynomials mul¬ 
tiplying and can also be obtained by the action of &J 5 s alone. Therefore the com¬ 
binations b\ n bi are redundant. 
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states. We can give an easy example of this fact. Let us consider the lowest angular 
momentum state with two fermions in the ground level and one in the first excited 
level which is described by the wavefunction (6.1.26). Some very simple algebra 
reveals that 

1 3 

$ = 212^3 + 223^1 + Z 31 Z 2 = -^2 a\b] $ 0) . 

1 7=1 

Notice that there is no analogue state involving the bf s for v ^ 1, since the 
corresponding wavefunction would be singular. 

To treat the problem systematically, let us recall that a basis in the ring of 
symmetric polynomials in 2 N variables ( a\ and b\ with / = 1, ...,iV) is given by 

N 

C tn = Y,A% n . (6-2.9) 

7=1 


where l, n are non-negative integers such that £+n < N . In particular the operators 

N N 

Cio = ^2 cl] and Coi = ^ b] describe completely the center of mass excitations 
7=1 7=1 


and were originally introduced in (Johnson and Canright 1990). More generally, 
the operators C^ n , which respect the statistics properties of the base state 


are step operators in energy and angular momentum, in fact 


[H,C en ] =iuC tn , 

[ j, c en ] = (n - i) C tn . 


(6.2.10a) 

(6.2.10b) 


Therefore Ct n can be used on the base state to produce exact energy and angular 
momentum anyonic eigenstates (Dunne et al 1992a; Grundberg et aL 1991b). 
However, since in general they contain derivatives, we are not guaranteed that this 
procedure produces regular eigenstates and particular care must be used to identify 
which particular C^ n ’s produce acceptable states. In any case it is clear that our 
procedure will lead to states of the form (6.1.24a) with P(^i,... zn, Zi, ... zn) 
being a polynomial function. Thus, for v ^ 1 it is not guaranteed that the C^ n ’s 
produce all anyonic eigenstates. 


To select the “safe” operators, let us begin with Co n 


N 


12 b\ n (with n < N). 


7=1 


Using (6.2.4), we obtain 



( 6 . 2 . 11 ) 


which is regular at all coincident points. This shows that all operators Con’s with 
n < N can be safely applied to the base state to generate regular wavefunctions, 
which have the following form 


N 

^({A0n})=n( C '0n) A °"^l (0) 

n=0 


( 6 . 2 . 12 ) 
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for arbitrary non-negative integers Ao n - 

From (6.2.9a) and (6.2.6) it is immediate to realize that the functions (6.2.11) 
are annihilated by the Hamiltonian H so that their energy is E = Nu/2. They 
correspond to configurations in which all particles are in the lowest Landau level 
with total angular momentum 


Won}) 


N(N-l) 
V 2 


N 

d” ^ } 7lA()n • 

n=0 


(6.2.13) 


In order to go to excited energy levels, we have to consider necessarily operators 
Ce n with l > 1; since these objects contain derivatives with respect to 2 /, they 
can lead to singular wavefunctions and a careful study is needed to select those 
combinations which preserve regularity. In this respect an asymmetry between 
0 < v < 1 and 1 < v < 2 seems to appear: one can surely apply one derivative dj 

to for 1 < v < 2 without getting a singular function for zjj —> 0; the same 
seems not to be true for 0 < v < 1. On the contrary, we will show that actually 

a/A\ 

there is no asymmetry and that C\m ( with m = 0,1,...) applied to * still gives 
a regular function for zjj —> 0 for 0 < v < 2. In fact, from (6.2.6) and (6.2.4), 

C'lm^ I (0) = ^2 {Zl ~ 9l) (Zl ~ Bl) m $ 0) 

/ 

= ^2(zi- 9/) z?4 0) (6.2.14) 

I 

= - E z ? d i$ 0) +E fr z ? - $ 0) > 

/ / 


where for simplicity we have set eB = 2. Only the first term is potentially singular; 
however 



(6.2.15) 


a/A\ * 

is regular because the factor in front of is a polynomial. 

Actually, Ci m can be repeatedly applied to (or to ^({Ao n })) without 
generating singularities. In fact, let us consider the state 





5 


and write Ci m; as 




= 53 (5/ - 9i) (;zi - a/) mi . 

7=1 


a/a\ 1 ^ 

The destruction operators aj either annihilate ; or cancel one of the a} in the 
C\m . (j > i). Therefore the potentially singular terms are again produced by 
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N 

S z ?9j (0 < n < mi) for which (6.2.15) holds. Therefore we have enlarged our 

7=1 

family of allowed states to include also 


N-i n 

in ({Ai m , A 0 ^}) = 

m=0 7=0 

= II II (Cnm) K ” II 

n=0 m=0 I<J 


(6.2.16a) 


where A nm are non-negative integers, with energy and angular momentum given 

by 


N-1 


E — UJ \ Aim + 


N 


(6.2.16b) 


m=0 


1 N-n 


j N(N-l) 

j = —- + J2 _ n ) A 


nm 


(6.2.16c) 


n=0 m=0 


The family (6.2.16a) can be further enlarged by observing that some of the 
states (6.2,12) are of the form 


ti a) = n izuT +2s , 


* = 1 , 2 , 


I<J 


Thus a string of operators C nimi C n 2 m 2 ... C nirni with £ n i < 2s can be applied 

3 = 1 

A / ^\ 

to without generating singularities, because such a string contains at most 

t 

derivatives of order 71 j with respect to zj. 

3=1 

Moreover, the same argument of (6.2.15) can be easily extended to prove that 

N 

X) Zj*dj +1 applied to JJ ( zjj) p gives fj ( zjj) p ~ 2k times a polynomial (notice 
7=1 I<J I<J 

that the exponent of f] {zuY is lowered by 2 k and not by 2 k + 1 as one may 

I<J 

naively expect). By the same reasoning that led to the states (6.2.16a), we can 
write the following large family of wavefunctions 

N N-n 

^i({A„ m }) = JJ n (Cnm) Anm J{{zij) v+2a , * = 0,1,2,... (6.2.17) 

n=0 m=0 I<J 

where the non-negative integers A nm are constrained by 

N N-n 

{ n }Anm — 2s 

n=0 m=0 


(6.2.18) 
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where {n} is the largest even number not greater than n. The energy and angular 
momentum of the states (6.2.17) are 

/JL N ^ N \ 

B = w EE nXnm + t 

\n= 1 m =0 ^ ) 

1\I( 1ST 1 \ -W N—n 

J =(v + 2s) - y^A m ~ n )^nm 

n =0 m =0 

In the constraint (6.2,18) we do not allow the left-hand side to be less than 2s 
because the corresponding states are already contained in (6,2.17) with s replaced 
by s — 1, s — 2,_We remark that the states (6.2,17) are not all linearly indepen¬ 

dent because of highly non-trivia! polynomial relations between products of the 
C nm ’s (Dunne et al. 1992a), It is also interesting to note that the left-hand side 
of (6.2.18) does not depend on n, but only on {n} and therefore it jumps by two 
steps; this confirms that there is no asymmetry between the cases 0 < v < 1 and 
1 < v < 2 . 

A similar discussion can be repeated starting from leading to the Type-II 
states 


(6.2.19a) 

(6.2.19b) 


N N—n 

^Il({Anm}) = II II (Cnrn) A " m , S = 0 , 1 , 2 , . . . ( 6 . 2 . 20 ) 

n=0 m=0 /< J 

where the non-negative integers A nm are constrained by 

N N-n 

= 2 s . (6.2.21) 

n=0 m=0 

The energy and angular momentum are given by 

£ = o,(f:X:"nA nm + ^ + (2-‘. + 2 ») Ar(W ~ 1> ) , (6.2.22a) 

\n=l m=0 / 

7V7Y at -I \ N N—n 

j = -(2 - i/ + 2 s) 3 } E ( m - n ) A ™ • ( 6 . 2 . 22 b) 

n=0 m=0 

So far we have considered symmetric polynomials in the operators a\ and b\ 
constructed from the symmetric step operators C nm . However there is another 
way of obtaining symmetric polynomials in a\ and namely multiplying an even 
number of antisymmetric polynomials. The latter can be written as determinants 
in a\ and b \, but as we have already stressed, not all such expressions produce 
regular wavefunctions and are allowed. A detailed analysis is therefore necessary to 
select “safe” operators, and in this way a further set of regular wavefunctions can 
be constructed. We refer the reader to the original literature (Dunne et al. 1992a) 
for a complete discussion of this procedure. Here instead we work out explicitly the 
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two-anyon case starting from the general states (6.2,17) and (6.2.20). Eq. (6.2.17) 
for s = 0 gives the wavefunctions 


V-i ({A nm }) = (C' 1 o) Alo (Co 1 ) Aoi (C 11 ) A "(C 02 ) A -( 2 i2) l/ • (6.2.23) 


The center of mass dependence can be easily factorized by defining, 


= a{ + <4 , 

a* = a | - <4 , 


B'=b\+b\ , 
b' = b\-b\ , 


(6.2.24) 


so that we can write 

C 10 =At ; C 0 i=B' ; C u = \ {A'b' + a+fet) ; Co2 = \ {B* 2 + b' 2 ) . 

(6.2.25) 

Using this notation, one can write a more convenient basis in the space spanned 
by the states (6.2.23), namely 


fa(N,M,n,m) = A' N B' M a' n b' m {z 12 f (6.2.26) 


where iV, M, n, m are non-negative integers, and m — n = 0,2,4,... 

Eq. (6.2.26) reproduces half of the well-known eigenfunctions for two anyons 
in a constant magnetic field. Indeed, after some elementary algebra we find 


A 

(iV, M, n, m) oc Z 


M-NtM-N 
N 


. _.2 \ _ \m—n+v j^m—n+v ( I .. |2 


l^l 2 M 


1*12 I 


(6.2.27) 

where Z = ~^(zi + 22 ) is proportional to the center of mass coordinate. The other 
eigenfunctions come from the Type-II states (6.2.20), 


V-II ({A nm }) = (C' 10 ) Alo (C'0i) Aoi (C 11 ) A "(C 2 0) A2O (*i2) 2 “ l/ • (6.2.28) 


Using the definitions (6.2.25) and 

c 20 = \ {a' 2 + a' 2 ) , (6.2.29) 

these may be written as 
i> u (N,M,n,m) = A' N B' M a' n b' m (z 12 ) 2 ~ 1 ' 

(6.2.30) 

for n — m — 0,2 ,... 

It is important to realize that in the particular case of two particles both 
the bosonic limit v —> 0 and the fermionic limit v —> 1 are smooth, so that 
the wavefunctions in (6.2,27) and (6.2.30) are the full solution of the two-anyon 
problem. In particular all the two-fermion wavefunctions, which are 16 


16 We write only the factors depending on the relative coordinate and understand those 
depending on the center of mass coordinate. 
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( 2 l 2 ) m_n+1 ^" n+1 ( 


eB 


1*12 I 


(6.2.31a) 


with m and n non-negative integers such that m — n — 0,2,, and 


(z 12 )”—+ 1 L”- m+1 (^l*12| 2 ) M 


m—n—1 j-m—n—1 
•^n+l 


eB 


1*121 


(6.2.31b) 


for m — n = 0, —2,..., are reproduced by taking v = 1 in (6.2.27) and (6.2.30). 

Furthermore, the whole spectrum of states for two anyons is given by (6.2.17) 
and (6.2.20) with s — 0. In fact for s = 1 we could for example also consider the 
state 


C # 2o(^12) 2+I/ = C*2 0 


However, from the identities 


det 


at at 


det 


t id 


= det 


= det 


2/ J 


a 


t 


Cio 


(* 12 )" 



a 


a 


(6.2.32) 


t\l 


t 

2/ J 


det 


det 


at at 


L \ * * 

it is easy to prove that 17 


2 C 20 


det 


id id 



t id 


1 b 


1 bl 


= det 


t\1 


2/ J 


= det 


2 Coi 

Cm C 02 

2 Coi 

Cio Cn 


= Cf 0 (2 Co2 - C 0 2 !) + (2Cn - C 0 iC 10 ) 


(6.2.33) 


This shows that the state (6.2.32) is a linear combination of the states (6.2.23). In 
general we see that among the independent symmetric polynomials C„ m relevant 
for N = 2, the only one not appearing in (6.2.23) is C 20 * However, from (6.2.17) and 
(6.2.18) it is clear that for any value of s, C 20 always appears in the combination 
(6.2.33) not to produce singularities. Thus, actually no new state is obtained for s > 
0. Identities like (6.2.33) indicate that the states (6.2.17) are not all independent 
of each other. Unfortunately, such identities are not easy to derive for N > 3; in 
(Dunne et al 1992a) the case N = 3 has been considered in detail, but the general 
case is still an open problem. 


17 Eq. (6.2.33) does not contradict the fact that the C nm ’s (m + n < N) are a minimal 
basis in the ring of symmetric polynomials; in general polynomial relations among these 
objects can be written, but it is not possible to solve them in order to express one C nm 
(m + n < N) as a polynomial in the others. 
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6.3 Closed-Form Eigenstates 


In this section we solve directly the eigenstate differential equations (6.1.16) to find 
solutions in a closed-form. These are of course contained within the states discussed 
in the previous section, but for the sake of clarity here we present their explicit 
functional form. This approach also illustrates the separation of the dependence 
on the center of mass coordinate in the N -body anyon problem and yields states 
which are the N -body generalization of those found previously by Y.-S. Wu (Wu 
1984b) in the three-body problem for the related system of anyons in a harmonic 
potential 18 . 

A, 

Our starting point is the “reduced” Hamiltonian H and the corresponding 
“reduced” wave function $ defined in (6.1.17). We shall seek Type-1 19 eigenstates 

A A 

ipi of H and of the angular momentum J c of the form (6.1.24), z.e. 

i>i = m*»r P ( Zl ^ (6.3.1) 

I<J 


so that the eigenstate conditions (6.1.16) on ^ become eigenstate conditions on 

A. 

P. Indeed, for ipi to be an eigenstate of J c , we must require that P be also an 

N 

eigenstate of the same operator J c , since the statistical prefactor n (zijY ™ 

KJ 

A. 

itself an eigenstate of J c , with eigenvalue vN(N — l)/2. Instead, for ipi to be an 
eigenstate of P, we see that P must satisfy a modified differential equation, which 
reads 



—di + cozj 
m 





Bi-Bj 

Zi - Zj 



E — uj 


N 

T 



(6.3.2) 


The general solution of this partial differential equation is not known, but some 
families of solutions can be found quite easily. For example we can make the 
following Ansatz (Dunne et al 1991a) 


P (zu ..., z N ; Zi ,..., z N ) = <S (d+) (zi, ..., z N ) <t> 



(6.3.3) 


where <S( d+ )(zi,... ,z^r) is a totally symmetric homogeneous analytic polynomial 
of degree in zj, and <f> is a function to be determined. Inserting (6.3.3) into 
(6.3.2) we obtain the following differential equation for </) 

18 Note that the harmonic potential problem involves the Hamiltonian (6.1.9) without the 
angular momentum term. But since H and J c commute, the eigenvalue problems are 
essentially the same. 

19 The steps discussed below for Type-I states may be repeated similarly for Type- 
II states, leading to analogous equations whose solutions are listed in (6.3.6b) and 
(6.3.37b). 
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y 4>"{y) +(d + +N+ ^N(N - 1) - yj <j>'{y) 


+ ff-f)^) = 0 (6-3.4) 


where y = (eB/2) ^ \zjf and the primes denote derivatives with respect to y. 

i=i 

This is a differential equation of generalized Laguerre type (see Appendix A), with 
solution 

Hy) = L d n ++N ~ 1+ ^ N(N ~ 1 \y) , (6.3.5) 

where n = ^ - y is a non-negative integer. 

Thus, restoring the exponential factor, our Ansatz (6.3.3) yields the following 
exact iV-anyon Type-I wavefunctions 


Tpi = X\{zijYS^{z u ... ,ZN) 


I<J 


d++N-l+%N(N-l) 

Jjtl 


eB 


£M 


1=1 


exp 


eB 


-tD«I 


1=1 


with energy and angular momentum given by 


E\ = cj In + 


N 


Ji = -N(N - 1) + d + , d 

An analogous Ansatz for Type-II states leads to 

foi = (zi,---,zn) 


= 0 , 1 , 2 , 
= 0 , 1 , 2 , 


* 1 1 


(6.3.6a) 


• 1 


I<J 


x L 


d—-\-^p^N(N—l)+N—l ( eB 


Em 


I-i 


exp 


eB 


-tEW 


J =1 


with 


{ N 2-v \ 

E u =uln+- + d. + -^-N(N-l)) , 

Ju = -^Y^N(N-l)-d. , 


n = 0,1,.. 

d_ = 0,1,.. 


(6.3.6b) 


More general solutions than those in (6.3.6) can be obtained with a more gen¬ 
eral Ansatz than (6.3.3) for the Type-I states, which may be achieved by factoring 
out the dependence on the center of mass coordinate proportional to (z\ +.. .+z^r). 
Such a factorization is of course quite natural because the “complicated” term 
2 in (6.3.2) depends only on the differences of the coordinates. To ac- 


i<J 

complish the required separation we first make the following (linear) change of 
variables 

zi, z \ y ..., zx —► Z,wi,w 2 ,... Wi\r-i (6.3.7) 
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where 


Z = 


(Z\ + Z 2 + . . . + Zjtf) 


y/N 

Wl = -j= (zi ~ z 2 ) , 

W2 = ~Tq ( 21 + 2 2 “ 2 z 3 ) 


(6.3.8) 


w N ~i = 


y/N(N - 1) 


(Z1+Z2 + ... + Zff-1 - (N - 1 )z N ) 


These formulas may be summarized succinctly as 


w 



(6.3.9a) 

N 

= E , 

/= i 

(6.3.9b) 


where are the N vectors (each one with N — 1 components indexed by i ) 
usually considered to represent the roots and weights of SU(N ) (see for example 
(Georgi 1982)). These vectors, together with their main properties, are listed below 


«(« = 


1 1 


V2' V6 ’Vl2 y/(m+l){m + 2y"' ’ y/(N-l)N 


( 2 ) _ _ 


1 1 


u v ' = 


V2’ Vfp Vl2’' ‘’ yV + 1 ){m + 2) ’ ‘ " ’ ^/(JV - 1)JV 


u ( 3 > = 



3’VT2’'” ’ yV + l)(m + 2)’"‘ ’ 


tt 


(m+1 > = 0 , 0 ,...,- 



m 


m 


+ 1’ \/(m + l)(m + 2) ’ ’ ^(JV-l)JV 


= ( 0 , 0 ,..., 0 ,..., 0 ,..., 0 ,- 



IV - 1 


IV 


(6.3.10) 

In terms of these vectors, the positive roots of SU(N) (there are ^N(N — 1) of 
these) are simply 
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r = uW-u^ ( I<J) , (6.3.11) 

and the fundamental weights of SU(N) (there are N — 1 of these) are 

i 

At = E it (/) , i = l,...,N -1 . (6.3.12) 

1=1 


Furthermore, we have 


N 

y =0 for all i = 1,..., iV — 1 ; 

i=i 

(6.3.13) 

N 

= fi* for all ij ; 

i=i 

(6.3.14) 

t^C-0 . — — . 

(6.3.15) 

The transformations (6.3.9) to the new set of variables Z, may be inverted using 
the vectors according to 

*' = VN Z + g“< 

(6.3.16) 

for I = 1,... N. Similarly, 


N-l 

di = r^dz + X) u i ■ 

V N i= 1 

(6.3.17) 


Using properties (6.3.13) and (6.3.14), it is easy to see that 


and 


N -1 


^ djdj = dzdz + ^2 ? 


i=i 


i=i 


N -1 


y: ^ w i®wi * 


J=1 


i=l 


(6.3.18) 

(6.3.19) 


Thus, in the new variables, the differential equation for P becomes 


- dz + w 

m 


\ _ N-* / 2 

+ y ^ Ujw i 





(6.3.20) 
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where the final term in the operator appearing on the left-hand side involves a 
sum over all the positive roots r of SU(N ). The angular momentum operator also 
separates according to 


N 

j = XI ( Zldl - Z*® 1 ) 

1=1 

N-1 

= ( Zd z ~ Zd z ) + ^2 ~ 

i=i 


These properties suggest the following separation of variables, 

P(z,z) = F(Z,Z)G(w,w) . 


(6.3.21) 


(6.3.22) 


Requiring P to be an eigenstate of J implies that F in (6.3.22) must be of the 
form 20 

\Z \. (6.3.23) 

With P separated as in (6.3.22) and (6.3.23), the differential equation (6.3.20) 
requires that the function / satisfy 


/ e B 

F(Z, Z) = Z a f (— 


x f" + (a + 1 — x) f + A/ = 0 , (6.3.24) 

where A is an arbitrary positive integer, x = (eB/2)\Z\ 2 and the primes denote 
derivatives with respect to x . This is once again the generalized Laguerre differ¬ 
ential equation (see Appendix A), giving 


F(Z,Z) = Z a L% 



(6.3.25) 


We remark that a must be an integer to ensure that F be single-valued, and that 
a may take also negative integer values, provided a > —A. in particular when a is 
negative, using the special properties of the generalized Laguerre polynomials, we 
can write 

z -\a\ L ~\a\ (^E\ Z \^ (X ZW4°l |o| (^l 2 ) • 

Therefore the Ansatz (6.3.23) leads to two types of solutions ((ei?/2)|Z| 2 ) 

and zH£W| a | ((ei?/2)|Z| 2 ) which may be succinctly summarized in (6.3.25). 

Having solved for the center of mass dependence, we now address the (more 
complicated!) dependence on the relative coordinates. With F given by (6.3.23) 
the differential equation (6.3.20) becomes the following differential equation for 
G(w,w): 


20 The in the argument of / is inserted for convenience. 
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7V-1 


E( — 

.5 v m / ra i 5 ;'' r ' 


r * d w 


w 


Consider the following Ansatz for G(w,w): 


G= IE-uj I ^ + A G . 

f 

(6.3.26) 



G(w , tn) = wf'wf* ... * g 


d Ar ~ 1 

E w 2 


(6.3.27) 


i=i 


Here the monomial factor depends only on the w% not on the w’s. Inserting this 
Ansatz into (6.3.26) we obtain a differential equation for the function g , 


yg 


”+(N-l + d+^N(N-l)-y) 5 , + (f-^" A ) 5 = ° 


(6.3.28) 


N -1 N -1 

where d = ^ and y = {eB/2) £ |u;;| 2 . In going from (6.3.26) and (6.3.27) 

i=i i=i 

to (6.3.28) we have used the fact that 


n rj = N . 

r> o 


(6.3.29) 


Hence 


|Ei w *i 2 = 2^Ei r - 

z 1=1 iV r>o 


w 


(6.3.30) 


and so 


E^ 


r ■ d w \ ( eB . 

g I — X, K 


N -1 


r • w i \ 2 
r> o / V i=i 


= ^JV(W-D 9' I ^ E to 


(6.3.31) 


i=i 


The differential equation (6.3.28) is once again of the generalized Laguerre type, 
and so we find 


_ d+N- 2 +%N(N-l) f eB_ . 

9 - L e-\n ~a ^ 2 ^ 


(6.3.32) 


Thus, with the Ansatz (6.3.27) we obtain the result that the factor G(w,w) in 
(6.3.22), describing the dependence of P on the relative coordinates, has the form 

/ N -1 

of a monomial in w multiplied by a Laguerre polynomial in ( (eB/2) £ \ w %? 

\ i=i 

However, this function G must also satisfy the appropriate symmetry condition, 
which is the requirement that P, and hence G, must be totally symmetric in the zj 
and zj coordinates. In terms of the new coordinates w , this means that G(w,w) 
must be a Weyl-invariant function of the w and tu. This is because the interchange 
zj zj induces (via (6.3.9b)) the following transformation in the w variables 


w 


a r (w) = w — (r * w)r 


(6.3.33) 
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where r = — u^ J \ Since the roots r are normalized with (r) 2 = 2, this 

transformation is simply the Weyl reflection of w in the hyperplane perpendicular 
to the root r = u^ — u^ J \ To verify (6.3.33) one needs the properties (6.3.15) of 
the u vectors. (Note also that there are \N(N — 1) interchanges zj *-*• zj and 
there are also \N{N — 1) Weyl reflections for SU(N ), one for each positive root.) 

This implies that G(w,w) must be invariant under all ^N(N — 1) Weyl re- 

N -1 \ 

(eB/2) ^2 \w{\ 2 ) is indeed Weyl- 

i= l / 


flections. For the Ansatz (6.3.27), the factor g 


N -1 

invariant because ^ \wi\ 2 is itself Wey 1-invariant. To see this, recall that the 

»=i 

weights in an irreducible representation of SU(N ) may be generated by repeat¬ 
edly reflecting the highest (dominant) weight A in all the positive roots. Therefore, 


l^-H 2 , (6-3.34) 

(where the sum is over all weights /x in the representation {A}) is Weyl-invariant 
since under the Weyl reflection (6.3.33) 


/x * w —> /x * (w — (r * w) r) 

= (/x — (r * /x) r) * w 
= c7r(fi ) * w . 

In particular, for the adjoint representation the weights are just the roots them- 

2 N ~ 1 

selves and so Y2 l r * w \ an d by (6.3.30) also Y2 Weyl-invariant. 

r> o »=i 

The monomial factor 

Af{dj(ty) = wf x wf* ... w (6.3.35) 

in (6.3.27) is not necessarily Weyl-invariant and must be explicitly symmetrized in 
the z variables. From the above discussion, they may be achieved by considering 

M^iw) = M W ( to)) (6.3.36) 

all r 


N-l 

which is a homogeneous polynomial of degree d = ^2 111 w i variables. Note 

»=i 

that in (6.3.36) the summation is over all roots of SU(N); however, vanishes 
if d\ is odd and we need only sum over the positive roots when d\ is even. 

Finally, collecting together the results (6.3.22), (6.3.25), (6.3.27), (6.3.32) and 
(6.3.35), and restoring the exponential factor, we arrive at a new family of Type-I 
eigenstates 

N / pn 

Tk = U (zurZ a Ll(—\Z \ 2 

I<J ' 1 

X *«<„) *£ exp (-f t | 2I | 2 
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with 

f N 

£?! = a; ( — + Ai + A2 
J\ = ^N(N -l) + a + d 

An analogous derivation for Type-II states yields 


(6.3.37a) 


^u = f[(zij) 2 - , 'Z a L a Xi ( 

i<J ' 


eB 


\Z\ 


x M^ d \w) • L 


d+N- 2+ 
A 2 


2 — 1 / 


N(N-i) / eB 


N -1 


£ K 


eB 


JV* 


exp 


-tEw 


1=1 


J =1 


with 


Eh = w I y + Ai + A 2 + a + d + ?^N(N - 1) 
Ju = ~ ■ 2 Z ^ N(N — 1) — a — d . 


(6.3.37b) 


These results (6.3.37) generalize to the iV-anyon case the results of Y.-S. Wu 
(Wu 1984b) for the three-anyon problem in a harmonic potential. Moreover, if the 
parameter a in (6.3.37) is negative, we obtain the “mixed” states introduced in 
(Dunne et aL 1991a,b). 

Finally, we note how the wavefunctions in (6.3.37a,b) are related to those in 
(6.3.6a,b). To make this comparison we use the following property of the general¬ 
ized Laguerre polynomials (see Appendix A) 


Li! +I * +1 ( X + y)= £ L{\(x)-L[l(y) , 

{Ai+A 2 =«} 


(6.3.38) 


and observe that 


N -1 


N 


i*i 2 + £ki 2 = £m 2 • ( 6 - 3 - 39 ) 

i=i i=i 

Then, if one takes the wavefunction in (6.3.37a) and sums over all Ai, A 2 such that 
Ai + A 2 = n we find 


N 


eB 


N 


H(z IJ rZ a M^(w)exp | 


I<J 


1=1 


E ^(f 

{Ai+A 2 =n} 


\zf) 


-l) ( eB 


N -1 


£ K 


1=1 


N 


eB 


N 


(6.3.40) 


= n (*ij) u Z a M w (w) exp ——^2 M 


i<J 


1=1 


x L 


a+d+%N(N-i)+N-i / eB 


N 


n 


£m 


1=1 


5 
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which is precisely the Type-I wavefunction in (6.3.6a) with the totally symmetric 
polynomial separated into center-of-mass and relative coordinates. Thus 

the wavefunctions in (6.3.6) are special linear combinations of the more general 
wavefunctions presented in (6.3.37). 

We conclude this section by observing that all the closed-form wavefunctions 
we have obtained by explicitly solving the differential equation (6.3.2) are special 
cases of the Type-I and Type-II states constructed in Section 2. For example all 
Type-I wavefunctions (6.3.37a) can be obtained by acting with suitable combina- 

A /Q\ 

tions of the step operators Co m , Cio and Cn on the base state ^ ; (Dunne et al. 
1992a). We leave as an exercise for the reader to verify this assertion, but we want 
to point out that the step operator approach of Section 2 is much more general 
than the analytic one. Indeed all Type-I (6.2.17) involving Ci m with m> 1 do not 
appear within the closed-form solution we have presented in this section, but they 
are perfectly good and acceptable eigenstates of our Hamiltonian. Furthermore, 
the algebraic representation of the anyon wavefunctions in terms of step operators 
acting on base states will be extremely useful in computing the canonical parti¬ 
tion functions and in analyzing the statistical mechanical properties of the anyon 
system, as we will see in Chapter 7. This analysis would be much more difficult 
using only the analytic methods related to the solutions of the differential equation 
(6.3.2). 


6*4 Perturbation Theory and the Three-Anyon Problem 


As we have already observed in Section 2, the set of Type-I and Type-II wave- 
functions in (6.2.17) and (6.2.20) is not complete, and the bosonic and fermionic 
limits (u —> 0,1) are not smooth due to missing states. To shed some light on 
the problem of these missing states, we adopt a perturbative approach: start¬ 
ing from an arbitrary bosonic or fermionic state, all of which are exactly known, 
we compute the perturbative corrections to its energy and wavefunction due to 
the non-trivial statistics (Chou 1991b; Me Cabe and Ouvry 1991; Karlhede and 
Westerberg 1991; Khare and Me Cabe 1991; Dasnieres de Veigy and Ouvry 1991; 
Sporre et al. 1992b; Chou et al. 1992). In particular we show that the change in 
energy, AE , to first order in perturbation theory, obeys the following inequalities 
(Karlhede and Westerberg 1991) 


N(N-1) 
V 2 


< AE < 0 


(6.4.1) 


where 


{ v for perturbations around bosons 
v — 1 for perturbations around fermions 


(6.4.2) 


and N is the number of particles. The Type-I states (6.2.17) are characterized by 
energy eigenvalues which are independent of v (see (6.2.19a)), and thus for them 


AE\ — 0 . 


(6.4.3) 
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On the other hand, the Type-II states (6.2.20) have energies linearly dependent 
on v (see (6.2.22a)), so that 


*** = 


(6.4.4) 


As we will see in the following, AE = 0 corresponds to the case when the 
relative angular momentum for all pairs of particles is non-negative, whereas 
AE = — v N(N — l)/2 corresponds to the case when the angular momentum of all 
pairs is negative. Intermediate possibilities when the relative angular momenta for 
different pairs of particles do not have a definite sign, correspond to missing states 
whose wavefunctions are not of the form (6.1.18). The exact analytic expression of 
such wavefunctions is not known at the moment, but some approximate analytic 
solutions to this problem have recently become available (Chin and Hu 1992). 

To prove (6.4.1) and apply perturbation theory, it is useful to write (see (6.3.1)) 


= IIC zijY P(zu...,z N ;z u ...z N ) , (6.4.5) 

I<J 

and to consider the reduced Hamiltonian on the single-valued function P . The lat¬ 
ter satisfies the differential equation (6.3.2) which we rewrite here for convenience 



—di + LJZj 

m 



m 7 ?A Zi ~ Zj 



E — LJ 


N 

2 



(6.4.6) 


To further simplify the formulas, from now on we set m = eB = 1 and consider 
P to be symmetric or antisymmetric depending on whether we perturb around a 
bosonic or a fermionic state. As is clear from its definition (6.4.2), the parameter v 
in (6.4.5) and (6.4.6) is the deviation from the bosonic or the fermionic statistics, 
and is taken to be the perturbation parameter. With these conventions (6.4.6) 
becomes 


(H 0 -i>K)P = 




(6.4.7) 


where 


H 0 


is the free Hamiltonian, and 


N 

X )(~ 2 0/ + zj)5j 

i=i 


(6.4.8a) 


72. = —2 ^ 7 Zij = — 

KJ 



Bi-dj 

Zl - ZJ 


(6.4.8b) 


is the perturbation. For definiteness and simplicity, we will consider only the case 
in which we perturb a system of fermions, so that according to our conventions P 
must be completely antisymmetric. The perturbation theory around the bosonic 
point is more subtle and requires particular care (see for example (Chou et al. 

1992)). 
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For N = 1 the eigenstates of the unperturbed Hamiltonian Ho are the Lan¬ 
dau wavefunctions defined in (6.1.12) (up to the exponential factor that must be 
removed), namely 


A » 




_ — _ z 3 L 3 (~\z 

tt2 i+i(n + j)\ n W 


They satisfy the orthonormality condition 



(6.4.9) 


(6.4.10) 


Note that the exponential factor which has been removed in the reduced wavefunc¬ 
tions, is restored inside the integral measure defining the inner product (6.4.10). 

Since we consider perturbations around the fermionic point, the unperturbed 
AT-body wavefunctions are antisymmetrized products of N Landau wavefunctions 
and can be written as Slater determinants, Le. 


ip n j(zi,...,z N ;zi,...,z N ) ocdet 


i) 


(*i) 


■>P J n\(ZN) 



(6.4.11) 


N N 

where n = ^ nj and j = ^ jj are the total quantum numbers of the N-body 

/= l i=i 

state. 

Since is built out of a finite number of polynomials, it is a polynomial itself 
and furthermore it has a definite angular momentum j . Due to the high degeneracy 
of the problem, it is not immediately clear which fermionic polynomials admit 
a regular anyonic continuation; however, exploiting the conservation of angular 
momentum, we can diagonalize the perturbation 72. in the basis of the definite 
angular momentum functions 

To compute the expectation value of 72. for and thus to find the first order 

corrections to the energy eigenvalues, it is sufficient to consider only one term of 
the sum defining 72., for example 


Tin = -—— (6.4.12) 

Zl - Z 2 

which depends only on the relative motion of particles 1 and 2. In fact, introducing 
Z = (l/\/2)(zi + z 2 ) and the relative coordinate z = Z\ — z^, we can rewrite 72-12 
as follows 

A 

n X2 = 2— . (6.4.13) 

z 

Since is a polynomial in zi,..., zjv &nd their complex conjugates, it is also a 
polynomial in z, Z, Z 3 ,..., zn and their complex conjugates; moreover the relative 
coordinate dependence can be conveniently factorized and we can write 
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*0n,j — ^ ^n° 0 (^5 «o»J*—Jo(^5 ^3? • • • ? , Zjv) (6.4.14) 

no do 

>s # 

where is the normalized Landau wavefunction (6.1.12) depending only on the 

relative coordinate of particles 1 and 2, and 4>n-n 0 j-j 0 Is some (not necessar¬ 
ily normalized) energy and angular momentum eigenfunction depending on the 

remaining degrees of freedom. If is normalized, from (6.4.10) it follows that 

E a a 

(^n—n 0 ,j—jol^n—n 0 ,j —Jo) = ^ • (6.4.15) 

no do 


Using the decomposition (6.4.14) and observing that the ft s are orthonormal be¬ 
cause they have different energies and angular momenta, we immediately find 


{ 1 &nj\'R'12\ 1 tpnj) — ^ ^Ptiq 1^12 \^Pno ) * (^n—n 0 ,j — Jo |^n—n 0 ,j — Jo ) 


(6.4.16) 


n 0 ,Jo 


The matrix element 1^-12 IV^q) can be easily computed using the properties of 
the generalized Laguerre polynomials reported in the Appendix. We leave the ex¬ 
plicit calculation as an exercise to the reader and we report here only the extremely 
simple result 

<^l^i2|^> = {? * j P*° Q (6.4.17) 

Inserting this into (6.4.16), we get 

= 2 y ^ y ^ i$n—rio,j—jo\4 ) Ti—Ti 0 ,j—j 0 ) * (6.4.18) 

no Jo<0 


Due to the symmetries of the problem (the full perturbation 72. is symmetric un- 

A » 

der particle exchanges, while is antisymmetric) the total first order energy 
correction AE is 


AE = {tp n j\(-i> n)\ip n j) = -PN(N - l)(Tp n j\K 12 \ii> n ,j) 


. N(N - 1) 

= —V --—:- - 



A A 

(^n—n 0 ,j—Jol^n—n 0 ,j—Jo) * 


(6.4.19) 


no Jo<0 


Each term in the sum appearing in the last line of (6.4.19) is positive definite, thus 


AE <0 


(6.4.20a) 


On the other hand the maximum value of this sum is 1 (see (6.4.15)), so that 


AE> ^ T zll 


This concludes the proof of (6.4.1). 

The bound (6.4.20a) is saturated (z.e. AE = 0) only if 


(6.4.20b) 
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A A 

ypn—n 0j j—j 0 \4>n—« 0 J—jo) = 0 (6.4.21) 

for all jo < 0. If this is true, the wavefunction is decomposed with respect to 
all pairs of particles in relative wavefunctions with non-negative angular momenta. 
As anticipated at the beginning of this section, this occurs precisely for the Type-I 
wavefunctions (6.2.17) constructed with the step operator approach. The other 
bound (6.4.20b) is saturated (z.e. AE = —i>N(N — l)/2) when 

{^n—noj—jol^n—noj—jo) = 1 (6.4.22) 

for all jo < 0. This means that ipnj contains relative wavefunctions for all pairs 
of particles whose angular momentum is always negative. As previously remarked, 
this is precisely the case of the regular Type-II wavefunctions (6.2.20). 

This analysis gives us a nice physical interpretation of the two families of 
anyonic wavefunctions which we have constructed using the step operators C„ m , 
and also furnishes some intuition on the nature of the missing states. The latter 
can indeed be characterized by saying that the relative angular momentum does 
not have a well defined sign for all pairs of particles. Also form this point of view, 
it is clear why the two-body problem (in which there is obviously only one pair) is 
completely solved by the Type-I and Type-II states of section 2 without missing 
states. 

We now specialize our discussion to the three-body case and compute in par¬ 
ticular the first order perturbative corrections to a fermionic wavefunction whose 
anyonic continuation is missing from our Type-I and Type-II states. To this aim 
it is convenient to use the Jacobi coordinates Z^w defined in (6.3.8), so that 
the differential equation (6.4.6) becomes 

( Hem. + &rel) P = (,E ~ |) P (6.4.23) 


where 


^c.m. 


tfrel 


(-2 d z + Z) d z , 

2 

^2 (- 2 ^ + m) d Wi 

1=1 



+ V3d W2 | d Wl - V3d W2 

Wi + V3w 2 W\ — y/ZW2 


(6.4.24) 


The center of mass dynamics governed by H Ctm . is trivial and decouples from the 
relative motion. From now on we will neglect the center of mass contributions 
both to the wavefunctions and to the energy spectrum, and will focus only on the 
problem of the relative motion. As we did before, we treat the ^-dependent terms 
of H re i in (6.4.24) as a perturbation, and write 

A A ^ 

H tei = H” el -2 vK(wi,w 2 ;wi,W2) . 


(6.4.25) 
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We are interested in particular in the perturbative corrections to the non degen¬ 
erate fermionic state described by 



-{w\w 2 - Wi w 2 ) 


(6.4.26) 


Despite its appearance this wavefunction is completely antisymmetric under all 
particle exchanges. In terms of the Jacobi coordinates the exchange of particles 1 
and 2 corresponds to 

W\ -► — W\ , 

w 2 —► w 2 ; 

the exchange of particles 1 and 3 corresponds to 


(6.4.27a) 



W\ — y/3w 2 



— y/3wi — w 2 
2 


(6.4.27b) 


and finally the exchange of particles 2 and 3 corresponds to 



Wi + v ^ 3^ 2 



y/3wi — W 2 

2 


(6.4.27c) 


It is immediate to verify that (6.4.26) is indeed antisymmetric under all transfor¬ 
mations (6.4.27). Up to a numerical factor, the wavefunction (6.4.26) is the same 
as that in (6.1.26) written in Jacobi coordinates, and in the closely related problem 
of three anyons in a harmonic oscillator potential it represents the non-degenerate 
fermionic ground state (Wu 1984b). 

Writing 

P = P°+uP 1 +v 2 P 2 + ... , 

n ,99 (6.4.28) 

E = E° + vE 1 + v 2 E 2 + ... , V ' 


and equating the coefficients of the different powers of v in the eigenstate equation 
for the relative motion 21 


[H? el -2un(wi,w 2 ;wi,w 2 )\ P = (E -l)P , 

we get 

(6.4.29a) 
(6.4.29b) 
(6.4.29c) 


(iJ r 0 el - E° + 1) P° = 0 , 

(iJ r ° el - E° + 1) P 1 - (2n( Wl ,w 2 ;w u w 2 ) + E 1 )P° =0 , 

(iJ r ° el - E° + 1) P 2 - (2K(w u w 2 ; w u w 2 ) + E 1 ) P 1 - E 2 P° = 0 


21 If also the the center of mass is considered, (E — 1) must be replaced by (E — 3/2) in 
the right hand side. 
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Using the explicit expressions of H® el and P°, we can easily solve (6.4.29a) and 
find 



The second equation allows to determine E 1 and also P 1 , but only up to terms 
proportional to P° because of (6.4.29a). We leave as an exercise to the reader to 
verify that 


P 1 


E 1 


^^wiW2 ln(wit/;i) 

— (V3w 2 — wi)(V3wi + W 2 ) ln[(\/3wi + W 2 ){'f$w\ + W 2 )} 

— (\/3u>2 + wi)(y/3wi + W2)ln[(\/3w2 — wi)(V3w 2 — tZ>i)] j 


3 

2 


(6.4.30) 


solve (6.4.29b). Thus, to first order we have (Karlhede and Westerberg 1991) 

E = 2-\v , 

2 (6.4.31) 

P = P° + i>P 1 . 


Unfortunately this solution cannot be extended in analytic form beyond the first 
order. We note that the energy eigenvalue in (6.4.31) has a negative ^-derivative; 
its value (—3/2) corresponds to that of the Type-II wavefunctions of Section 2 
(see (6.2.20) for N = 3). One is therefore tempted to identify (6.4.31) with the 
perturbative expansion of some Type-II states. However this cannot be the case 
because, as we have already remarked, the state described by P° is not contained 
in the fermionic limit of (6.2.20). This means that the anyonic continuation of P° 
has higher perturbative corrections and its energy depends non-linearly on v. This 
expectation is indeed confirmed by recent numerical investigations (Sporre et al 
1991a; Murphy et al 1991; Sporre et al 1991b). 


6.A Appendix 

In this appendix we list some of the main properties of the generalized Laguerre 
polynomials which are useful to derive the results presented in the previous sec¬ 
tions. The generalized Laguerre polynomial, U n (x)^ satisfies the differential equa¬ 
tion 

x f n + (j + 1 — x) f + n / = 0 . (6.A.1) 

The Rodriguez formula for L? n is 

Li(x) = (e -x £ n+J ) , (6.A.2) 

and L £ is a polynomial for j > — n, n 6 Z + . Some special cases are 
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4 

L{ 


= 1 , 

= j + 1 - x , 


= L n , the ordinary Laguerre polynomial , 


r—n _ 


(-i) 

n! 


n 


X 


n 


They form an orthogonal set for each j > — 1, 



00 ... (n+iY 

= -— -~r^-S nvn 


o 


n! 


An explicit expansion for L 3 n is 




(6.A.3) 


(6.A.4) 


The U n satisfy a two-term recurrence relation 

nL 3 n {x) = (2 n + j - 1 - x) L^^x) - (n + j - 1) L 3 n _ 2 (x) (6.A.5) 

for n = 2,3,4,.... When the upper index is a negative integer, we have 

L-Hx) = (-xy i -^^-Li_ j (x) . (6.A.6) 

71 . 

The polynomials L 3 n satisfy the following differential-difference relations 

±Li = -Li+_\ , (6.A.7a) 

x-^-Li = nLi-(n + j)Li_ 1 , (6.A.7b) 

Finally, the generalized Laguerre polynomials satisfy the following addition for¬ 
mula, 

L » 1+ ... +ak+ k-i {xi+ +Xk)= £ L? l '(x l )...L%(x k ) . (6.A.8) 

{»! + ...+ifc=n} 



7. Statistical Mechanics of Anyons 


In this chapter we are going to study some statistical mechanical properties of 
systems of many anyons, and in particular we will discuss the partition functions, 
the virial coefficients, the magnetization and the magnetic moment. However only 
very partial results can be obtained in this context, because the exact solution 
of a gas of anyons is not known. In fact, in contrast to the bosonic or fermionic 
case where the statistics is implemented by hand on the many body Hilbert space 
by constructing completely symmetric or antisymmetric products of single particle 
wavefunctions, for anyons the complicated boundary conditions for the interchange 
of any two particles require the knowledge of the complete many-body configura¬ 
tions. As we have pointed out in the previous chapter, only the two-body problem 
is exactly soluble for anyons, and hence only the two-body partition function can 
be computed exactly. Since the thermodynamic limit cannot be performed, one 
has to resort to approximate or alternative methods to study the statistical me¬ 
chanics of anyons. For example if the thermodynamic functions are analytic in 
the particle density, it is well-known that the low density, or equivalently the high 
temperature limit, of a (free) gas can be investigated using the virial expansion 
(see for instance (Huang 1987; Reichl 1980)) in which the equation of state that 
relates the pressure P, the temperature T and the density /?, can be expressed as 
follows 

P — p &B T (l + &2 P "I" a 3 P* "I" * * *) (7*1) 

where k-Q is the Boltzmann constant and a 2 ,a 3 ,... are the so-called virial coef¬ 
ficients. For bosons and fermions, (7.1) is known to be consistent and one can 
assume that a similar expansion holds also in the case of particles with fractional 
statistics. 

When a 2 ,fl 3 ,... are all vanishing, the equation of state is that of a classical 
ideal gas (P = pk-Q T). The second virial coefficient a*i depends only on the two- 
body quantum correlations, and it is positive for bosons and negative for fermions, 
signaling the presence of a statistical attraction or repulsion respectively. Since the 
two-body problem is completely soluble for anyons, the second virial coefficient can 
be computed exactly also for intermediate statistics (Arovas et al . 1985), and it 
turns out that it continuously interpolates between the negative bosonic value and 
the positive fermionic one. The third virial coefficient a$ depends on the quantum 
two- and three-body correlations. It is known for bosons and fermions, but not for 
anyons since the three-anyon problem is not solved. However, even in the case of 
intermediate statistics, it is possible to derive some qualitative features of a$ using 
semiclassical arguments (Bhaduri et al. 1991) or, as we will see, using the Type-I 
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and Type-II states we constructed in the previous chapter (Dunne et aL 1992b). 
For the higher virial coefficients instead, only very preliminary and perturbative 
results are available at the moment (Comtet et aL 1991). 

Remarkably enough, the exact solution of the two body problem allows to 
extract a few interesting features of the anyon system, like the presence or absence 
of spontaneous magnetization, or the presence or absence of parity violating effects, 
or the relation between the canonical magnetic moment and the fractional orbital 
angular momentum that we introduced in Chapter 5. However it would be even 
more remarkable and interesting to extend these results to the thermodynamic 
limit and deduce macroscopic effects that, at least in principle, could be detected. 
Lacking for the time being such possibility, we will content ourselves to compute 
the canonical partition functions and derive some consequences using standard 
techniques of statistical mechanics. 

7.1 Partition Functions 

The canonical AT-body partition function, Z^, for a system of free particles is 
clearly divergent because each energy level is infinitely degenerate. Thus, in or¬ 
der to get a finite Z one has to break the degeneracy of the energy spectrum. 
This can be done, for instance, by confining the particles in a finite volume or, 
as we will do, by introducing a confining interaction, like a harmonic potential 
with frequency a;, which serves as a regulator. Since we will be interested in the 
magnetic properties of the anyon gas, we introduce also an external magnetic field 
B as we did in the previous chapter. The coupling to a magnetic field alone is not 
enough to regularize the theory and remove the degeneracy of the energy levels. 
Indeed, even though the presence of the magnetic field organizes the spectrum in 
discrete Landau levels spaced by the cyclotron frequency lj c = eB/m , each level 
remains infinitely degenerate in an infinite volume, thus yielding divergent parti¬ 
tion functions. Therefore hereinafter we will consider a system of N particles in 
the infinite plane interacting both with a harmonic potential and with a magnetic 
field; sometimes the latter will be switched off to simplify the situation and render 
the formulas more manageable. 

For bosons and fermions it is well-known how to compute the N- body par¬ 
tition functions Z#. First one solves the single-particle problem and then forms 
all iV-body wavefunctions by constructing all Slater permanents (for bosons) or 
determinants (for fermions) with the single-particle wavefunctions. In this way 
the energy levels and their degeneracies follow immediately, and one can easily 
compute Zn. For example, for particles of mass m interacting in the plane only 
with a harmonic potential of frequency a;, the single-particle wavefunctions are (cf 
( 6 . 1 . 12 )) 

~ z j Li(mu\z\ 2 )e-h m ^ (7.1.1) 

where n = 0,1,2,..., j = — n, — n + 1,... and L 3 n is the generalized polynomial of 
degree n (see Appendix 6.A). The energy corresponding to (7.1.1) is 

E 3 n =u (2n + j + l) , 


(7.1.2) 
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while the canonical one-body partition function is 


OO 

Z\ = y t e _/3£ " 

n=0 

j=-n 


(7.1.3) 


where, as usual, j3 = 1/(&b T). If we define j* = j + n so that j' = 0,1,... and 
= u(n + j* + 1), we can easily evaluate the summations in (7.1.3) and get 


—P E l! — e -P w ( n +j , + l ) 

n=0 n=0 

j'=0 j '—0 



e~P w _ 1 

(1 - e~P “) 2 ~ 4 s inh 2 


(7.1.4) 


Obviously Z i is the same for bosons, for fermions and for anyons because no 
statistical effect is present in the one-body problem. 

Let us now compute the two-body partition function Z 2 . For simplicity we 
first consider the case of bosons. The two-body bosonic wavefunctions are totally 
symmetric products of the single-particle wavefunctions (7.1.1). If (ni,j{) and 
(^ 2 ? J 2 ) are respectively the labels of the two component factors, the total energy 
is 

E 3 n\H = U (m + j[ +n 2 + j' 2 + 2) , (7.1.5) 

and the partition function is 


Z bos _ 

2 = 



-P K 


3 1 ,3 2 
1 > n 2 




(7.1.6) 


where the symbol ^ denotes the sum over all symmetric states meaning 

J 2 }s 

that the state characterized by the ordered set {ni, j[j7i2jj 2 } must be identified 
with the one characterized by { 712 , j 2 ,^ 1 To calculate Z 2 it is convenient to 
make the following observation: If we sum independently over all labels m, 712 
and J 2 , we clearly overcount the bosonic states because we do not implement the 
identifications due to symmetry; however if we add to this sum the contribution 
of the diagonal states for which n\ = 712 and j[ = j 2 , we get exactly twice the 
wanted bosonic states, that is 



+ E = 


E 


«1 

3*1 


n 2 


«1— n 2 
■t _ ■! 


{«iJ 5;«2 ii 2 }s 


(7.1.7) 


The bosonic partition function, Z 2 oa , then follows immediately from (7.1.7) and 
(7.1.5), and we get 
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Z bos 
2 


1 

2 


^ ^ ^ ^ e — /? a, ( n i+Ji+« 2 +J 2 + 2 ) _|_ ^ ^ (2ni+2j[+2) 


1 

2 


ni = 0 ri2=0 

ii=o j2=o 

0 -2/3o> 


ni=0 

J>0 


—2/3 w 


+ 


(_(1 - e -^) 4 (1 - e " 2 ^) 2 J 

cosh(/? u>) 


(7.1.8) 


8sinh 2 sinh 2 (/?u>) 

The fermionic two-body partition function 


^ er 


E 


-/3 e„\: 


i[j’ 2 

n< 


(7.1.9) 


{ni,ib n 2 »i£)a 


where the symbol denotes the sum over all antisymmetric states, can 

{ni,j';n 2 ,j2}a. 

be computed in an analogous way. A moment thought reveals that 



- E = 


E 


(7.1.10) 


n i n .2 

•/ •/ 
^2 


ni=n 2 

•/ _ •/ 

•?1“ -?2 


{ni,ij;n 2 ,i2}a 


The diagonal states for which two quantum numbers are the same, are clearly 
impossible according to the Pauli exclusion principle, and thus must be removed. 
fYom (7.1.10), (7.1.9) and (7.1.5), it is easy to get 


Z fer 
2 


1 

2 


E E /3aj(n 1 +jJ+n 2 +j2+2) 

ni=0 n 2 =0 
J[=0 J2=0 



/3 a; (2ni+2jl+2) 


1 

2 


e -2(3<*> 

(1 - e - ^ w ) 4 

1 


e -2(3<*> 

(1 — e -2 ^^) 2 


8sinh 2 sinh 2 (/?u>) 


(7.1.11) 


Let us now turn to the three-body partition functions for bosons and fermions. 
The bosonic three-body wavefunctions are totally symmetric products of three 
single-particle wavefunctions (7.1.1), and hence can be labeled by (^2,^2); 

(^3> .7*3 )■ By taking into account the symmetry imposed by Bose statistics, it is not 
difficult to prove that, in obvious notations, 


+ 3 EE+ 2 E = 6 

Tli 712 71 3 Tli —712 713 7li=7l2=7l3 

•* •/ •/ ■/ _ •/ •/ ■/ _ ■> _ •/ 



E 

{nij[;n2j^}s 


(7.1.12) 


Therefore, since 
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/ ■/ •/ 


Eh\'?n 2 ,n 3 — u> (n x + j[ +ri2 + j*2 + nz + + 3) , 


(7.1.13) 


the partition function Z^ os is given by 


ghos _ ^ 


l 

6 


OO OO OO 

E E E 

Til =0 712 = 0 713=0 

j[= 0 J*2 = 0 J*3=0 


—/3 (t> (7ii + j[ + 7i2 +j 2 + n 3 +ia+3) 


oo oo 


+ 3 E E 


—/3 a; ( 2 ni+ 2 jJ+Ti 3 +j 3 + 3 ) 


oo 


+*E 


—/3 a; (37i! +3 j[ +3) 


7l!=0 713=0 

j[=0 j ' s =0 


71 i =0 

j[=0 


1 

6 


e -3(3u> 


(1 _ e~ 0 w ) 


+ 3 


e -3(3u> 


(1 - e -pu>y (1 - e -2po>) 


2 


+ 2 


e ~3/3 u> 


(1 — e -3 ^ w ) 


cosh(3/?u>) + 2 cosh 2 
32 sinh 2 sinh 2 (£u;) sinh 2 

(7.1.14) 

In the case of fermions we must exclude all states in which at least two quantum 
numbrs of the constituent factors are the same; simple combinatorial considerations 
lead to the following result 



-3 E E+ 2 E - 


6 


E 


(7.1.15) 


Til 712 71 3 

•/ •/ •/ 

3l 3 2 0z 


711 = 712 713 

•/ _ •/ ■/ 

0z 


711=712=713 

•/ _ •/ _ •/ 

3 i — 32 — 3 z 


{niJ[;n 2 1 a 


Using (7.1.13), we get 


^3 er 


1 

6 


OO OO OO 

E E E 


—/3 a; (ni+jfJ+7i2+i2+ n 3+i3+^) 


Til 

ii 


:0 ti 2 =0 713=0 

0 j'=o i'=o 


cx? cx? 


00 


-3 V 5>- (27ii+2jJ+7i3+j3+3) _|_ 2 ^ ^ e —/3 o;(3tii+3jI+3) 


71 i =0 713 = 0 

j l=o is=o 


1 

6 


—3/3a; 


—3/3a; 


(1 — e~P w ) 


-3 


(1 - e -pu>y (1 -e-W") 


71 i =0 

il=° 


+ 2 


e —3pu> 


(1 - e~ 3 P w ) 


1 +2 cosh 2 


32 sinh 2 sinh 2 (/?u>) sinh 2 (^f^) 

(7.1.16) 

This method can be easily applied to compute the bosonic and fermionic parti¬ 
tion functions for any number N of particles, but the resulting formulas are more 
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and more cumbersome when N increases. However, such a procedure cannot be 
applied to anyons of intermediate statistics. In fact, in this case, the N -body wave- 
functions cannot be constructed from products of single-particle wave functions, 
because one has to represent the braid group and not simply the permutation 
group. An alternative approach is therefore needed. This is provided by the step 
operators introduced in Chapter 6 (Dunne et al. 1992b). The idea is to work in the 
anyon gauge in which case the Hamiltonian is identical for bosonic, fermionic or 
anyonic systems, but the many-body wavefunctions have different symmetry prop¬ 
erties for different statistics. The Hamiltonian for N particles interacting with an 
external uniform magnetic field B and a harmonic force of frequency a;, is (in 
complex coordinates) 


N 


2 _ * e 2 B 2 , , 2 , 1 2 


H = J2[--d I d I + . 

\ m 8m 

7=1 x 


Zl | + -mu) \Z! 


eB 

2m 


N 


( zid! - z T di ) , 


7=1 


(7.1.17) 

which is a simple generalization of (6.1.9). We are interested in the eigenvalue 
problem H 'ip = E 'ip for multi-valued functions For convenience let us define 


Q* = 4ur + 


B 2 


2 _ a . .2 i ^ ^ a . .2 i . .2 


m 


= 4u) + u) 


c ’ 


and 


i> = exp l^l 2 ^ ^ ’ 


A 

so that the “reduced” Hamiltonian acting on is 

N£2 / 2-1 1 

H = —— + ^2 f ~—didi + -w_z/<9/ + -a;+z/d/ 

7=1 ' 


(7.1.18) 


(7.1.19) 


(7.1.20) 


where u>± = i? ± u) c . In analogy with (6.2.3a) the Hamiltonian H can be written 
as follows 



Nfi 

2 



N 

^ a\ai + u)- 

7=1 



(7.1.21) 


where the operators 


a 


t 

7 





aj = di 


bi = dj 


(7.1.22) 


satisfy the Heisenberg algebra |a/,ajj = |^ 7 ,&jj = <S/j, with all other commuta¬ 
tors being zero. This problem can be analyzed in complete analogy with the one 
discussed in Section 6.2. The symmetric step operators 



97 


N 


c nm = J>} n 4 m 


(n + m < N ) 


(7.1.23) 


7=1 


act on the on the base states ^ = Yl { z ijYi and V'n = FI ( 2 /j) * > which 

/<J /<J 

entirely encode the anyonic statistics (here we take 0 < u < 2), and the following 
two families of regular anyonic wavefunctions can be constructed 


1 N-n 


*=lin (Cnmf™ $ 


7.(0) 


(7.1.24a) 


n=0 m=0 


for Type-I, and 


1 N—m 


( 0 ) 


■4>u =n n ^ v i 


(7.1.24b) 


m=0 n=0 


for Type-II. Here A nm are non-negative integers. For N = 2 the states in (7.1.24), 
completely exhaust the space of wavefunctions for anyons of arbitrary statistics. 
The energy corresponding to (7.1.24a) and (7.1.24b) is respectively 


El ({A nm }) =y!2 + ijV(JV - l)u-v 


1 N—m 


+ 9 ^2 ^2 Xnm ( nu ++ muj -) > 


(7.1.25a) 


m=0 n=0 


and 


Eu ({A„ m }) =jO + \n(N - l)w+(2 - v) 


1 N—m 


^ A MV 

+ - ^2 S Xnm ( nw++ mu> - ) 


(7.1.25b) 


m=0 n=0 


With the anyonic states represented as in (7.1.24), it is immediate to compute the 
corresponding partition functions. In fact, the contribution of the Type-I states 
(7.1.25a) to the partition function is 


4>= E « 

{m } 


-0£l({A„ m }) 


1 N—n oo 


(7.1.26) 




nm 


n=0 m=0 A nm =0 


The summations over A nm can be easily done and, with trivial manipulations, one 
gets 


7 ( 1 )_ 

z jv - 


e -f(3(2- 


N(N-l)w-is 


J_TT _ l _ 

22N e=i sinh sinh [f (w_ + - l)o>+)j 

(7.1.27a) 
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The corresponding contribution zffi of the Type-II states (7.1.24b) to the parti¬ 
tion function can be computed in the same way or simply obtained from (7.1.27a) 
with the substitutions u>+ <-+ u>_, v <-+ 2 — u, yielding 


* i N 

v(H) _ N(N-l)«>+(2-v) 1 TT 

Z N - e 2 4 o2N Xi . 


2<1N t= i sinh sinh |^§ (u>+ + (^ - l)a;_)J 

(7.1.27b) 


Hence the total partition function for our boson-based anyons is given by 


Z n (u,B) = 


(ii) 


2 in 


N(N—1) u>_ (1—i/) 


TV 

n 


t=i sinh ^§a;_^ sinh |^f (a+ (^ - l)u>+) 


TV 


+ c -f JV(JV-1)W + (1-I/) JJ _ 


fci sinh (f u>+^ sinh [§ (a/+ 4- (* - l)u>_)] 

(7.1.28) 

Note that this partition function enjoys the property Zn(u, B) = Z ^(2 — z/, — B) 
for any JV and Z2(0,i?) = Z^^B) for the special case of two anyons. Moreover, 
we observe that Z\ is obviously independent of z/, and that Z 2 in (7.1.28) is the 
complete partition function for two anyons of arbitrary statistics, and smoothly 
tends to the bosonic and fermionic ones for v —> 0,2 and v —> 1, respectively. For 
N > 3, Zn(u, B) does not reduce to the known bosonic and fermionic partition 
functions for v = 0 and v = 1, since the states (7.1.24) do not cover the entire 
spectrum. 

Let us exemplify this for the special case of B = 0. In this case u>+ = u;_ = 2u> 
and (7.1.28) gives 


Z 2 (z/,0) 

^(M) 


cosh [(1 — z/)/?u>] 

8 sinh 2 sinh 2 (/?u;) 

cosh [3(1 — v)f}u\ 

32 sinh 2 sinh 2 (/?u>) sinh 2 


(7.1.29) 


By comparing these with (7.1.8) and (7.1.11), we see that Z 2 (0,0) = Z\° % and 
Z 2 (1,0) = Z| er . These qualities do not hold for Z${v,Q). However, let us note 
that the difference with respect to the exact results (7.1.14) and (7.1.16) is the 
same for v = 0,1. Indeed, 


Z 3 bos - Z 3 (0,0) = 


cosh 2 (^) 

16 sinh 2 sinh 2 (/?u>) sinh 2 



Z 3 (1,0) . 


(7.1.30) 
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Furthermore, we note that for N = 3, the coefficient of the contribution at first 
order in z/, i.e. 0)| 0 , coincides exactly with the one computed with nu¬ 

merical methods in perturbation theory around bosonic statistics (McCabe and 
Ouvry 1991). 

Having established that (7.1.28) is the complete partition function for two 
anyons in a magnetic field B and a confining harmonic potential of frequency a;, 
in the next sections we will compute the virial coefficients and the magnetization 
properties. 


7.2 Virial Coefficients 

The low density or equivalently the high temperature limit of a gas of particles can 
be described using the virial expansion, in which the equation of state is expressed 
as a series in the particle density p (see (7.1)). The coefficients in this expansion 
are known as virial coefficients. To set up the notations, let us first recall the 
basic definition of the first few of these. In the low density approximation, the 
gran canonical partition function Q can be written as a cluster expansion (see for 
instance (Huang 1987; Reichl 1980)) 


Q = exp 


oo 

A be z l 

i =1 


(7.2.1) 


where be are the so-called cluster integrals (Ursell 1927; Pais and Uhlenbeck 1959), 
^ is the fugacity and A is the area of the system (we obviously limit our consid¬ 
erations to the case of two dimensions; in higher dimensions A would be replaced 
by the volume V). On the other hand, using standard techniques of statistical 
mechanics, the pressure P and the density p are given by 


P = k B T 



7 lnQ 

A 

4 

\lnQ 

A 


A,T 


Upon inserting the cluster expansion (7.2.1) into (7.2.2), we get 


(7.2.2) 


oo 

P=k B T^b t z l 

e=i 

oo 

p = ^21 be z e , 
i=i 


(7.2.3) 


so that the virial expansion of the equation of state (7.1) becomes 


oo 



e=i 


oo 

J2*bez e 

i=\ 


1 + a 2 


OO 


y^/b e z e 


+ 


oo 

Y/bez* 


i=i 



(7.2.4) 
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If we now expand both sides of (7.2.4) and equate the coefficients of equal powers 
of we easily obtain the expressions for the virial coefficients in terms of the 
cluster integrals; the first two are 

b 2 

a 2- J)2 

<z 3 = 4al - 2^| . 

We can make further progress if we use the definition of the gran canonical partition 
function 

oo 

Q=J2 zN Z n (7.2.6) 

N =0 

to express the cluster integrals bg in terms of the TV-body canonical partition 
functions Zff. From the comparison of (7.2.1) with (7.2.6), we obtain 


(7.2.5a) 

(7.2.5b) 



b 2 

b 3 


Z i 
A ’ 

2 Z 2 - Z\ 


2 A 


1 


3 Z$ — 3 Z 2 Z\ -f Zi 



(7.2.7a) 

(7.2.7b) 

(7.2.7c) 


Actually all these expressions are meaningful only in the infinite area limit, so that 
particular care must be used in their evaluation. 

In two spatial dimensions, one has 



(7.2.8) 


where At is the thermal wavelength of a particle of mass m and is given by At = 
(27r/J/m) 1 / 2 . Therefore from (7.2.7a) we have 


1 



Zi 

A 


(7.2.9) 


Using the expression of Z\ given in (7.1.4) and taking the high temperature limit 
(i.e. P —> 0), we are lead to the identification 


U) 


2 




2tt 1 
Pm A 


(7.2.10) 


Thus the harmonic frequency u> can be thought as being related to the inverse of 
the area A, and hence the infinite area limit (A —► 0) becomes equivalent to the 
zero harmonic frequency limit (u> —> 0). This is of course to be expected because 
putting the system in a harmonic potential confines the particles to move in a 
finite region; removing the harmonic force ( Le . letting u> —> 0) amounts to remove 
any restriction on the space (i.e. A —> 0). 



101 


Let us now evaluate the cluster integrals and the first viriaJ coefficients in the 
u> —> 0 limit. For simplicity we will use the anyon partition functions Zn in (7.1.28) 
with B = 0; the case with a non-vanishing magnetic field is left as an exercise for 
the reader. 

Recalling that the harmonic regularization requires to multiply be by a normal¬ 
ization factor l d ! 2 where d is the number of space dimensions 22 , the second virial 
coefficient turns out to be 



(7.2.11) 


in agreement with the original derivation (Arovas et al. 1985) which is based on 
different methods (cf also (Comtet et al. 1989; Blum et al. 1990; Sen 1991a)). We 
remark that both Z\ and 2Z 2 (z/, 0)/Zi are divergent in the limit u> —> 0 - in fact 
they are 0(l/u> 2 ) - but the divergence cancels in the combination 



z 2 M) 


appearing in (7.2.11). The second virial coefficient a 2 (z/) interpolates continuously 
between the bosonic value 


a 2 os = <z 2 (0) = - i , (7.2.12) 

and the fermionic one 

<4 er = a 2 (l) = iA| . (7.2.13) 

This is an obvious consequence of the fact that the two-body partition functions 
is exact for arbitrary statistics v . 

Let us now turn to the third virial coefficient which is given by 


° 3 ( v ) = 4<z 2 (z/) 2 — 2^1 

= 4a 2 (z/) 2 — 2 Ai lim 

a >—*-0 



5»M) 


-3Z 2 (i/,0) + Z 2 


(7.2.14) 


22 For two particles the harmonic potential is 

W 2 (z? + z\) = 2a; 2 Z 2 + ^ - z\) 

where Z = (l/\/2) (z i 4- 2 2 ) is the center of mass coordinate. Thus, in the two-body 
problem u>c. m . = 2 w 2 , and in general in the ^-body problem u>c. m . = £ w 2 . Since it is 
cJc.rn. which is actually related to the area according to (7.2.10), the two-dimensional 
cluster integral bg must be multiplied by L 
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If we substitute in (7.2.14) the explicit expressions of the partition functions given 
in (7.1.29) and compute the u> —> 0 limit, we find a divergent result. However, 
remarkably enough, the coefficients of the divergences turn out to be independent 
ofv , and thus the difference a 3 (u) — a 3 {y 9 ) is well-defined and finite. Furthermore, 
a simple calculation shows that 

— a 3 (u 9 ) = 0 . (7.2.15) 


This result is the consequence of highly non-trivial cancellations and was first 
pointed out in (Dunne et al. 1992b). 

We remark that the third virial coefficient is finite for bosons and fermions; 
indeed one has 23 

^° 8 = <4 er = ^4 . (7.2.16) 

This result can be easily obtained if we use the complete bosonic or fermionic par¬ 
tition functions given in (7.1.14) and (7.1.16) respectively. As previously noted, 
these are not the smooth limit of Z 3 (v, 0 ) for v —> 0 , 1 , and in fact the differ¬ 
ence given in (7.1.30) precisely cancels the divergences in the cluster integrals and 
renders the third virial coefficient finite. 

The presence of divergences in a 3 for fractional statistics could signal either 
the failure of the virial expansion due to long-range statistical interactions, or the 
presence of “missing states” which are not included in (7.1.24). This last possibility 
seems to be suggested by recent numerical (Sporre et al. 1991a; Murphy et al 
1991) and perturbative (McCabe and Ouvry 1991; Comtet et al. 1991; Khare and 
McCabe 1991) calculations. To shed some light on this problem, it is therefore 
necessary to investigate in more detail the third virial coefficient a 3 (v). To this 
purpose, let us turn back to the three-anyon problem, described by the reduced 
Hamiltonian (see (7.1.20) for B = 0) 



3w + ^ 


1=1 


— djdi + u) zj di + u> zj dj 
m 


(7.2.17) 


As we have seen in Chapter 6 , H can be nicely split into a center of mass part and 
a relative part by introducing the Jacobi coordinates (cf (6.3.8)) 






(^1 + +23) 


{z\-z 2 ) , 


(Zi +z 2 - 2z 3 ) 


(7.2.18) 


Indeed we have 


H — Hem. 4" -^rel 


23 It is well-known that all odd virial coefficients, a 2 k+u are equal for bosons and fermions 
(Huang 1987). 
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where 


H c .m. — w — dzdz + wZ 9^ + wZ c^J , 

JU / 2 - 

H te i = 2 u> + >-9 W . 9 Wi +w«)i <9 Wt . + w w,- 9 Wi 

" V m 

t=l N 


(7.2.19) 


fYom now on, we focus only on the relative problem and neglect the trivial center 
of mass dynamics. In analogy with (7.1.22) we introduce the operators 


1 

'9 

II 

1 B 
o Wi , 

mu) 

4 

ii 

-to 

II 

s 

1 

i p. 

Ow\ i 

mu) 

Bi = d Wi 


(7.2.20) 


for i = 1,2, so that H re \ becomes 


H te l = U> 


2 + Y. ( A i + B i B <) 

i=l V 


The operators (7.2.20) satisfy the Heisenberg algebra 

[Ai, At] = [B h B]] = Sij , 


(7.2.21) 


(7.2.22) 


with all other commutators being zero, and are energy step operators. In fact, using 
(7.2.21) and (7.2.22), it is easy to prove that A,- and Bi lower the eigenvalues of 

-Hrei by u>, i. e. 


[-^reb Aj] — U) Ai , 


[Hteh Bi] = —U) Bi , 


(7.2.23a) 


whereas a\ and B\ raise the eigenvalues by w, i.e. 


[H teh A]} = u> A\ , [H tel , B}} = u B] 
Following (Sen 1991b,c), we define the operator 

Q = A\b 2 -A\ Bi 


(7.2.23b) 


(7.2.24) 


which is completely antisymmetric under particle exchanges, and commutes with 
H re !. Being antisymmetric, Q acts on symmetric wavefunctions to produce anti¬ 
symmetric ones with the same energy, and vice versa. Thus Q changes the statis¬ 
tics mapping bosons into fermions or vice versa, and for this reason it has been 
called a supersymmetry operator (Sen 1991b,c). More generally, Q connects staes 
of equal energy in two anyonic theories with statistics v and u + 1. In order to keep 
the statistics within the range [0,1], we can follow the action of Q with a parity 
transformation P which makes the replacements Wi «-+ Wi in the wavefunctions, 
corresponding to (xi,y») —> (#*,— yi) in terms of the cartesian coordinates. Since 
P changes the sign of the statistical parameter, the combined operator 
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Q = PQ (7.2.25) 

maps anyonic states of statistics v to anyonic states of statistics 1 — v without 
changing their energy. 

Let us now consider a wavefunction ^{u) describing the relative motion of three 

^ A» 

anyons of statistics u with a given energy Efy). Under the action of Q, in principle 
there are three possibilities: 

^ A 

i) The wavefunction Q'tp(v) is a regular wavefunction for three anyons of statis- 

A 

tics 1 — v and energy Efy); 

*+ A 

ii) The wavefunction Q ip(u) identically vanishes; 

A? A 

Hi) The wavefunction Q^{u) is singular. 

A detailed analysis is therefore necessary to see which one of these possibilities 
actually occurs. We refer the reader to the original literature (Sen 1991b,c) for a 
complete discussion; here we limit ourselves to summarize the results and see their 
consequences on the third virial coefficient. It turns out that if ^(u) is a regular 

A A 

wavefunction, then case Hi) never occurs, that is Q 'tp(v) is always an acceptable 
wavefunction (eventually vanishing) of a new anyonic theory with statistics 1 — v. 
Furthermore, the operator Q annihilates only all the Type-I and Type-II states 
which we constructed in Chapter 6 and which contributed to the partition function 
Z 3 {y , 0) in (7.1.29). The missing states (if they exist) are instead mirror symmetric, 

A 

that is for any regular missing wavefunction ^{u) there is a corresponding regular 

A 

missing wavefunction ^(1 — ^)* This property is confirmed also by recent numerical 
calculations (Sporre et aL 1991a; Murphy et al . 1991). 

Let us now assume that the missing states exist and construct the full three- 
body partition function according to 

Z 3 (u) = 

A 

{all states tp] 

= Z 3 (u, 0)+ J2 e-P E ^ . 

{missing states i p'} 

By construction we have 

Z 3 {u = 0) = Z 3 bos , 

and 

Z 3 {v = 1) = Z { 3 et , (7.2.27b) 

since the missing states are defined to make the bosonic and fermionic limits 
smooth. Taking into account the mirror symmetry of the missing states, we easily 
obtain 

Z z (u) - Z 3 (l - u) = Z 3 {v, 0) - Z 3 ( 1 - z/,0) , (7.2.28) 

which shows that only the regular Type-I and Type-II states contribute to the 

A? AT 

difference Z 3 (u) — Z 3 (l — u). Let us now define 


(7.2.27a) 


(7.2.26) 
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a${v) = 4a 2 (z/) 2 — 2 X\ lim 

(jj— 1-0 


3^1 -3Z 2 (u,0) + Z 2 1 


(7.2.29) 


which is the same expression as in (7.2.14) with Z$(v, 0) replaced by the full 
partition function Z 3 (z/), and then compute 


a>z{v)-a z (\ -u) . 


Using (7.2.28), we have 


a 3 (z/) - a 3 (l - v) = 4a 2 (z/) 2 -4a 2 (l - u) 


2 At lim 

1 w ->0 


Zi(v) - Z 3 ( 1 - u) 


-3[Z 2 (z/,0)-Z 2 (1-z/,0)] 


(7.2.30) 


= <z 3 (z/)-a 3 (l-z/) . 

From (7.2.15) with u f = 1 — z/, it immediately follows that 

a 3 M = a 3 (l-z/) . (7.2.31) 

The mirror symmetry of the third virial coefficient expressed by (7.2.31) is a rig¬ 
orous result which was first pointed out in (Sen 1991c). In particular for v = 0, 
we have the well-known property 

a*°° = a 3 (0) = a 3 (l) = a { 3 et = . (7.2.32) 

Unfortunately, the mirror symmetry (7.2.31) does not allow to deduce the value of 
a 3 (z/) for any v\ however, using (7.2.31) together with the perturbative calculations 
of (McCabe and Ouvry 1991), it can be conjectured that actualy the third virial 
coefficient does not depend on the statistics v . A proof (or a disproof) of this 
conjecture is still missing. 

For the higher virial coefficients, only partial pertubative results are available 
at the moment but no definite conclusions can be drawn (McCabe et al. 1991). 


7.3 Magnetic Moment 


In this section we are going to analyze some magnetic properties of the anyon gas. 
Let us consider a two-body system and define the magnetization according to the 
standard formula 

1 Q 

M 2 = -— log Z 2 • (7.3.1) 

As a consequence of the properties of Z 2 given after (7.1.28), it is easy to see that 
M 2 satisfies 
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M 2 (v,B) = —M 2 (2 — v , —B) 


(7.3.2) 


There are two regimes in which one can extract simple and interesting results. 
The first is the regime of high magnetic fields , in which B sets the biggest scale in 
the problem, namely 

m _ com „ _ 

« 1 , —— « 1 . (7.3:3) 


e/3B 


eB 


The first inequality tells us that the thermal excitations whose energies are VP 
are negligible with respect to the magnetic excitations whose energies are ~ co c = 
eB/m . Instead the second equality in (7.3.3) implies that the excitations induced 
by the harmonic force whose energies are ~ co are small compared to the magnetic 
excitations. In this regime M 2 reduces to a double series in the small parameters 
(7.3.3), whose first terms are given by 





B > 0 
B < 0 


(7.3.4) 


For u) = 0 the magnetic moment M 2 is independent of statistics; the second term 
on the right-hand side of (7.3.4) represents a finite temperature correction. As 
expected, the statistical dependence appears in the term proportional to it is 
in fact the harmonic potential that breaks the infinite degeneracy of the Landau 
levels of two anyons. The last term in (7.3.4) can also be reinterpreted as a “finite 
area correction” if we take into account the identification of co 2 with the inverse 
of the area as shown in (7.2.10). With such identification the last term of (7.3.4) 
represents a correction inversely proportional to the total flux BA through the 
sample. 

We remark that M 2 in (7.3.4) satisfies 


M 2 {—B, v = 1) = —M 2 (B, v = 1) 
M 2 (-B, v = 0,2) = -M 2 (B, v = 0,2) 


(7.3.5) 


which show that parity is not broken for bosons and fermions; instead for all other 
values of v parity is broken. 

The second regime which leads to an interesting result is that of low magnetic 
field 

^=^«i 

u) mu) 

This is the regime in which one can take the limit B —> 0 for finite values of u> and 
/?, and is the relevant one to study spontaneous magnetization effects. Indeed for 
low magnetic field, M 2 is given by a series in the small parameter (7.3.69) with 
functions of flu) as coefficients; more precisely we have 



M 2 (B, v) = jth[/?u>(l — z/)] coth(/?u>) — (1 — z/)| + O 


eB 


mu 


(7.3.7) 


For intermediate statistics (i.e. for v 7 ^ 0,1,2) there is a parity breaking sponta¬ 
neous magnetic moment p given by 
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fi = M 2 (B = 0,i/) . (7.3.8) 

This spontaneous magnetic moment vanishes in the limit of high temperature, 

lim fi = 0 , (7.3.9) 

0-0 

whereas at zero temperature it reduces to 

ill, = lim 0 = (7.3.10) 

0—oo 2 777. 

where L u is the orbital fractional angular momentum defined in Chapter 5 (see 
in particular (5.85) and (5.86) for N = 2). Notice that the gyromagnetic factor 
for fiz, is 1, as usual for an orbital magnetic moment. This is to be contrasted 
with the spin magnetic moment fis computed in (Kogan 1991; Stern 1991) for 
which the gyromagnetic factor is 2, as should be expected. Finally, we remark 
that the spontaneous magnetic moment fi vanishes for non-confined anyons (i.e. 
when to = 0) at any finite temperature. It would be interesting and important 
to extend these results to systems of many anyons, and see if the spontaneous 
magnetization we have evidentiated for two anyons, survives in the thermodynamic 
limit. Another interesting open problem is the study of the possible existence of 
a Curie-like critical temperature and a corresponding phase transition; however 
to make progress in this direction, it is necessary first to solve the many-anyon 
problem. 



8. Anyons and the Fractional Quantum Hall 

Effect 


The only known physical objects which can be described as anyons are the quasi¬ 
particle and quasi-hole excitations of planar systems of electrons exhibiting the 
fractional quantum Hall effect (QHE) (for a review see for instance (Prange and 
Girvin 1990)). Actually most of the great interest that anyonic theories have at¬ 
tracted in the past few years derives precisely from their relevance to a better 
understanding of the fractional QHE (Halperin 1984), in conjunction with several 
claims that anyons can provide also a non-standard explanation of the mechanism 
of high temperature superconductivity (Chen et al. 1989). It is far beyond the 
scope of these lecture notes to treat these issues in a systematic and adequate 
way, and therefore we refer the reader to the many good reviews already existing 
in the literature, for example (Wen and Zee 1989b; Arovas 1989; Lykken et al. 
1991). However, for the sake of completeness we think necessary to spend some 
time on at least one of these physical applications in order to convey the idea 
that anyons are not just mathematical fantasies. Since recent experiments have 
cast some shadow on the relevance of fractional statistics to the observed high 
temperature superconductivity (Lyons et al. 1990; Kiefl et al. 1990; Spielman et 
al. 1990), here we will concentrate only on the application of anyons to the theory 
of the fractional QHE. Of course our presentation will be schematic and not at all 
exhaustive. Furthermore, we will concentrate more on the formal aspects than on 
the condensed matter issues. 

Let us recall that the QHE is observed in two dimensional systems of electrons 
at very low temperatures (a few degrees Kelvin or even less) and in very strong 
magnetic fields (~ 10 Tesla or even more) orthogonal to the plane where the parti¬ 
cles move. The electrons are usually trapped in a thin layer at the interface between 
two different semiconductors or between a semiconductor and an insulator. The 
low temperature and the strong magnetic field freeze the motion along the direc¬ 
tion perpendicular to the layer, and thus the relevant dynamics takes place on a 
plane; moreover the effects due to the electron spin can be neglected in a first 
approximation. The QHE is characterized by the fact that the Hall conductance, 
<jjj, is quantized in units of e 2 /h (—e is the electron charge and h is the Planck 
constant) 

os.=v e -^ . ( 8 . 1 ) 

The quantum number v can be an integer (integer QHE) (von Klitzing et al. 1980) 
or a simple fraction (fractional QHE) (Tsui et al. 1982). Quantization rules are 
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usually associated to symmetry considerations, like for example the quantization 
of spin which is related to the algebra of rotations. The quantization rule of the 
QHE seems instead to have a different origin, at least in the integer case; it is a 
topological quantization (Thouless et al 1982; Avron et al. 1983; Niu et al. 1985; 
Kohmoto 1985), like the one of the magnetic monopole charge (Dirac 1931). How¬ 
ever, contrary to the latter which has not been observed so far, the quantization 
of an is a well-established fact. Actually, what is astonishing and fascinating is 
that (8.1) is experimentally observed with extremely high accuracy (the precision 
is roughly 10 -7 ~ 10 -8 for the integer QHE, and 10 -4 ~ 10 -5 for the fractional 
QHE). 

The external magnetic field organizes the energy spectrum of the electrons into 
Landau levels (see for instance (Landau and Lifschitz 1977) and Chapter 6) and 
forces the particles to fill such levels from bottom up. A quantity which plays a 
central role in the QHE is the filling factor v 24 , defined as the number of electrons, 
iV, per number of Landau levels available. Each Landau level is highly degenerate 
(see (6.1.13)) but for samples of finite area A, this degeneracy is finite. The number 
of available levels is in fact given by 


where is the magnetic length 


l 

Therefore the filling factor is 

v = - 

j 

where p = N/A is the electron density. Eq. (8.4a) can also be written in a different 
way using the explicit definition of the magnetic length given in (8.3), namely 




) — 


[he 

vTb 


N 


i/27 xt\ 


= 2*§P 


( 8 . 2 ) 


(8.3) 


(8.4a) 


N 

ABe/hc 



(8.4b) 


where <f> = A B is the magnetic flux through the area A, and <f >o = h c/e is the flux 
unit. It is not difficult to prove that the filling factor (8.4) is the quantum number 
appearing in (8.1). Here we sketch a possible “heuristic” proof of this fact. Let us 
suppose that the electrons move in the (x, t/)-plane with a velocity v x . Since the 
magnetic field is directed along the z-axis, the electrons feel a Lorentz force in the 
^-direction 

F y ~ -v x B , (8.5) 

c 

which is then compensated by an emerging electric field E y , such that 


24 


We will show in a moment that the filling factor is the quantum number appearing in 
(8.1); this is why we used the same symbol. 



Ill 


e Ey ~ -v x B . (8.6) 

c 

The transverse Hall conductance is defined as 

= ®xy ~ ^=T ( 8 * 7 ) 

where j x ~ epv x is the electric current in the x-direction. Using (8.6), we obtain 

epv x ec 
aH ~ v x B/c ~ B P 

which is a well-known classical result. Using in (8.8) the quantum mechanical 
relation (8.4) to express p in terms of the filling factor v , we get 

e 2 

an ~ v n 

as in (8.1). We refer the reader to (Prange and Girvin 1990) for a rigorous deriva¬ 
tion of this formula and a discussion of the fundamental role of impurities. 

When v is integer (integer QHE) there is an integer number of Landau levels 
completely filled, and the phenomenon becomes essentially a direct manifestation 
of the Landau quantization for “non-interacting” electrons in a magnetic field. 
Furthermore, the system is incompressible and develops a gap. The situation is 
more complicated and certainly less obvious in the fractional QHE that results 
from the condensation of the two-dimensional electron system into a new type of 
collective ground state driven by the Coulomb repulsion. This is incompressible 
and exhibits a gap like in the integer QHE, despite the fact that only a fraction of 
Landau levels is filled. 

When the fraction is v = 1/m, with m an odd integer, this ground state can 
be described very accurately by the variational wavefunction proposed by R.B. 
Laughlin (Laughlin 1983) 

„ -^£m 2 

JJ(*-*;) ra e 4€ ° ' 

I<J 

where zj is the complex coordinate for the /-th electron and M m is a normalization 
factor. Notice that the wavefunction (8.9) has precisely the same form which we 
found in Chapter 4 for particles of statistics 

v — m . (8.10) 

Since m is an odd integer, -0m is totally antisymmetric, and so it describes ordinary 

fermions (or better “super”-fermions if m = 3,5,...). 

Before exploring how anyons may enter the stage of the fractional QHE, let 

us examine some properties of the Laughlin wavefunction (8.9). First of all, the 

prefactor J~[ (zj — zj) m is purely analytic, which means that all particles are in 

I<J 

the lowest Landau level (see Chapter 6). However, is not simply the product of 
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single particle wavefunctions, but is a complicated superposition of such products. 
The prefactor is also of the Jastrow type: it has a zero of order m at coincident 
points (zj = 2 j), showing that electrons tend very strongly to repel each other in a 
way that is appropriate to minimize the Coulomb interaction. In this sense they are 
called sometimes “super”-electrons (with no reference at all to supersymmetry!). 
If zj goes around zj by an angle A(p, the wavefunction acquires a phase e imA<p ^ as 
if each particle carried m units of flux. 

Many of the phenomenological consequences of the Laughlin wavefunction de¬ 
rive from the so called plasma analogy. This is based on the interpretation of the 
quantum probability density, |as a classical Boltzmann distribution with 
potential energy V m , according to 

IV’ml 2 = e~ Vm (8.11) 


where 


V m = -2m ^ ln I z i ~ Z J 

KJ 



( 8 . 12 ) 


The logarithmic terms of V m clearly derive from the prefactor of (8.9), while 
the last term originates from the exponential 25 . It is well-known that a two- 
dimensional one-component plasma of particles with charge q moving in a neutral¬ 
izing background of density />, is characterized by a potential 



(8.13) 


where the first term is due to the Coulomb interaction among the particles, and 
the second term is due to the interaction of the particles with the neutralizing 
background. Eq. (8.13) has the same form of (8.12) and so we are naturally led to 
regard V m as the potential of a plasma of particles with charge 

q = v2m 

in a neutralizing background of density 

_ 1 _ _L 

Variational arguments suggest the identification of p in (8.15) with the electron 
density of the original problem, and therefore from (8.4) we conclude that the 
wavefunction describes a (circular) droplet of fluid with filling factor 


o 


1 _ 

m 


(8.14) 


(8.15) 


v = 27t^q P = — * (8.16) 

772 

There are very strong arguments which support the identification of in 
(8.9) as the wavefunction for the ground state of the fractional QHE at filling 


25 


In (8.12) we have neglected the normalization factor which would give rise only to 
a trivial constant shift of Vm. 
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v = 1/m, but essentially the justification is that has overlap one (with very 
high accuracy) with the numerical solutions for the QHE computed for various 
kinds of repulsive potentials (Laughlin 1983). This shows that captures the 
essential and universal features of the fractional QHE since it is largely insensitive 
to the specific form of the repulsive interactions among electrons. 

Any deviations from the density p in (8.15) are accommodated in the system 
as localized quasi-particle or quasi-hole excitations across a gap. These excitations 
turn out to have fractional charge and fractional statistics; thus they are anyons! 
To prove this, let us begin by recalling how such excitations can be obtained. Let us 
consider the following ideal situation: in the state described by (8.9) we introduce 
an infinitesimally thin flux-tube at a point z a , and then turn on adiabaticaily the 
flux <j> from zero to a final value of one unit ( i.e . <j> = 0 —► <f> = ±<f >o = c/e), in 
such a way that the system remains an (instantaneous) eigenstate of the changing 
Hamiltonian. Due to Faraday law, the variation of the flux from <j> = 0 to <j> = ±<f >o 
gives rise to a (circular) electric field around the point z a . The particles will then 
flow inwards or outwards (depending on the sign of <£), and a net positive or 
negative charge will accumulate around z a . However, since the change of flux by 
one quantum <j >o can be compensated with a gauge transformation (Wu and Yang 
1975), the final state can be considered an excited state of the original Hamiltonian. 
The excitation is a quasi-particle or a quasi-hole and its charge turns out to be 
e* = ^=e/m. In other words, one can say that a quasi-hole (or a quasi-particle) is 
formed in the incompressible fluid described by as a bubble of size such that 
a fraction 1/m of an electron is removed (or added). 

For definiteness let us begin by considering the case of a quasi-hole. It is not 
so difficult to guess that the wavefunction for the state at filling v = 1/m with a 
quasi-hole at the point is of the form 

= M+Y[(zi - z a ) 4> m (8.17) 

I 

where A/+ is a normalization factor and is the ground state wavefunction (8.9). 
Indeed, when the /-th electron travels along a path surrounding other particles, it 
sees m units of flux at all zj with J ^ I where the other electrons are located, but 
when it travels around the quasi-hole in it sees only one unit of flux. Thus, as 
expected, the quasi-hole behaves like 1/m of an electron and in particular carries 
1/m of its charge (of course with opposite sign!). These conclusions can also be 
obtained by applying the plasma analogy directly to (8.17) (Laughlin 1983). 

The wavefunction describing a quasi-particle excitation above the ground state 
■0 m is more complicated and in a certain sense less obvious than that of the quasi¬ 
hole. However, also in this case it is possible to write the wavefunction in analytic 
form (Laughlin 1983), namely 

= N- n (*8£ - *«) i’m . (8.18) 

Here £ 0 is the magnetic length (8.3) and it is understood that all derivative opera¬ 
tors do not act on the exponential factor of Considerations similar to those we 
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have mentioned above for the quasi-hole, allow to conclude that the quasi-particle 
has a charge e* = —e/m (notice that the electron charge is —e). 

The charge of the excitations of the fractional QHE can be computed also in a 
more direct but more formal way using the concept of Berry phase (Berry 1984), 
an approach originally advocated by D.P. Arovas, J.R. SchriefFer and F. Wilczek 
(Arovas et al 1984). This method is very elegant and will allow us to compute also 
the statistics of the excitations and conclude that they are indeed anyons (cf also 
(Forte 1991)). For simplicity we begin by considering the case of the quasi-holes 
wavefunction (8.17). In the Berry phase approach, the charge of the quasi-hole is 
determined by calculating the change in the phase of the wavefunction as the 
quasi-hole position z a is adiabaticaily moved around a closed loop encircling a flux 
<f >. This phase is then compared and identified with the Aharonov-Bohm phase (see 
(3.2.44)), and the charge is determined accordingly. Before entering the details of 
this calculation, let us briefly recall the concept of Berry phase (extensive reviews 
as well as original research papers on this subject can be found in (Shapere and 

Wiczek 1989)). 

Let us consider a quantum system described by a time-dependent Hamilto¬ 
nian H(t). If we assume that at a given time to the system is described by an 
instantaneous eigenstate ^ a (to), satisfying 


H(t) 1p a (t) = E a (t) 1p a (t) , 


(8.19) 


then in the adiabatic approximation (see for instance (Schiff 1968)) the system at 
time t\ will be represented by 


*M*i) = exp |-^/ dt E a (t) -i ip a (ti) 


( 8 . 20 ) 


where the “unusual” phase 7 is 


ICC 


- -* /“ 

Jto 


dt . 


( 8 . 21 ) 


More generally, if the system at time to is a superposition of instantaneous eigen¬ 
states V>a(t), i-e. 

*(<o) = J3<*aV’a(<o) ’ ( 8 - 22 ) 


a 


then in the adiabatic approximation the wavefunction at time i 1 is 

&(ti ) = ]C exp [ - ^/ dtE a (t)~ i a a i> a {ti) 


(8.23) 


The phase 7 ^ in (8.21) is a geometric phase whose relevance in physics has been 
recognized only recently (Berry 1984). In fact it is usually unobservable (Schiff 
1968), but in some cases, when the wavefunctions are not single-valued, it has 
interesting and observable consequences. Anyons, whose wavefunctions are multi¬ 
valued (see Chapter 4), are therefore the obvious candidates to test the importance 
of Berry phase. 



115 


Let us now return to our problem and compute the Berry phase for the wave- 
function (8.17) under the assumption that the quasi-hole is slowly transported 
around a closed loop T so that z a becomes a time-dependent parameter and the 
adiabatic approximation can be used. Since the entire time-dependence of the 
problem comes solely from z a , using the explicit expression (8.17), we easily ob¬ 
tain 

d (t) = E J t N 2 ' - *«(*))] (t) , (8.24) 


dt 


and 


7 


= ~ 1 Jt dt ^ ^ _ Za ^ ^ 


m 


(8.25) 


The integral over t in (8.25) is evaluated from the initial time to when the quasi¬ 
hole begins its tour, to the final time ti when it returns to its original position 
after a non trivial loop T traversed anticlockwise. We can rewrite (8.25) in a more 
transparent way by observing that the electron density in the state , is (see 

(5.78)) 

p( z ) = (V&*"(0I Yl 6(2) ( z ~ iV’m*" (*)> (8-26) 

I 


so that 


7 


= -i f dt f d 2 z -^[ln (z-z a (t))] p(z) 

J to J 


= -‘/ 


d 2 



*0 


u dt dz<x ® 1 

dt z a (t) — z 


p( z ) 


(8.27) 


d 2 



dz a 


z a -z 


p( z ) 


Under the reasonable assumption that the density p(z) is a regular function (see 
for example (Arovas 1989)), it is easy to realize using the residue theorem that all 
points z outside T give a vanishing contribution to the contour integral, and thus 
are left with 

= —i / d 2 z <p dz a - p(z) 

J<r Jr z <x ~ z (8.28) 


7 


= 27 t f d 2 z p(z) = 2irNr . 

J<r 


The symbol f <r d 2 z denotes the integral over the region inside the loop T, and Nr 
the average number of electrons that are in there. This concludes the calculation 
of the Berry phase for the quasi-hole wavefunction (8.17). 

In Chapter 3, we have seen that if a particle of charge q is moved along a closed 
loop T encircling a flux <j>r , then its wave function acquires a phase 


exp 


—l 


. q<t>r 


he 


(8.29) 


(see (3.2.44)) due to the Aharonov-Bohm effect (Aharonov and Bohm 1959). If we 
identify the Berry phase 
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exp [-i 7 ] 

with the Aharonov-Bohm phase (8.29), we have 

7 = IpL = 271 -Nr 
n c 


(8.30) 


(8.31) 


On the other hand, using (8.4b) and (8.16), the flux <j>p turns out to be 


<f>r = 


Np<f> 0 hcNp 


m 


v 


(8.32) 


inserting this into (8.31) and solving for < 7 , we finally obtain 


hcN P 

<f> P 



(8.33) 


Thus we have proven that the quasi-holes have charge q = e* = e/m. 

The calculation of the charge of a quasi-particle proceeds along the same lines, 
but particular care must be used in treating the derivative operators which appear 
in the wavefunction (8.18). Formally we can write 


I 




Tpm Za (t) > 


(8.34) 


where we have to remember that the derivatives d/dzj do not act on the expo¬ 
nential factor inside Upon using (8.34) into (8.21), the Berry phase for a 

quasi-particle becomes 



We stress that (8.34) and (8.35) are only formal expressions whose precise meaning 
needs to be clarified and specified in the following. This can be accomplished by 
defining the multi-particle generalization of the Bargmann-Fock inner product (see 
for instance (Girvin and Jach 1983)) according to 

<<5|0 ( {M «S7 ;2r) ) = 

1 j | 2 

jY[d 2 zje~ W ° ^ " g({ Zl })0^{2ll^--,zjyjf({ Zl }) . 

1 (8.36) 

Here / and g are two arbitrary analytic functions of many variables, and O is an 
operator, functional of 21q d/dzj and 2 /, acting on the space of analytic functions. 
If such operator is normal-ordered, that is if all derivatives d/dzj are put on the 
left of all complex coordinates 2 /, then it is not difficult to prove the following 
identity (Girvin and Jach 1983) 
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«<?l O ({2^;*/}) l/» = «<?l 0({zj;zj})\f}} . (8.37) 

Indeed, any derivative 2 l^d/dzi acting from the right on the measure factor 



of (8.36) yields simply z/, and therefore inside the Bargmann-Fock inner product 
we can safely perform the replacement 



Using the definition (8.36) and remembering that the derivatives d/dzj do not act 
on the exponential factors, it is easy to realize that the integrand of (8.35) can be 
rewritten as follows 



where 'ip rn Za is the wavefunction of the quasi-particle without the exponential fac¬ 
tor, that is 

415-El^l 2 

Vw‘“(0 = e °' VC‘“(0 • (8.40) 

Since the formal operator 

^ ! ['“<2^ - M<)) 

is normal-ordered according to our previous definition, upon using (8.38), the Berry 
phase becomes 

7 = -i f* 1 dt i - urn\fa-m • (8.4i) 

As in the case of a quasi-hole, 7 can be written in a more transparent way using 
the electron density />(z), so that 





(8-42) 
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Comparing this with the Aharonov-Bohm phase (8.29), we deduce that the quasi¬ 
particle have charge 



e 

ra 


(8.43) 


as expected. 

As we have seen, the calculation of the charge of the excitations above the 
ground state presents essentially no difficulty, and both the plasma analogy and 
the Berry phase approach give the same results. On the contrary, the issue of what 
statistics these excitations obey is more subtle, and indeed all possible statistics 
(fermionic (Laughlin 1983), bosonic (Haldane 1983) and anyonic (Halperin 1984)) 
have been proposed. In the following we are going to show that the quasi-particles 
and quasi-holes are anyons of statistics v = 1 /m; the calculation is once again 
based on the Berry phase 7 (Arovas et ah 1984). For simplicity and to avoid 
repetitions, we will examine only the case of quasi-holes. 

Let us then consider a state characterized by the presence of two quasi-hole 
excitations, one located in the point and the other located in the point zp. If 
we neglect any interaction between the two quasi-holes, the wavefunction turns 
out to be 

i>m a ,+Z0 = N a p n ( 2 -r - 2 «) ( z * - z p) V’m ( 8 . 44 ) 

I 


which is clearly a straighforward generalization of (8.17). In (8.44) is a nor¬ 
malization factor that is irrelevant for our present discussion. Now let us assume 
that the quasi-hole in is adiabatically moved on a full loop F, while the quasi¬ 
hole in zp is kept fixed. During this process the wavefunction ipti 0 ** 20 acquires 
a Berry phase 7 given in (8.21) that can be calculated as before. Skipping the 
intermediate steps which are exactly the same as in (8.24-27), we have 




<01-^ 


+ Z a ,+ Z0 

m 




(8.45) 


where p(z) is the electron density in the two quasi-hole state (8.44), namely 


p(z) = (^ Za ' +Z0 (t) I X> (2) (* - * 1 ) I• (8-46) 

I 


If in its motion z a does not encircle zp , the integral in (8.45) has the same value 
as in (8.28) so that 

7 = 27t Nr . (8.47) 

Of course this case corresponds to the situation in which the two quasi-holes are 
not exchanged. Indeed, as far as the quasi-hole configuration space is concerned, 
the loop r is homotopically trivial, according to our general discussion of Chapter 
2, because z a does not encircle zp . However, if the loop r traversed by contains 
zp and it is homotopically non-trivial, we have 
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7 = 27r J d 2 z p(z) = 27r ^ . 


(8.48) 


In this case in fact we have to take into account that a quasi-hole is present inside 
r and hence the number of electrons Np has to be diminished by the fraction 1/m 
that is needed to build the quasi-hole. Comparing (8.48) and (8.47) we see that 
when z a encircles zp the wavefunctions picks up an extra phase 


exp [—i Ay] = exp 



(8.49) 


which can be interpreted as a statistical effect. Indeed in this case the two quasi¬ 
holes are exchanged twice, and from ( 2 . 1 ) with A<p = 2n we conclude that their 
statistics is 



(8.50) 


When m = 1 , i.e. when the lowest Landau level is completely filled (see (8.16)), 
the quasi-holes are fermions, but in general for m = 3,5,... the quasi-holes are 
anyons of fractional statistics. A similar result can be obtained also for the quasi¬ 
particles using the replacement (8.38) inside the expectation values which express 
the Berry phases. 

The derivation of (8.50) clearly shows that the statistics v is directly related 
to the fraction of electrons which forms a quasi-hole or a quasi-particle. As we 
have mentioned, this fraction is a well-defined number as a consequence of the 
incompressibility of the ground state (Laughlin 1983); in particular it does 
not depend on the details of the loop T but only on its homotopy class so that 
(8.49) can be interpreted as a statistical effect. Furthermore, comparing (8.16) 
and (8.50), we observe that the excitations above the Laughlin ground state are 
characterized by a statistics which is numerically equal to the filling factor of 
Vw This observation will be helpful in deriving the so-called hierarchy scheme 
(Halperin 1984) which is a systematic way to understand the observed fractions 
in the QHE and is based on the explicit use of anyons. 

We have seen that the quasi-hole excitations above ipm have charge e* = e/m 
and statistics v = 1 /m; therefore it is natural to expect that the effective wave- 
function for two such quasi-holes with coordinates z a and zp should contain terms 
like 

V'm, ~ (*a - ^) mi e'^ (kl +M ) (8.51) 


where 


mi =-(- 2pi 

m 


(8.52) 


with pi a positive integer. In fact, under the exchange of the two quasi-holes the 
factor (z a — zp) mi produces the required phase 


exp[i7rmi] = exp 


i7r(— + 2pi) 

= exp 

r. ii 

17 r — 

m 


m 
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The integer p\ has no effect on the statistics since the latter is defined only modulo 
2 and when p\ ^ 0, the quasi-holes of if> mi are called “super-anyons” in analogy 
with the electrons of 'tpm which are sometimes called “super-fermions” when m = 

3,5,_The exponential term in (8.51) is typical for objects of charge e* = e/m, 

interacting with a magnetic field field B; indeed the magnetic length for the quasi¬ 
holes with charge e* is _ 


£ o = )j^B = '/™ l o 

where Iq is the electron magnetic length (see (8.3)). The effective wavefunction 
(8.51) for the quasi-holes should result from the wavefunction (8.44) after the 
integration of all coordinates zj of the underlying electron system, and can be 
easily generalized to an arbitrary number of quasi-holes according to 



V’m, = ^4 n ( Za - Z p) mie 

ol<(3 



(8.54) 


This is clearly the straighforward generalization of the Laughlin wavefunction 'tpm 
which described an incompressible system of electrons with charge —e, statistics 
v = 1 and filling factor v = 1/m. Using the same arguments, it is therefore easy 
to conclude that in (8.54) describes an incompressible system of quasi-holes 
with charge e* = e/m and statistics v = 1/m. Let us now compute the new filling 
factor by applying the plasma analogy to . We interpret |^+ I | 2 as a classical 
Boltzmann distribution with potential 


v mi = -2mi ln 



(8.55) 


Comparing this expression with (8.13) we immediately find that the density of the 
neutralizing background is 

P* = - • (8.56) 

27 r£$ mrrii 

As for the electrons, variational arguments again indicate that p* in (8.56) has to 
be identified with the density of quasi-holes; therefore in an area 27 t£q there are 

2 _ _ 1 


2*/$ p m = 


mm\ 


quasi-holes carrying a total charge 


?i = 


m Vmmi 


(8.57) 


This charge has to be removed from the underlying electron system; hence the 
electron filling factor changes from v = 1/m to 


V\ = v — 

_ 

m 


<h 

e 


m* mi 


(8.58) 
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Notice that —qi/e is precisely the number of electrons required to produce a charge 
qi in the area 2^1^. Elementary algebra leads to 



2pi 

1 + 2pi 777 


1 


771 + 



(8.59) 


We stress once again that this is the electron filling factor resulting from the 
formation of an incompressible system of quasi-hole excitations above the ground 
state at filling 1/771. It is easy to see that if 771 = 3 and p = 1, then 

2 


which is one of the experimentally observed fractions in the QHE. 

Notice however that the wavefunction 1 with m = 3 and pi = 1 does 
not describe the electron state at filling v\ = 2/7 since it is a wavefunction for 
quasi-holes! However it describes a new incompressible quasi-hole fluid with its 
definite filling factor which in turn is related to the effective filling factor of the 
parent electron system. Moreover the fractional statistics of the quasi-holes is in 
correspondence with this effective filling factor in the same way as the ground state 
filling factor v = 1 /m (with m an odd integer) was related to the fermi statistics 
of the electrons. 

This construction and these observations can be further generalized and ex¬ 
tended in a hierarchy scheme (Halperin 1984; Haldane 1983) by considering exci¬ 
tations over excitations in a recursive way. For example one can construct second- 
generation quasi-hole excitations from the quasi-hole wavefunction ip + li . These 
are characterized by a charge e** = —e*/m\ = —e/(mmi), and by a statistics 
v\ = 1 /77ii • Thus their effective wavefunction would have the form 


V’mt = -V++ II ( Za - Z p) m2e 

OL<(3 


1 


4m m\t 


3 

0 


Ei*«i a 


(8.60) 


where 

777-2 =-("2p2 (8.61) 

777 1 

with P 2 being a positive integer. By repeating the same steps as before, one may 
compute the net filling factor for the electron system: this is given by the original 
one, v = 1 / 777 , diminished by the proper amounts necessary to produce the two 
generations of quasi-holes, namely 


1 1 

v 2 = v 1 - 

777 777 i 777 777i m 2 


1__I_ 

777 TO 1 777 1 J 777 2 777^ 7772 


(8.62) 


Elementary algebra leads to 
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(8.63) 


Clearly this hierarchy construction can be further iterated, and when also the 
quasi-particle excitations are taken into account, the generic filling factor is rep¬ 
resented by a continued fraction (Haldane 1983) 



771 + 


a i 


2p\ + 


a 2 


2p 2 + 


<*3 


2pz + 


(8.64) 


Here a* = +1 if the i-th generation consists of quasi-holes (with respect to their 
parent fluid) and a* = —1 if the i-th generation consists of quasi-particles. In this 
way all observed fractions are reproduced. Unfortunately however, this hierarchy 
construction predicts also fractions which are not experimentally observed. To 
avoid this problem, other approaches have been proposed (Jain 1989a,b, 1990) 
and we refer the reader to the original literature for their discussion since they do 
not make an explicit use of the idea of anyons. 


9. Anyons and Conformal Field Theories 


In the past few years, two-dimensional conformal field theories (CFT) (Belavin et 
al 1984) have attracted an enormous interest in many different contexts. In fact, 
CFT describe the scaling limit at the critical point of many statistical mechanical 
systems in two dimensions (for reviews see for instance (Ginsparg 1990; Cardy 
1990)); they are relevant in string theory because they define the background ge¬ 
ometry in which the string propagates (see for example (Green et al 1987)); they 
are essential in the non-perturbative formulation of two-dimensional gravity pro¬ 
vided by matrix models (for reviews see (Bilal 1990; Ginsparg 1991)), and last but 
not least, CFT have a very rich and elegant mathematical structure (Moore and 
Seiberg 1989). The latter derives essentially from the fact that CFT are quantum 
theories which carry a non trivial representation of the two-dimensional conformal 
algebra. This is known also as Virasoro algebra, and reads as follows 

£ 

[L n i L m ] = (^ ^) L n _|_ m “b ^n+m,0 i (9.1) 

where n, m are integers and c is the so-called central charge. For a discussion 
of the properties of the Virasoro algebra and of its representations, we refer the 
reader to the extensive literature on the subject (see for example (Itzykson et al 
1988) where the most relevant papers on CFT are collected). Here we simply recall 
that in a CFT the iV-point correlators of fields describing local observables, are 
real functions of the complex coordinates where the fields are located, but can 
be decomposed into a sum of factorized terms, each one being a product of a 
purely holomorphic factor times a purely anti-holomorphic one. These factors are 
called conformal blocks and carry non-trivial representations of the braid group 
Bn (Moore and Seiberg 1989). 

The appearance of the braid group naturally leads to conjecture some corre¬ 
spondence between anyons and CFT. A first signal of such correspondence already 
appeared in Chapter 3, where we pointed out the presence in anyonic systems of 
the conformal symmetry generated by L±\ and Lq, which indeed close the finite 
sub-algebra of (9.1). 

In this chapter we will elaborate on the formal analogy existing between anyons 
and CFT, and show that the anyonic wavefunctions described in Chapter 6 - in 
particular those which are relevant for the fractional quantum Hall effect - can 
be represented as conformal blocks of a particular two-dimensional CFT. The 
literature on this subject is quite extensive (Fubini 1991; Stone 1991a; Fubini and 
Liitken 1991; Moore and Read 1991; Cristofano et al 1991a,b,c; Balatsky 1991; 
Dunne et al 1991a,b; Balatsky and Stone 1991; Nagao 1992a,b; Ting and Lai 
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1992; Cappelli et aL 1992) but still many important and fundamental issues have 
to be elucidated, in order to really understand the significance of this analogy. In 
particular it would be nice to exploit the vast knowledge reached in CFT to make 
further progress in the solution of the anyonic theories, and also to understand 
the incompressibility of the Laughlin wavefunctions at rational filling factors (see 
Chapter 8) as a direct consequence of the conformal symmetry. 

To be specific, let us consider the case of N identical and indistinguishable 
particles interacting with an external magnetic field and confined in the lowest 
Landau level. The single-particle wavefunctions are given by (6.1.12) with n = 0, 


t.e. 



where £ 0 is the magnetic length (cf (8.3)) and j = 0,1,2,... denotes the angular 
momentum of From now on, we will remove the exponential factor from the 
wavefunctions and absorb it into the integration measure whenever we have to 
compute matrix elements or overlap integrals. Thus, the “reduced” wavefunctions 
for the lowest Landau level are simply the powers zK To discuss the iV-particle 
states we use the formalism of second quantization under the assumption that the 
dynamical variables of the single-particle states suffice also to describe a collection 
of many identical particles, even when these interact with an external field or 
among themselves (including the case of statistical interactions). Thus, for each 
single-particle mode with j = 1,2,... 26 , we introduce a creation operator dj, 
and an annihilation operator dj, satisfying the Heisenberg algebra 



with all other possible commutators being zero. Then, we define a bosonic quantum 
operator Q(z), functional of the single-particle modes, as follows 


a 1 


Q(z) = Qo(z)+i z 3 

m 3 


(9.4) 


where 

dj = d_j for j < 0 . (9*5) 

Qo(z) in (9.4) is a zero-mode part which we specify in a moment. Since dj is 
a creation operator for j > 0 and an annihilation operator for j < 0, the Fock 
vacuum |0) is such that 


at |0> = 0 

and 

<0| at = 0 


for j < 0 , 

(9.6a) 

for j > 0 . 

(9.6b) 


26 The state with j = 0 has the lowest angular momentum and can be regarded as the 
“vacuum” from which all other states are generated by step-operators. Thus, the state 
with j = 0 is treated differently from all other single particle states. 
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On the other hand, Q(z) is meant to be a quantum bosonic scalar field in two 
dimensions, and hence its propagator should be 

(0| Q(z) Q(w) |0) ~ ln(z - w) . (9.7) 

This requirement fixes Qo(z) to be 

Qo(z) = x — ip In z (9.8) 

where x and p are zero-modes which satisfy 

[x , P] = i , (9.9) 

and 

p|0) = 0 , (0| x = 0 . (9.10) 

Eq. (9.10) implies that p can be considered an annihilation operator, and x a 
creation operator on the Fock vacuum |0). 

The bosonic field Q(z) given by (9.4) and (9.8), is used in the operator for¬ 
malism of string theory to describe the string coordinates (for a review see (Green 
et al. 1986)), and also defines a CFT with central charge c = 1 whose Virasoro 
generators are given by 


L n = \ ^ ■ ■ (9.H) 

j =—oo 

In this formula, do = dj = p, and the symbol : : denotes the normal ordering in 
which, as usual, the creation operators are put on the left and the annihilation 
operators on the right. 

Now we wish to show how to use Q(z) to recover the multi-anyon wavefunctions. 

To this end, we have to keep in mind that the braid group properties of the anyonic 

wavefunctions cannot be enforced by hand, contrary to the bosonic and fermionic 

cases where the appropriate symmetrizationsor antisymmetrizations can always be 

performed without any problem. Thus, the information about the anyonic statistics 

has to be built in the second quantized formalism itself. However, the field Q(z) 

defined by (9.4) and (9.8) is bosonic in nature and cannot directly be the relevant 

operator to discuss anyonic statistics. Furthermore, a simple product of N fields 

Q(zj) is certainly inadequate to reproduce the anyonic wavefunctions 27 which, 

as we have seen in Chapter 6, are characterized by a prefactor n (zi - zjY* 

I<J 

These difficulties can be easily overcome if we introduce the so-called Fubini- 
Veneziano vertex operators which are essentially the exponentials of Q(z) (Fubini 
and Veneziano 1970). More precisely, to describe anyons of statistics v we define 
the vertex operator 

V u {z) = : ? (9. 12 ) 

27 To avoid repetitions, we limit ourselves to discuss Type-I wavefunctions. The Type-II 
wavefunctions can be treated in a similar way if z *-+ z and v «-* 2 — v. 
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where : : denotes the normal ordering also of the zero-modes (i.e. x is put on the 
left and p on the right). If we insert in (9.12) the definition of Q(z) and perform 
explicitly the normal ordering, we find 


V„(z) = exp 



(9.13) 

The Fubini-Veneziano vertex operators have been extensively used in string the¬ 
ory to compute string scattering amplitudes (Green et aL 1987), and in CFT to 
compute conformal blocks (Ginsparg 1990). Here, following (Fubini 1991; Fubini 
and Liitken 1991; Dunne et aL 1991a,b) we show that the vertex operators can be 
used also to compute wavefunctions for anyons in the lowest Landau level. 

First we note that by a simple reordering, the product of N vertex operators 
can be rewritten as 



n v v( z i )=n - zj ) 

7=1 I<J 


exp 



N 

x exp(iv^iVx) JJ Zj”** exp 

7=1 



(9.14) 


This has the immediate consequence that the vertex operators are good candidates 
for anyonic operators in that they satisfy anyonic exchange properties. In fact, one 
has 

V„(z 1 )V„(z 2 ) = V„(z 2 )V„(z 1 )e**''* , (9.15) 

where the signs in the exponent correspond to the two topologically distinct ways 
of interchanging z\ and z<i (t.e. clockwise and counterclockwise). Moreover, the 
vertex operator (9.12) is a primary field of weight v j2 in the CFT whose Virasoro 
generators are given in (9.11); this means that 


[L n , V v (z)\ = z n 



+ £(» + !) 



(9.16) 


for any n. In particular, (9.16) for n — 0 implies that the vertex operator behaves 
as a fractional spin operator under rotations, namely 

e iLo9 V„(z) e ~ lL °0 = V„(ze iff ) , (9.17) 

where the rotation generator Lq has the oscillator representation given in (9.11) 
and $ is the rotation angle. From (9.17) we can read the value of the spin that is 

s= V - , (9.18) 

in agreement with the spin-statistics connection. 
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Now we wish to show that the expectation values of products of vertex oper¬ 
ators, which in CFT represent conformal blocks, coincide with the reduced anyon 
wavefunctions. To this end, let us define the state |p) and its dual ((p| such that 


p\p)=p i p) , 

aj- Ip) = 0 for 
((p| aj- = 0 for 

((p|p') = 6(p - p') 


((p\ P = ((p\ P , 

j< 0 , 

j > 0 , 


(9.19) 


Using (9.9), it is straightforward to prove that 


(H) |p) = |p + i) , ((p| (ix) = {{ P - 1 | . ( 9 . 20 ) 

Clearly the state |p = 0) is the Fock vacuum |0) (cf (9.10)), but the state ((p = 0| 
is not the dual vacuum (0|. In fact, the latter is annihilated by x, while ((p = 0| is 
annihilated by p. However, we can easily establish the precise relationship between 
(0| and the states ((p| if we interpret x and p as a standard pair of canonically 
conjugate operators 28 . Let us define the state (x| such that 

(x| X = {x\ X , 

then using the Fourier transformation and going to the momentum representation, 
we can write 

<x| = X>| e- 1 '* (9.21) 

V 

where the sum is over all the allowed values of the momentum p. Putting x — 0 

v 

in this equation, we find 

<0| = £«P| • (9.22) 

V 

Thus, the dual Fock vacuum (0| is a combination of all dual states with definite 
momentum (this is of course a natural consequence of (9.9) and (9.10)). Let us 
now compute the following expectation value 

«pi n m*/) io>. 

7=1 


Using (9.14) and (9.19), we obtain 

N 

«pi n v -<*'> i°>= n (ZJ - zjY 6(p-N yfi) (9.23) 

7=1 I<J 

28 We recall that in string theory x and p are respectively the position and the momentum 
of the center of mass of the string. 
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where the 5-function arises from momentum conservation. If we make use of (9.22), 
our result can be written also as follows 

N 

(oi n v ^) i°)= n (** - z jy ■ (9.24) 

i=i i<j 


Notice that momentum conservation selects the momentum state N\ from 
the dual vacuum (0|. Thus, the vacuum expectation value of a product of vertex 
operators is precisely the statistical prefactor which characterizes anyonic wave- 
functions. 

We can proceed even further and use this operator formalism to reproduce 
the entire wavefunctions for anyons in the lowest Landau level, and not just their 
statistical prefactors. In the dual Fock space of the Q field, let us consider the 
state 

« a ‘>i=<°i 0 ?) A W’• • - ■ <9 - 25) 

where {A&} is an ordered set of non-negative integers, and let us denote by d its 
level, which is defined by 

d = Y^k\ k . (9.26) 

k 

A straightforward computation then shows that 


N 

«A fc }i n^)i°>=< } ( 21 ,..Mn^-^r, (9.27) 

1=1 I<J 

where -®{li}( z i, * * * z n) is the following symmetric homogeneous polynomial of 
degree d 



Apart from the exponential factor that has been removed, (9.27) is precisely the 
generic wavefunction for anyons in the lowest Landau level with angular momen¬ 
tum 

■j N(N — 1) J 

J = v—^ - L +d . 

Therefore, the insertion of the oscillators a& in the correlation function (9.24) is 

seen to increase the angular momentum of the base state Yl ( z i ~ z jY * I n fact, 

I<J 

like in string theory, ({Afc}| is a higher angular momentum state. 

It is interesting to consider some applications of this formalism to describe 
the wavefunctions of the fractional quantum Hall effect that we discussed in the 
previous chapter. In particular, it is immediate to realize that the relevant part of 
the Laughlin wavefunction ipm given in (8.9) can be written as 
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N 

V’m = (0| nte) 1°) > (9-29) 

7=1 

with m being an odd integer. Thus, the vertex operator formalism provides a 
naturally factorized expression for the Laughlin wavefunction in which the elec¬ 
tron located at zj is represented by the fermionic vertex operator Vm(zi)* Similar 
considerations apply also to the wavefunction for quasi-hole excitations, whose 
relevant part is (cf (8.17)) 



N N 

= J[( Za ~ zi ) T[( z i- z j) m » 

7=1 I<J 


(9.30) 


where z a is the quasi-hole position. In fact, it is not difficult to see that 'ip£ l Za can 
be written as follows 



N 

= <0| Vl(2: a ) n ^/) |0> 

m 


7=1 


(9.31) 


This expression identifies the anyonic vertex operator Vj_ (z a ) as the one describing 

m 

a quasi-hole excitation at the point z a . Notice that using (9.15), it is immediate to 
see that the statistics of the quasi-hole is v = 1/m, in agreement with the Berry 
phase derivation which we described in Chapter 8 (cf (8.50)). 

The situation is more complicated for the quasi-particle wavefunctions which 
are not holomorphic even after the removal of the exponential factor (cf (8.18)). 
Therefore, to give a representation of these wavefunctions in terms of vertex op¬ 
erators it is necessary to enlarge the previous discussion and include also the 
anti-holomorphic modes. This procedure is well-known in string theory (Virasoro 
1969; Shapiro 1970) where it is used to describe the closed string (Green et aL 
1987). In particular, one introduces a new set of operators, dj, d J 5 x and p, such 
that 



(9.32) 


These new operators commute with those we have previously introduced. Then, 
one defines a new bosonic field Q(z) which is anti-holomorphic and, correspond¬ 
ingly, a new vertex operator 


V v {z) = : : . (9.33) 

Assuming that the states |0) and (0| are Fock vacua also for the new oscillators - 
«.e. that (9.6) and (9.10) hold also for a], d J5 x and p - it is easy to derive that 

N N N 

<oi n v -( z i) ■■ n : i°>=n ^ • 

7=1 7=1 I<J 


(9.34) 
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The normal ordering on the product of the anti-holomorphic vertex operators 
reduces to 1 their contribution to the vacuum expectation value. Therefore, their 
use in (9.34) looks like a cumbersome and unnecessary complication. However, 
the doubling of the operators and (9.34) allow us to represent the quasi-particle 
wavefunctions of the fractional quantum Hall effect (Fubini and Liitken 1991). To 
see this, let us put v — m and insert in (9.34) a vertex operator Vi_(z a ); after 

m 

some simple algebra we get 

N N N N 

<0|Vj_(2 a ) J] Vvi*!) ■ n *<*/> •' = II(^r • (9.35) 

m 7=1 7=1 7=1 I<J 

The right hand side of (9.35) is weakly equivalent to the reduced quasi-particle 

a. _ 

wavefunction ^»" Zo given in (8.18). This means that the two expressions are not 
necessarily identical but yield the same matrix elements which are expressed as 
Bargmann-Fock inner products (see (8.36)). Indeed, as we discussed in Chapter 8, 
inside these inner products one can safely perform the replacements 

" Sl • 

so that (9.35) when projected onto the lowest Landau level, becomes equivalent 

to V'm*" (Girvin and Jach 1983). With this interpretation, the vertex operator 
Vj_(4) can be regarded as creating a quasi-particle at the point z a , and upon 

m 

using the obvious extension of (9.15), it is immediate to realize that the quasi¬ 
particle statistics is v = 1/m in agreement with the Berry phase calculation of the 
previous chapter. Furthermore, since Vj (z a ) represents a quasi-hole and V\_(z a ) 

m m 

represents a quasi-particle, the particle-hole duality (Laughlin 1983) is manifest in 
this vertex operator formalism suggested by string theory and CFT. 

We conclude this chapter with a few comments. In analogy with string theory 
where the Fubini-Veneziano vertex operators can be used to compute multi-loop 
scattering amplitudes on Riemann surfaces of higher genus (see for instance (Di 
Vecchia et al. 1989)), here too one can think of using the vertex operators (9.12) 
and (9.33) to compute the anyon wavefunctions on surfaces of non-trivial topology. 
Indeed, in (Cristofano et al. 1991a,b,c) this has been done for the torus, and 
the genus one generalization of the Laughlin wavefunctions (Haldane and Rezayi 
1985) - as well as their degeneracy - have been obtained in the operator formalism. 
Preliminary results for higher genus Riemann surfaces are also available (Maharana 
and Panda 1991). 

The discussion we have presented here is limited to the case in which all anyons 
are in the lowest Landau level. Despite the extensive literature on the subject 
(Fubini 1991; Stone 1991a; Fubini and Liitken 1991; Cristofano et al. 1991a,b,c; 
Balatsky 1991; Dunne et al. 1991a,b; Balatsky and Stone 1991; Nagao 1992a,b; 
Ting and Lai 1992), it is not known at the moment how to generalize the CFT 
construction with vertex operators to the case in which some particles are in 
excited states. We think that this is a necessary step towards a full understanding 
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of the analogy between the conformal blocks of CFT and the anyon wavefunctions. 
In (Moore and Read 1991) this analogy has been studied in some more detail, 
and the ultimate reason for its validity has been indicated in the relation of the 
two-dimensional CFT on one hand, and of the anyon theory on the other, with 
the Chern-Simons field theory. In fact gauge theories in (2+1) dimensions with 
Chern-Simons action, are topological field theories which reproduce the braiding 
and fusing properties of two-dimensional CFT (Witten 1989; Moore and Seiberg 
1989). On the other hand, Chern-Simons gauge fields are intimately connected to 
anyons as we discussed at length in Chapters 3 and 5. This circle of ideas needs 
to be analyzed and studied in more detail, and perhaps the best place to do so 
is provided by the gapless edge excitations of the fractional quantum Hall effect 
(Halperin 1981) where non-abelian generalization of fractional statistics and Kac- 
Moody algebras seem to play a crucial role (Wen 1990b; Stone 1991b; Frohlich 
and Kerler 1991). Finally, it would be quite interesting to understand the role and 
the meaning of the Virasoro algebra (9.1) in the context of anyon theories, and to 
clarify the significance of its non-linear extensions (w-algebras) which have been 
recently identified in quantum Hall systems (Cappelli et ah 1992). 

These issues are still under current active investigations, and I hope that the 
material and the discussions presented in this chapter as well as in the rest of this 
book, might be of some help in these developments. 
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